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Photoelectric observations are reported for 187 star clusters, 
mostly globular, in the Galaxy, Magellanic Clouds and the M31 
group of galaxies. All were observed in photographic and visual 
light, P and V, and 117 in the infrared, 7. Globular clusters in 

he Galaxy and Clouds were measured through a series of aper- 

ures up to 25’ diameter, to obtain total magnitudes and diameters 
containing 0.9 the total light. For comparison of integrated colors, 
28 galactic open clusters were observed with apertures large 
enough to include most of the cluster members. Some space- 
reddened F- and G-type supergiants having six-color photometry 
were measured to obtain total/selective absorption ratios of 
Ay = (2.90.2) Eyp-v)= (2.31.4) E(y—z. These ratios, with color 
excesses estimated in two ways from spectral types, were used to 
compute corrected total magnitudes, distance moduli and linear 
diameters for the galactic globular clusters. The principal results 
te: (1) The galactic globular and open clusters are generally well 
Separated in the plot of (P-V) vs (V-Z); (2) Linear diameters 
with 0.9 total light range from 20 to 50 parsecs as My ranges 


from —6.8 to —9.6, but the scatter is so large that the correlation 
is not strong; (3) Galactic globulars appear to be systematically 
bluer than M31 clusters, by about 0.2 mag. in (P-V), notwith- 
standing uncertain allowance for space reddening; (4) M31 clusters 
well outside the main spiral structure have an intrinsic color range 
of 0.4 mag.; (5) Except for a few relatively blue objects apparently 
like some in M33, the M31 clusters seen over the spiral have 
colors in the range from (P-V)=+0.50 to +1.94 mag., with the 
reddest being faintest; (6) From eight M31 clusters that are 
brightest and reddest in the V vs (P-V) plot, it was found that 
Ay/E,pe-y)=2.50+0.14, which is not regarded as significant of 
different absorbing matter in M31 than in the Galaxy; (7) Com- 
parison of magnitude-frequency histograms gave estimates of 
distance moduli ranging from 23.5 to 24.0 for M31 and 19 for the 
Magellanic Clouds, with all values uncertain by 0.5 mag.; (8) The 
galactic center distance was estimated in two different ways at 
12.5 and 12.0 kpc, with an uncertainty of about 1.5 kpc, on the 
assumption that for RR Lyrae variables Mp=0.0. 


I. INTRODUCTION 
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HERE is a large literature relating to globular 
cluster photometric characteristics, such as mag- 
Faces, colors, and diameters. Visual, photographic, 
and photoelectric methods have been used, but in sum 
here still is not available a comprehensive and homo- 
eous set of data based on modern standards and 
procedures. Thus, the previous most extensive programs 
for the measurement of total magnitudes are those 
carried out by W. H. Christie (1940), who measured 
by the schraffierkasette method the photographic mag- 
des of 68 galactic globular clusters, and by Seyfert 
Nassau (1945), who measured on in-focus Schmidt 
escope plates the photographic magnitudes of 212 
ebulous objects” in and near the Andromeda nebula 
(M31). Also, the most comprehensive investigation of 
their colors is that by Stebbins and Whitford (1936), 
ho measured photoelectrically 69 galactic globular 


* Lick Observatory Bulletin, No. 569. A preliminary summary 
presented at the NSF-sponsored symposium, “The differences 
lobular clusters,” Am. Astron. Soc. 103rd meeting, 
st 31, 1959, Toronto, Canada; Astron. J. 64, 428 (1959). 
_TNow at Kitt Peak National Observatory, Tucson, Arizona. 


clusters. Likewise, the literature (Shapley and Sayer 
1935, Mowbray 1946, others in H. B. Sawyer 1947) 
contains numerous diameters obtained by visual in- 
spection, microphotometer tracing, and iris photometry 
of photographic plates. So far as we are aware, however, 
no extensive series of globular cluster diameter deter- 
minations by photoelectric techniques has been pub- 
lished. 

The present paper gives the results from photoelectric 
photometry of star clusters, mostly globular, in the 
Galaxy, Magellanic Clouds, Fornax dwarf system, 
and the M31 group. A total of 187 clusters was observed 
in two colors, 117 of these in three colors, and 77 
through a series of graded apertures. For each cluster 
so observed, the measurements were made with an 
aperture size or range sufficient to yield a reasonably 
close approach to integrated color, total magnitude and 
a diameter within which 0.9 of the total light is con- 
tained. In most cases, only a small extrapolation was 
necessary in order to obtain the total magnitudes and 
the specified diameters. These were, of course, limited 
chiefly by the interfering light from the surrounding 
field and sky. The diameter results for 10 Magellanic 
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TABLE I. Filter and cell combinations. 
Filters (balsam-cemented) 
Observations Cell Blue Yellow Infrared 
Small-field photometer 
3-color, with separate 1P21 4 0.4C30608+1 mm BG12 2 mm GG 14+ glass 
infrared 
CE25A/B, e W88A>-+-cover glasses 


3-color, simultaneous Lallemand, II 
coated 
1-color, separate infrared, Lallemand, I 


for M31 clusters 


Large-field photometer 
2mm GG13+1 mm BG12 


Dumont 6291, No. 1 


CE23A/B, a 
Lallemand, I 


3-color, with separate 
infrared 


3-color, simultaneous 
+2 mm GGi3 


2mm BG 2341 mm BGi2 2mm Get +2 mm BG18 
+2 mm WGI cryolite 


2 mm BG23+1 mm BGi2 


2mm BG21+2 mm RG8) | 

| 

2mm RG2 commercially 
coated 


2mm GG14 
W88A-+t cover glasses 


2mm GG1442 mm BG18 C3966+W89-+cover glass 
not cemented | 


® Corning Glass number, 40% standard thickness. 


b W refers to number of Wratten gelatine filter, furnished by Eastman Kodak Company. 


Cloud clusters are reported here for the first time with 
the kind permission of S. C. B. Gascoigne, who, with 
Kron (1952), had previously published colors and total 
magnitudes obtained with the Mt. Stromlo Observa- 
tory’s Reynolds reflector and a series of four apertures. 
Mainly for comparison purposes, color data also are 
given for 28 galactic open clusters. 


Il. INSTRUMENTATION 


All extragalactic clusters were measured in the blue 
and yellow with a 1P21 multiplier, and at Mount 
Hamilton the Crossley reflector was used exclusively 
for the measurements. Clusters in the Galaxy were 
measured with a stellar photometer, but the larger ones 
were observed in addition by means of a special, large- 
field photometer (Kron, Greeby, and Willson 1956). 
In the early part of the program, 35 galactic globular 
clusters were measured in the blue and yellow with a 
1P21 multiplier phototube and in the infrared with a 
CE-25 photocell on the stellar photometer, or with a 
DuMont 6291 multiplier phototube (blue-yellow) and 
a CE-23 photocell (infrared) on the large-field photom- 
eter. In December, 1955, two superior infrared-sensitive 
receivers were put into operation: one, a large cathode 
7-stage multiplier used with the large-field photometer, 
and the other a 12-stage multiplier with effective 
cathode size restricted by means of an electronic lens, 
used with the stellar photometer (Kron 1958). Both of 
these multipliers were generously furnished by Pro- 
fessor A. Lallemand, who developed them in his labora- 
tory at the Paris Observatory. The new multipliers, 
employed with a revised filter system, made possible 
the measurement of all three colors, P, V, and J, in one 
operation, which was impossible with the infrared- 
sensitive photocells owing to insufficient sensitivity in 
the blue spectral region. All photometric receivers, with 


the exception of the CE-23 receiver, were refrigerated 
with dry ice. 

In accordance with the practice of Johnson and 
Morgan (1951), the ultraviolet transmission of all blue 
filters was limited by a 2-mm thick piece of Schott 
GG13 glass; none of our blue filters had a red or infrared 
leak. The photometric system was stabilized by careful 
matching of filters with receivers, so that it was in 
practice unaffected by the changes in receivers made 
during the program. Thus the filter and receiver com- 
bination was the important unit, rather than only the 
filter; Table I lists the various filter-cell combinations. 
The stellar photometer is equipped with focal plane 
apertures graded in diameter from 4” to 5/5. The large- 
field photometer has apertures ranging from 3°7 to 
24/8; there is an overlap with the stellar photometer of 
two apertures. The large-field photometer is also 
equipped with graded occulting disk fittings, which 
permit transforming the circular-hole apertures into 
annular apertures. Photometer sensitivity was moni- 
tored with a radium-activated light source during all of 
the observing. 


Ill. THE PHOTOMETRIC SYSTEM AND STANDARDS 


A large part of the published photometry of star 
clusters is on the International system; we therefore 
adopted, as our photometric system in the blue and 
yellow, a photoelectric approximation to International. 
The only measurable systematic deviations resulted 
from use of an ultraviolet-free blue color. Although this 
circumstance causes some departure from International 
photographic, we feel that the use of an ultraviolet-free 
blue color is justified for reasons given by H. L. Johnson 
(1952). Our standard magnitude system is based upon 
observations of North Polar Sequence stars, and it is a 
close approximation to the International photovisual 
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TABLE II. Photometric standards. 
_—== 
| North Polar Sequence 
Star Vi (P-V) (V-1) Sp. T Star V (P-V) (V-D Sp. T 
[PS 1 4™39 —0"06 —0™21 A2 C8 
me 2 5.26 (0), —0.30 B9 
3 5.61 +0.12 +0.06 — FOIIT s COS i ee a Ore GOV, 
4 5.83 22012 0.00 A3p B 7.30 +0.10 +0.05 A7p 
5 6.47 —0.01 ai AS © 6.98 +0. 84 +0. 82 KOU 
D 6.85 +1.20 Sei) 15 K31II 
6 7.13 +0.09 : A3 E 8.21 +0.21 +0.21 FOI 
7 7.53 =O, ih SO Hi BO 
Pi 8.10 +0.29 +0.23 F4III F 7.28 +1.07 +1.01 K2III 
10 9.07 +0.16 oo 9 
13 10.32 +0.22 A7 D10 
€ 7.80 =0107 =). iv! B9 
16 11.23 +0.42 ey #5 E 7.71 +0.38 +0.40 FS 
19 12.21 +0.47 ez = F 8.27 —0.08 —0.18 B9 
2r 6.47 +1.46 +2.14 MI:III H 7 12 a7 eat KO 
3r 7.54 +1.30 +1.32 K2iil Me 7 {5- 41.28 A303 Mb 
4r 8.26 +0.97 +0.89 GOIII pers 
BE Piccme seta 2 e197)’ K-00 ° aes sate ot TA Ms2 
6r 9.27 dei Ml +1.24 GOI C12 
7r 9.88 +1.08 out G8IV 
8r 10.46 +0.97 Aa 6.78 —0.10 —0.24 AO 
12r 12.61 Seq es GC 7.42 +0.04 —0.04 A3 
D 6.65 SS +1.10 KOIII 
2s 6.29 +0.25 +0.17 F21I1 E 7.17 +0.88 +0.88 KOMI 
3s 6.36 +0.30 +0.32 F4V F 8.27 +0.22 +0.16 AOp 
4s 9.83 +0.48 +0.39 FSIIL 
6s 10.68 +0.71 +0.66 25 G 8.36 +0.21 2b FOp 
I 8.13 ATL TY? = teat K3III 
C6 L 9.20 +0.40 =e FSIII 
M 8.63 aio SEARO K3III 
A 6.93 —0.09 So Onliy A3p N 9.74 +0.42 ae KSIII 
» +0.43 +0.45 GOV 
.06 +0.33 +0.35 F5V 
2 a dana aes ae 6) 9.21 +0.99 +1.05 KOI 
K 7.53 +1.54 +1.93 MQIII 
M 8.41 ily +1.16 K31 


* Mag. var: Kron, White, and Gascoigne (1953). 


system as defined by Seares and Joyner (1945). The old 
Tnternational photovisual magnitudes are, in fact, quite 
Boa, and the use of 24 stars permitted a very satis- 
factory photoelectric copy of the system. The 24 NPS 
Stars listed in Table II were chosen to include a generous 
assortment of bright stars easy to find, with a good 
range in color. In the selection an effort was made to 
obtain nonvariable red stars and those of intermediate 
color, which are rather rare but important in establish- 
ing linearity of color and magnitude transformations. 
In addition, these 24 stars include the nine observed by 
Stebbins, Whitford, and Johnson (1950), although one 
of the nine is NPS 2r, a red star of doubtful light 
constancy. We designate our photoelectric copy of 
International photovisual by V, as originally suggested 
by Kron and J. L. Smith (1951). With (P-V) as abscissa, 
Fig. 1, upper panel, shows by filled circles the deviations, 
pv—V, between our magnitudes V and the International 
magnitudes, pv as given by Seares and Joyner (1945); 
those with crosses represent the nine NPS stars observed 
by Stebbins, Whitford, and Johnson. The two most 
discrepant points (in parentheses) correspond to NPS 
2r and 12r; if these are disregarded, a horizontal straight 


line, indicative of no color equation, may reasonably be 
drawn through the remaining points. Figure 1, lower 
panel, shows with open symbols the differences, AV, 
between our V magnitudes, and those of Johnson and 
Morgan (circles) and of Eggen (triangles). This plot 
indicates that magnitudes on our system are linearly 
related to the other two, within the range of (P-V) 
from about —0.4 to +1.0 mag., according to the fol- 
lowing formulae: 


Vim=Vxm—0.075(P-V)—0™027, 
Ve=Vxm—0.115(P-V)—0™026; 


where V xm refers to the magnitudes of this paper, Vim 
to those of Johnson and Morgan, and Vx to those of 
Eggen. There is some suggestion that redder objects 
may depart systematically from a linear transformation. 

We designate our blue color index by (P-V) but with 
the difference (from Kron and Smith 1951) that here 
the P magnitude is on an ultraviolet-free basis. For the 
zero point and scale of (P-V), we use the same approxi- 
mation to the International CI employed by Eggen 
(1955). With this definition, (P-V) can be reduced to 
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Johnson and Morgan’s (B-V) system, with a precision 
of +0.03 mag., simply by adding 0.10 mag. to (P-V) 
for the range (P-V)=—0.4 to +1.0 mag. If a more 
exact conversion is required over this range, and if the 
range is extended to include redder objects, then the 
following two linear relationships may be used: 


(B-V)=+0™10+0.96(P-V), 

for (P-V)= —0.4 to +1.0 mag., 
0™00+1.06(P-V), 

for (P-V)=+1.0 to +1.5 mag. 


(B-V)= 


Our decision to use a third color, 7, which gives a 
second color index, (V-Z), was based on the hope of 
deriving color excesses for clusters according to the 
method devised by W. Becker (1938). Also, the success 
attained with the use of this method on stars by 
Johnson and Morgan (1953), who used a three-color 
system employing the ultraviolet, at first led us to 
make a number of ultraviolet observations of globular 
clusters. As may be expected from their color indices, 
however, many globular clusters are too faint in the 
ultraviolet for satisfactory large-aperture photometry. 
The situation is especially bad for the interesting and 
greatly reddened clusters at low declinations around the 
galactic center. Although even the heavily reddened, 
high-luminosity clusters actually have enough ultra- 
violet light to detect, we found in practice that the 
amount of skylight admitted through the large aper- 
tures reduced the ultraviolet cluster-sky ratio. too much 
for satisfactory measurement. Likewise, with the 
Crossley reflector, we were unable to make any meas- 
urements of acceptable precision in the ultraviolet for 
the clusters in M31. On the other hand, experiments 
with the infrared indicated that the practical advan- 
tages in measuring in this spectral region for a third 
color more than made up for the large loss in quantum 
efficiency, which resulted from abandoning the Sb-Cs 
photosurface of the 1P21 in favor of he) infrared- 
sensitive CsO— Ag photosurface. 

The infrared index (V-Z) used here is based on an 
infrared photometric system like the infrared part of 
the R, I system by Kron, White, and Gascoigne (1953). 
Thus if infrared magnitudes J are computed from the 
V magnitudes and (V-J) color indices of the present 
paper, they will be on the same magnitude system as 
the infrared magnitudes J computed from (R-J) color 
indices. 

All of the photometric standards used in this work 
are listed in Table II, which contains more than the 
data for the NPS standards. The awkwardness of direct 
work on the North Polar Sequence, the primary 
standard, persuaded us to establish secondary standards 
in four Hagvard C and D regions. It should be noted 
that I magnitudes derived from these standards differ 
slightly from those obtained from the earlier ones 
(Kron, White, and Gascoigne 1953), because the 
present data in Table II include some additional obser- 
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_supergiants for which recent six-color data were avail- 


vations. In establishing these working standards, many, 
transfers were made to the pole, as well as many cros ‘ 
transfers among the Harvard C and D regions. The 
spectral types with luminosity classes for NPS stars| 
are by Keenan (1940), types for the C region stars are) 
by Nassau and MacRae (1955), and for D10 from the: 
Henry Draper Catalogue. q 

In addition to establishing secondary standards basetl 
upon the North Polar Sequence, we have measur 
some stars for calibration of our three-color system asi 
means for estimating color excess. Among these are a] 
number of stars of known spectral type and presumed 
small reddening, taken from the listing of Johnson and 
Morgan (1953). In order to evaluate the ratio of total 
to-selective absorption on our three-color system for 
objects of intrinsic color more nearly comparable to) 
clusters, we also observed a number of F- and G-type) 


able. Some of these stars had previously been found to’ 
be much reddened in six-color photometry by Kron 
(1958). All these data are included in Table III, which 
is arranged to indicate how group means were formed! 
to define the main sequence, and the giant and dwart 
branches for late-type stars. Since the observing} 
schedule did not provide time enough to determine’ 
new P, V, I data for a complete calibration, a number 
of colors were obtained by conversion from previou 
work. Thus-some (P-V) colors are from the (B- 
colors by Johnson and Morgan (1953), or from 
(P-V)x colors by Eggen (1955). 


4 
4 

A) 
IV. OBSERVATIONAL TECHNIQUE FOR CLUSTERS 


Before a cluster was observed, a photograph of the 
object and its surroundings was examined. Represen- 
tative areas for the measurement of background light 
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6 #9 .2) Ot 2046 Om cme 
(P-V) 


Fic. 1. Differences between North Polar Sequence star magnitudes 
determined by different observers. 7 
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‘TaBLe III, Photometric data for calibration of three-color system. 
Star Sp. T. (P-V) (V-D) Star Sp. T. (P-V) (V-D) 
Main Sequence Giants 
a Leo B8 V —0™24 —0"32 o Tau G8 III +082 +0273 
134 Tau B9 IV —0. 188 2027 P G8 III 0.88 
# Cet B9 III ae —0.26 Ruut sine yah 
a Peg B9V = (13 = 026 Means G8 (2 0.85 0.76 
; C Hya AO V —0. 13 —0.25 @) of a 
8 Tau KO III +0.92 +0.78 
| 109 Vir AO V =O) 14 =0,49 Tau KO III 0.928 0.78 
| begs a iN. —0.118 —0.16 pe AC ay . i 
my! | y Op —0.07 —0.12 Means KO (2 0.92 0.78 
\ ¢ Aqr A1V 0.11 —~0.19 2 e a 
iy HR 875 A1V —0.028 —0.10 6 Oph K2 II +1.04 +0.96 
io — Ser K2 II 1.118, 1.01 
. | Means AO (10) = 042 =n 24 ae See's ae ci 
;) Means K2 (2 1.08 0.98 
4 B Ser A2IV (noe? 0.12 + 
p B Leo A3 V —0.02">  —0.08 6 Cnc K4 II +1.40 4-1),35 
i \ Gem A3V —0.01 —0.08 y Dra K5 Ill q' Ade!) 148 
| ¢ Lep A3 V —0.00* —0.08 a Tau K5 Il 4+1.443b5 +1.54 
i B Eri A3 IIT +0.01 —0.05 = 
i Means K5 (3) +1.43 +1.46 
A a Oph AS III +0.04 +0.04 
. a Aql A7 IV-V +0.128>  +0.08 ee 
Means A4 (7) +0.02 —0.04 e Eri K2V +0.829> — +0.79 
HR 753A K3V +0.91 +0.91 
eo FO UI 40.22» 40.23 HD 156026 dKS +1.03> ze 1s 
vy Vir FO V +0,25> +0.26 +2° 3312 K7V +1.22 +1.47 
+0° 2989 MO V eis +1.66 
| Means FO (2) +0.24 +0.24 
| —3° 1123 M1V +1.33 +2.05 
| a! Lib F5I1V +0.32 +0.30 +2° 348 M2V +1.35 +2.00 
| x Ori F6 V +0.36 +0.35 +1° 4774 M2 V +1.38 +1.92 
| Sel F6V 40.4085 40.37 +1° 2447 M2V +1.40 +218 
| oPse F7V 40.42 40.41 +15° 4733 M2.5V SKY) +2.00 
‘| B Vir F8 V +0.45 +0.44 aca aca 
| _ ES ar Means M2 (5) +1.36 4-203 
q Means F6 (5) +0.39 +0.37 
| | % ae ! F- and G-type supergiants 
| n Boo +0.502: +0.51 
i) & Ser GO V 40.525> 40.48 HD Sp. Ts, (@V) | V-Dy Let 
t 4362 GOIb +1701 +0798  +0™93 
] Means GO (2) 0.51 +0.50 3092 F6Ib +0.80 +0.94 +0.77 
| ' So ei aati sbi val desley 
| x Cet G5 V +0.59 +0.55 150181 ~G21b \ -F0.897  -10.81 2 5)-10. 55 
| 70 Vir G5 V +0.63 +0.60 165553 F9Ib +0.95 +1.04 +41.03 
r Cet G8 V +0.64 +0.65 
| am 172365 F9Ib +0.69 +0.73  +0.30 
| Means G6 (3) +0.62 +0.60 174104 GOIb +0.66 +0.63 +0.02 
| 187203 GOIb +0.87 -+0.88  +0.66 
B Aql G8 IV +0.77 +0.75 JO4020 = <GOLbae, --es6n yates ee 22s 
o? Eri A KOV +0.74> +0.76 204867 GOIb +0.75 +0.68  +0.27 
107 Psc K1 V +0.74 +0.73 
———— 209750  G2Ib +0.88  +0.75 -+0.59 
Means KO (3) +0.75 +0.75 H19135).  GOilb 1-007 910/99") *-61.02 


—— 


« (P-V) from (B-V) determined by Johnson and Morgan (1953). 
» (P-V) from (P-V)® determined by Eggen (1955). 


were selected, and a preliminary decision on aperture 
sizes was made. For most clusters in the Galaxy, and 
for some extragalactic clusters, two or more areas were 
employed for measurement of background light. An 
inconsistency between the measures of more than one 
background area was always investigated, and care was 
taken to ensure that the photometric results were 
reasonably free from systematic errors due to unusually 
unfavorable background conditions. Because of bright 
Stars in the field, exploratory measures on a few clusters 


} 
) 


r 


were rejected; for some others, large-aperture photom- 
etry was not attempted. After elimination of such cases, 
most of the remaining globular clusters in Helen B. 
Sawyer’s (1947) list, and north of 6=—38°, were 
measured. Clusters in the Magellanic Clouds were ob- 
served in 1951 by Gascoigne and Kron by means of the 
same technique, with Harvard Region C12 as a 
standard (Eggen 1951). 

No attempt was made to measure diameters for the 
clusters in the M31 group, although this may be a 
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good possibility with more powerful telescopic equip- 
ment. All our measurements in the M31 group were 
made through a focal-plane aperture of 1377 diameter, 
and the standard stars were measured through the 
same aperture. Measurements were not made on nights 
of bad seeing, as experience indicated that appreciable 
light was lost from these clusters when stellar images 
appeared to be much larger than 2’ in diameter. Good 
seeing was also required for reliable, visual centering 
of the fainter clusters in the aperture. 

The photometry of M31 clusters through an aperture 
so large as 13"7 required close examination of the field 
around each object. Stars fainter than 17.5 mag. could 
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Tas_e IV. Magnitudes and related diameters of 67 galactic globular clusters. 


not always be seen reliably with the Crossley even| 
during the very best seeing; yet, a star of 17.5 mag. in| 
the same field of a 16.0-mag. cluster could cause the} 
cluster measurement to be wrong by 0.25 mag. In order} 
to minimize errors of this type, fields around the M31} 
clusters were inspected on enlarged prints, and on dia- } 
positives made from 100-inch reflector plates, kindly} 
furnished by Dr. W. Baade. The comparison sky fields 
were selected as close.as possible to the cluster, usually } 
only about-1.5 aperture diameters away, in order to} 
minimize errors from background brightness gradients } 
in the spiral. y 

M31 is currently estimated as about 60 times more} 


Aperture 0/33 0'49 0/75 esa 4 1/65 2/44 3/61 5138 6/95 
Diameter 8.29 9.95 LLG: 13.4 16.6 20.0 24.8 
NGC M Cone. } 
V magnitudes* through above apertures* Ve n do.9 n Class 
2419 127908 12726 11767 11727 10791 10765» .-- tes <10™9 5 cee te a 
4147 11.97 11.21 10.90 10:66’, 10.46, 10740) «10732 tee 10.26 2 Re i os) (a,b) 
4590 68 tee 10.13 9.42 9.01 8.81 8.71 8.46 8.35 8.31 3 6.4 3 ba 
5024 53 10.53 9.75 9.19 8.69 8.33 8.08 8.00 7.82 ley 2 8.3 1 a,b | 
VAL Ca RED, _— 7.78 a nese ace 0 
5053 ane or cae oe fe ae 10.64 9.9 3 oe 
10.36 10.14 10.05 /. ve ane wes . 
5272 3 9.44 8.72 8.00 da 5 (fea) 6.86 6.64 6.58 (eit 8 ORS me ba 
653) (65494 16c45. 6 43 01S iy 
5466 ao tee tee ve tee 9.74 9.56 9.2 3 Sf 
9.39 ewer pier ae sate noe 
5634 Bets ee tere Lot ee OURAby eee se Rens <9.8 2 ote 
5694 11:52; 11-34 10:95. 10560 105445) 105405 210527; see tes 10.25 2 Bs 2 
5824 10.19 — 10.04 9.78 9.52 9.36 9.25 9.14 9.07 tee 9.04 4 See) p. 
5897 vee tee 12.30 21.36 8 10x61) 10305 9.40 8.95 8.80 (8:4) +2) (sire 
5904 5 9.43 8.55 7.70 7.26 6.87 6.59 6.29 6.17 5-93 4+ 10.7 3 
5986 10.47 9.72 ON Ga oho 8.25 8.03 Teste ee Hos <2 Oye | 
6093 80 9.28 8.82 8.47 8.14 7.91 7.68 Peot 7.45 (eel) 8) ROS 2) 
6121 4 tee tee 9.29°' 8.74 7.83 7.34 6.89 6.65 5.901 2) 225Gumeee 
6.45 6.57 6.18 6.17 . 6.08 tee 
6144 oe 11.46 10.70 9.92 9.62 9.29 9.51 9.00; 4s. "eaten 
Chit OF 127169 11355 10560 9.98 9.51 9.19 8.75 8.53 8.8 ~3  pIQESeae2 
8.45° 8.40 tee 8.31 ee tee tee 
6205 13 10506) 9220) e840) 7.69 7.05 6.58 PE = oi) SP Si aan eS 3 
600 =) OOD NY 5 0008 S506 ee OZ t renee ve 
6218 12 11.49 10.66 ae MOS RSH 7.85 7.36 eves} 6:72 135 72135 aeee 
Scho ; wate 6.89 shake Sunt areie y 
6229 11749 10582) 10551 OO ons t 9.65 9.57 9.48 tee 0 Ade 2 Fe oO 2 
6254 10 ee 10.98 9.95 9.20 8.50 7.97 7.48 7.22 7.05 6645 (32 "1632 ee 
: tee 6.86 6.79 6.73 tee 6.70 oe 
6266 62 tee ao sae UAW Hoo 7.06 6.89 ae OL 3 SiS 
6273 19 10.15 9530" 8507, 10) 7.64 7.34 Gals 7.06 6.88 4 Ons 
6284 10.84 10.53 10.14 9.81 9.53 9.39 9,22 9.05 BeOS, 2 er 2 ca 
6287 tee 1100 salle 25 Od Oma el Onni7 9.90 9.74 9.48 9.48 3 SeSae b 
6293 11.00 10.50 10.17 9.53 9.14 8.87 8.69 8.59 i432 On 2 irae (om 
6304 1210) 142110534) Ons 9.46 9.14 8.88 8.68 Ses oS SOR b 
6316 see NO OOO! MOSES 9.64 9.40 8.90% Ora i soins b 
Se UPABTAL GIAERE isis LSS) RY TIO). tee 10.91 3 ASO 13) b 
aN 9 10.41 9.77 9.25 8.76 8.43 8.20 7.99 7.92 7.79 4 oy 8) bg 
92 9.06 8.45 8.00 7.56 7.24 7.00 6.81 6.71 6.53 2) 1253 eae b 
aD 6.64 Ae 6.53 deh nas srs 
poe 11925 1837 104 O70 O43 en On22 9.86 10.10 2 4.7 2 b 
tee 12,09 11539. 10.90' 10.62) 10,28 9.98 9.43 - 9.80 2 pil 2 b 
aie 10.84 10.04 9.52 9.14 8.87 8.68 8.50 8.40 3 Guay 2 b 
eo , ina ae wee Sn 9.85 at ace ae <9.8 1 ae On 
4 11.67 10.68 9.84 9.08 8.57 8.16 7.89 7.78 756. Vase OeS ames a 
bas TAO sed 7 60 BSG ae She 


| servation of M31 clusters through an aperture of 1377 
is equivalent to observing globular clusters near the 
galactic center through an aperture of approximately 
as many minutes. Since the diameters containing 0.9 of 
the light of the galactic center objects do not exceed 
10’, the total magnitudes for the M31 clusters do not 
need to be corrected for aperture diameter before being 
compared with those for the galactic center globular 
clusters. 

At an early stage in the photometry of the galactic 
globular clusters, measurements in each of the three 
_colors were made through several apertures for 48 

clusters. Since these data indicated that there was no 
‘appreciable dependence of diameter on integrated 

color, the program was completed by making color 
measurements through only one aperture that admitted 
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TABLE 1V.—Continued. 
| Aperture 0/33 0/49 0/75 TAZ 1/65 2144 ~ - 3/61 5/38 6/95 
Diameter &.29 9,95 TO} 13.4 16.6 20.0 24.8 
; Conc. 
M V magnitudes* through above apertures® V; nN dng nm Class 
tee leon 90) 13233 12.71 11,83» 11.275 10.61% <a, 3 tee tee tee 
iva 10.85" 10.30 10.14 10.01 9.82 9.61 9. 33> 9.40 2 Race) 2 (c) 
12.26 11.90 11.49 11.16 10.94 10.69 10.42 3 Hell 2 c 
10.24 9.70 9.32 9.10 8.89 8.76 8.53 8.63 3 4.1 2 b 
10.91 10.42 9.98 9.78 9.63 OEOe) tee tee 9.54 3 3.0 2 a 
see tee cee Ziel Oe Soe OLAZ SOS 9.96 9.67 5 Sissy 8} b 
1ivtad 10.58 9.77 9.31 8.91 Soo 4 o@RIe ado (8.3) 3 (S24) 2) (b) 
11.59 10.91 10.09 9.43 8.92 8.63 8.43 On22 4 See 2 a 
TAR tel AOS? 9.84 9.47 9.20 8.96 8.63» 8.66 4) 6.6 1 b 
tee tee 9.20 8.84 8.58 8.47 8.40 tee 8.35 2 4,2 2 a 
28 see tee =a 8.00 7.68 7.40 7.18 7.10 6.95 3 9.1 3 (b) 
69 10.20 9.46 8.85 8.44 8.16 7.96 7.82 tee 7.67 2 6.8 2 b 
is 106A 10504 9.63 9.44 9.32 9.22 9.17 2 4.3 2 a 
10.49 10.13 9.77 9.36 9.18 9.05 S572" ecb 8.92 1 4.2 1 b 
22 10.45 9.56 8.67 7.84 (alls) 6.50 6.06 tee 5.09 2 26.2 2 b 
DoS oo GES: 5.43 Dost Dos) DB 
70 10.10 95612 59503 8.85 8.56 8.35 8.20 tee 8.17 2 Baal 2 (b) 
11.82 11.01 10.26 9.57 9.16 8.85 8.50 8.38 8.13 4 1253 3 b 
54 9.04 8.67 8.38 8.10 7.93 7.84 Uo oon LAGS AES ae? Cc 
10.82 9.89 9.26 8.61 8.13 7.76 ToS PAL 2 Use 2 b 
HANS ~ SUES TORWeE © IOE27 9.88 9.57 9.35 9.13 3 8.9 2 b 
56 tee 10.81 9.94 9.45 9.05 8.71 9.52 8.21 3 10.1 2 b 
55 : 11.23 10.16 9.14 8.47 ati 7.11 6.30 2 Pil il 2 a 
58 6.66 6.60 6.52 na 5.46 Siete 
| 71 12509" 11-03. ~~ 10.53 9.82 9.28 8.82 8.64 8.46 8.27 2 HO. 2 (a) 
| di08 8.39 SO HAO coe soe 309 
6864 75 9.89 9.65 9.32 9.06 8.92 8.78 8.70 8.60 2 4.9 2 c 
6934 10.74 10.18 9.73 9.47 9.28 9.10» ---. 9.12 3 AC) 2 a 
)\6981 72 122 tod 10079" 10. 23 9.86 9.66 9.50 tee 9.44 2 6.4 2 a 
7006 HPO som edie 20e 1 045 10.37 > 10.76 10.72 tee 10.72 2 3.0 2 b 
7078 15 8.53 8.00 7.56 7.23 6.96 6.80 6.63 6.54 6.36 7 9.4 2 c 
6.50 6.46 560 ee oe ea Res 
| 
“70892 Priesist 7.99 7.46 7.08 6.83 6.66 6.61 G55 2a OS eed ee 
i Sah 6.57 BAP 6.57 Ba ee Bete 
| 7099 30 9.94 9.42 8.85 8.47 8.12 7.93 7.76 - 7.66 7.60 D 6.8 2 b 
a Magnitudes in italics correspond with the apertures in italics at the top of the same column. 
+ One or more bright field stars in the aperture or nearby. 


about 75% of the total light. There was a resultant 
gain in accuracy because of the improved contrast 
with sky background, as compared with apertures that 
admitted a larger percentage. Finally, all the diameter 
measurements were repeated through a colorless filter 
to gain light, and through six to twelve more apertures 
of different sizes to gain observational weight. 

Magnitude and color observations were corrected 
for atmospheric extinction with mean extinction coef- 
ficients determined from stellar observations with ap- 
propriate cells and filters. Clusters far south were 
observed close to the meridian. 


V. TOTAL MAGNITUDES AND DIAMETERS 


The determination of galactic globular cluster total 
magnitudes and angular diameters involved the photom- 
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Relative Magnitude Difference 


5Logd 


Fic. 2. Light concentration curves used to estimate total mag- 
nitudes and diameters containing 0.9 total light. 


etry of surface brightnesses that decrease gradually into 
the field-star and sky background. Ideally, the obser- 
vations should be pushed to the diameter where the 
photometer is unable to detect any residual cluster 
light in the presence of the interfering field-star and sky 
light. This point of diminishing or no returns depends 
upon such factors as night sky brightness, atmospheric 
extinction, and field-star density. Consequently, cluster 
total magnitudes and diameters are bound to be less 
precise than magnitudes for stars or diameters for 
sharply defined objects. In an attempt to allow for these 
factors, we have used a large range of apertures, several 
colors, and two photometers. The most important 
innovation, however, probably was the design and use 
of one of the photometers with apertures large enough 
to reach, in nearly every large object observed, the place 
where cluster light was almost indistinguishable from 
background light. In a few cases, because of brighter 
field stars, or of cluster low surface brightness, the ob- 
servations could not be made as close to this limit as 
for the large majority of clusters. These cases are noted 
by the superscript letter b in Table IV, which gives the 
individual magnitudes measured through the various 
apertures, and the total magnitudes and diameters 
derived as described in the following. 

For each cluster the magnitudes measured through 
the different apertures were plotted as ordinates against 
5 logd as abscissa, where d is the angular diameter in 
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* 5053 and 5466), too faint (NGC 6426), or too close t¢ 
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minutes of arc. This procedure gave a convenient scal 
on which to compare the cluster “curves of shape} 
regardless of distance. For it is clear that a plot of mag 
nitude for a particular cluster against the paramete| 
5 logd will apply to all other clusters having the sam} 
brightness distribution function, subject only to trans 
lation to allow for differences in linear size, distancg 
and interstellar absorption. With the intercompariso) 
of the cluster curves of shape thus facilitated, man 
superpositions of such curves on transparent papé 
indicated that a fair proportion, about 80%, of thi 
observed clusters could be placed in three categories 0 
central concentration. These categories and typical exam 
ples are: a (low), NGC 6205; b(medium), NGC 5272 
and c (high), NGC 7078. The exceptional 20%, fo} 
which concentration classes are not given or ar 
assumed (in parentheses), are very open objects (NG 


bright field stars (NGC 5634). Next, mean curves wer 
drawn through the individual curves grouped for each 
category, and they were extrapolated to a horizonta} 
tangent at their upper ends, as shown in Fig. 2. Finally) 
the mean curves were traced on transparent paper 
superimposed on the individual cluster plots of 
against 5 logd, and used to estimate the total magni 
tudes, V,;, and diameters including 0.9 of the to 
light, do.9, listed in Table IV. | 
As one would expect from the mean curves of Fig. 2. 
the derivation from their upper ends of total magnitudes 
did not involve serious extrapolation, because of the use 
of large apertures. For any particular cluster, the extra- 
polation actually made to obtain the total magnitude; 
V., is given by the difference between V; and the V 
magnitude measured through the largest aperture used 
without interference from bright field stars. If V—V, 
= Am, then the data in Table IV give (Am)=0.13 mag., 
with a range from 0.0 to 0.4 mag; there are only 12 
clusters for which Am>0.2 mag. The case is different, 
however, for the derivation of limiting diameters cor- 
responding to the total magnitudes. Since the curves 
approach parallelism with the 5logd axis, a small 
change in the cluster magnitude is accompanied by a 
much larger change in 5 logd. Although an attempt was 
made to obtain limiting diameters, these proved to be 
so large and uncertain, with extrapolation factors for 
diameters averaging 4.1 and ranging up to nearly 10, 
that clearly some other diameter definition should be 
used. Some trials indicated that a fair compromise 
between interpolation and extrapolation could be 
obtained with a diameter that includes 0.9 of the total 
light. With this specification, the ratio r, between do.9 
and the diameter of the largest aperture used without 
field-star interference, averages to give (r)=1.2. For 
27 clusters r<1.0, and for 33 clusters r>1.0; only four 
diameters have r>2.0, and of these the two largest 
have r=2.3 and 3.4 for NGC 5897 and 6544, respec- 
tively. The first value goes with the largest extrapola- 


P. 
'| tion (Am=0.4 mag.) made to obtain a total magnitude, 
|while the second is a marginal case involving field-star 
‘interference. 


VI. ERRORS IN THE PHOTOMETRY 


Compared to the errors inherent in the cluster data, 
‘those in the selected stellar standards probably are 
Inappreciable, because the point-source character and 
“greater relative brightness practically eliminate sky 
"background as a serious source of systematic error. On 
the other hand, the cluster observations may not be 
so free of systematic error from sky light: the galactic 
‘globulars because of their large size, and the extra- 

alactic clusters because of their faintness. In both 
cluster groups, however, the internal errors of measure- 
“ment seem satisfactorily small, as indicated by the 
average deviation, AD, from the mean of n nights of 
| observation, listed with the data in Tables VI and VII. 
These average deviations infrequently exceed 0.1 mag., 
and many are less than 0.05 mag. Nevertheless, there 
are several potential sources of systematic error that 
“require discussion, separately for the galactic and 
_ extragalactic clusters because somewhat different pro- 
_ cedures were followed in their observation. 


| Galactic Globular Clusters 


For practical reasons involving photometer sensi- 
tivity, an indirect method was used initially to deter- 
mine the (V-Z) colors. This method was to measure the 
I magnitudes with the infrared-sensitive photocell 
‘independently of the V magnitudes, which were meas- 
_ured by necessity with the 1P21. This procedure was 
adopted only after experience had shown good agree- 
ment between (V-Z) colors obtained either differentially 
on the same night, or by differences.of V and J mag- 
‘nitudes measured on different nights. Later in the 
program the Lallemand multiplier phototube became 
available, and then all three magnitudes, P, V, and I 
"were measured successively the same night for a given 
cluster. 

‘e To minimize systematic errors from field stars in the 
large apertures, sky measurements were made for each 
‘cluster on an adjacent field that appeared visually 
“normal in star density. When it was difficult to make a 

“selection, two or more fields were observed, and their 

“intensity was averaged to correct the cluster obser- 
vations. In some cases, when moderately bright field 

tars could not be avoided, they, too, were observed 

d their light was subtracted from that of the cluster; 

ut, when such stars were so bright that they outshone 

the cluster, observations of the latter were not at- 
tempted. We wish to emphasize that the observations 
for the sparsest clusters, like NGC 5053, could be con- 
iderably affected by the background sky light, but 
because there are no other independent measurements 
of similar kind, we can only hope that our data for poor 
‘clusters are of the correct order of magnitude: total 
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Fic. 3. Comparison of infrared color indices measured with extra- 
transparent red filter (above) and regular one (below). 


magnitudes probably reliable to 0.5 mag., and diameters 
to 25%. Lastly, and for the same reason—no other 
comparable photoelectric data—it is impractical to 
evaluate in a detailed numerical manner the possible 
systematic errors in the diameters. For clusters of 
moderate-to-high surface brightness and located in 
uniform star fields, the night-to-night values of the 
diameters usually were accordant to about 5%. Since 
this degree of internal consistency seems satisfactory, 
we therefore conclude that the diameters tabulated in 
Table IV probably are reliable to 10% or better, except 
for the most open clusters, as mentioned in the fore- 


going. 
Galactic Open Clusters 
Compared with representative globular clusters, these 


clusters are of considerably larger apparent size, and of 
much lower mean surface brightness. In addition, they 


TaBLeE V. Comparison of infrared filters. 


Star Regular Light Colors 
or Filter Filter Tables II, VI 
Cluster (v—7z) (v—t) —1™50 (V-I) 
C8A —0715 +0711 +049 
B —0.53 —0.26 +0.05 
(3 +0.15 +0.44 +0.82 
D +0.42 +0.74 +1.15 
F +0.31 +0.61 +1.01 
NGC 5024 +0.10 +0.35 +0.71 
5272 +0.12 +0.41 +0.79 
5904 +0.12 +0.42 +0.80 
6205 +0.09 +0.38 +0.71 
6254 +0.38 +0.70 +1.14 
6341 +0.10 +0.38 +0.73 
6402 +0.78 +1.16 +1.56 
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TABLE Vla. Colors of star clusters in the Galaxy. Globular clusters. 
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NGC M (P=V) Dy (V-D) AD 14 NGC M (P-V) AD x2 (V2) RADI Te 
24198 0™64 5 iS 0™92 7 + 6402 14 ey 4 5 1256 3 4 
41478 0.48 0 2 0.67 1 3 6426 0.85 18 3 1.18 8 3 
45908 68 0.59 6 3 0.81 2 3 6640> 1.81 6 3 eel 0 Zz 
5024* 53 0.50 0 2 0.71 3 2 6453» 1.04: 2 2 1.90: 6 2 
5053* 0.55 14 3 0.76 13 3 6517 1.67 1 3 2EOL 2 2 
52728 3 0.56 4 3 0.79 2 3 6522» 1.02 8 3 1.38 1 SE e 
5466? 0.64 5 3 0.67 2 3 6528» play 3 3 1.56 3 3 
56348 0.56 2 2 0.87: 22S 6539 1.67 5 4 223 8 3 
56948 0.61 0 2 0.87 3 3 6544> LD, 6 3 1.85 5 3 
58248 0.63 8 4 0.94 7 3 6553» 1.39 1 2 Pye fl 4 z 
58978 0.62 0 2 0.86 2 2 6569» a Us| 6 2 1.56 3 Zz 
59048 5 0.63 1 | 0.80 2 f 6624» 1.00 7 fe 1.34 eZ 2 
5986 0.75 2 2 1.12 4 2 6626 28 0.97 2 3 1.28 2 3 FF 
6093 80 0.76 2 4 0.95 4 4 6637» 69 0.87 5 2 1.18 6 2) 
6121 4 0.92 1 2 1.34 2 z, 6638» 0.96 1 2 1.46 1 2 
6144 0.82 5 6 0.98 6 5 6652» 0.78 ses 1 1.02 oes 1 
61718 107 0.97 6 5 1.34 5 3 6656 Hyp 0.86 1 2 1.42 5 2 fi 
6205* 13 0.57 2 6 0.71 3 - 6681 70 0.59 2 2 0.87 0 2% 
62188 12 0.77 1 = 1.02 2 4 6712 1.04 3 4 Lio? 7 4. 
62298 0.64 1 2 0.88 3 3 6715 54 0.70 3 2 1.00 1 2 3 
62548 10 0.79 1 3 1.14 1 2 6723 0.61 1 2 0.92 4 2 
6266» 62 0.95 5 3 1.60 11 3 6760 1.55 2 3 2.12 1 2 
6273» 19 0.90 3 + 122 1 2 6779 56 0.74 4 3 1.04 2 2 
6284> 0.91 1 2 1.20 1 2 6809> 55 0.55 0 2 0.83 5 2 
6287 113 2 3 1.73 1 2 6838 71 0.95 1 2 1.34 9 2 
6293> 0.86 1 2 1.16 2 Z 6864" 75 0.71 0 2 1.09 “ 2 
6304> ates 2 3 1.68 2 2 69348 0.62 1 3 0.84 3 3 
6316» Ae/ 4 3 1.81 1 2 6981 72 0.57 0 2 0.81 5) 2 
6325» 192 2 3 2.14 2 2 7006 0.61 1 Zi 0.84 5 2 
6333> 0.87 3 4 Par 3 3 70788 15 0.59 1 2 0.70 4 2 
6341? 92 0.53 0 2 0.73 3 2 7089* 2 0.57 1 2 0.73 5 2 
6342» 1.20 2 2 1.42 2 2 70998 30 0.48 0 2 0.72 1 2 
6355> 1.41 3 2 1.76 5 2 
6356 1.06 os 3 1.38 1 2 

6401» Lilia 1 IRL} Re a0 1 


* Clusters with |b|> 20°. “ 
b Galactic center group. 


are usually found close to the galactic plane and in rich 
star fields. A single bright field star, if measured along 
with the cluster and not properly compensated by an 
equivalent object in one of the background measure- 
ments, could easily seriously affect the resultant meas- 
urement on the cluster. It must, therefore, be recognized 
that these photometric data for the open clusters prob- 
ably are more subject to systematic error than those for 
the other objects observed in this program. 


Extragalactic Clusters 


For these objects, too, we are convinced that the 
most serious source of systematic error probably is the 
failure of a sample of comparison sky to be accurately 
representative of that for a cluster, particularly in parts 
of systems having steep structural gradients or high- 
intensity levels of light. The error from intrusive field 
stars, however, is expected to be small for clusters 16 
mag. or brighter, because a 19.5-mag. field star would 
produce an error of only 0.05 mag., and such a star 


j 


could hardly have been overlooked in the examination 
of the survey plates for comparison sky fields. Never- 
theless, the results for the fainter M31 clusters may in 
some instances have larger systematic errors. For ex- 
ample, the object B327 was measured as 18.9 mag. in 
the blue, so that a field star of photographic magnitude 
21.5 in a comparison sky field could cause an error of 
0.1 mag. An attempt to minimize this possibility was 
made by observing several sky fields and checking their 
agreement. Also, for clusters such as H89 in fields having 
a steep background-light gradient, the sky brightness 
was measured on several sides of the cluster, and the 
cluster measurement was corrected by linear inter- 
polation. This procedure yielded results for which 
systematic errors from sky background are unlikely to 
be greater than 0.1 to 0.2 mag. 

In addition to the foregoing type of error, there was 
a possibility that the infrared measurements of the 
M31 clusters could be affected by a systematic error 
of instrumental nature. Following some preliminary 
work, it was found desirable to gain more light by use 
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Tas Le VIb. Colors of star clusters in the Galaxy. Open clusters. 


NGC,iC M (P-V) AD (V-I) AD an Aper. NGC,JICM (2-V) AD (V-I) AD xn Aper 
457 0™47 0 0™62 Se 1646 2099 37 0mP45 2 O52 4 4 2478 
559 0.69 291536 8 4 9.95 2682 67 0.68 8 0.71 Me 3} bees 
581 103 0.30 3 0.88 Ph 98) 8.29 AZo 25 0.48 -- 0.81 ---- 1 24.8 
637 0.25 4 0.30 Stet 5.38 6633 0.20 1 0.03 8 2 24.8 
663 0.63 ON rS1 s 10. 2 24.8 6811 OSS 0.80 4 2 1056 
869 h Per 0.29 2030) Soe 24S 6866 Masil il 0.38 SB ES 
884 x Per 0.54 Cele Ss 3 24.8 6940 0.49 1 1.20 4 2 24.8 

1805 0.43 10 0.63 6 4 24.8 7092. 39 —0.14 4 -0.16 15 2 24.8 
Tr2 0.40 OpmeOnS2e 9 12°) °4, 24.8 7209 Ost 4 0.79 8 2 24.8 
1027 0.41 6 0.54 2) 4> 24.8 7243 —Nalsy SG SOe08 Sap 248 
1502 0.38 4 0.50 Deitee 1676 7380 Ox 2 eel 0.44 Nay © ills) 
1528 0.25 3. 0.47 OM ae e248 7789 0.86 1 1.05 AZO 50 
pF 61662 0.41 S002 100) 35 220.0 7790 OF49 TS 0.64 14 2 20.0 
1912 38 O25 0 0.44 4 2 24.8 
1960 36 —0.02 4 0.03 -10 4 24.8 


of a more transparent filter, i.e., one with a broader 
transmission band. For this it was assumed that the 
reduction to the standard system was the same as that 
for the less transparent filter used for the brighter 
objects. As a preliminary test of this assumption, two 
different infrared filters were used, and accordant mag- 
nitudes were obtained. Since this check was necessary 
but not sufficient, we also used the Lal II infrared- 
sensitive multiplier with the stellar photometer provided 
with the GG 11+ BG 18 yellow filter, the BG 21+RG 8 
regular infrared filter used on bright objects, and finally 
the RG 2 extra-transparent filter with which the J 
magnitudes of the M31 clusters were obtained. We 
then measured a number of bright objects with the 
yellow filter and with both the normal and lighter 
infrared filters. These “natural’’ photoelectric colors, 
(v-7), were then compared with others on the standard 
(V-I) system, with the results given in Table V and 
shown in Fig. 3. The upper set of points is for the more 
transparent J filter, the lower for the regular one. From 
the fact that observations made with either filter give 
points that fall on a single straight line for both stars 
and clusters—as indicated in the figure—we conclude 
that the J magnitudes for the fainter M31 clusters 
probably have no large (v-7) to (V-Z) transformation 
errors due to the use of a more transparent filter than 
that used for the brighter galactic globular clusters. 


VII. OBSERVATIONAL MATERIAL 


The observations made at Mount Hamilton were 
begun in July, 1953, and were completed in August, 
1957; those at Mount Stromlo were obtained during 
September through November, 1951. The principal 
results are collected in Table VI for the galactic objects 
and in Table VII for the extragalactic. 

In Table VI, the column headingsare self-explanatory ; 
Table VIb contains three-color data for 28 open 
clusters included mainly for purposes of comparison 
with globular clusters. In most cases, these open 


clusters were observed with photometric apertures suf- 
ficiently large to give a satisfactory integrated color of 
a majority of the cluster members. 

In Table VII, the first column contains identifying 
letters and numbers that mean the following: The letter 
H, in M31 with arabic numbers and in NGC 205 with 
roman numerals, designates clusters found by Hubble 
(1932). The letter B denotes those found by Baade, 
who communicated the identifications for clusters with 
numbers of 200 to 300 to Nassau and Seyfert (1945), 
and who also kindly sent us a marked chart with later 
discoveries, to which he assigned numbers of 300 and 
larger, and to which we gave numbers greater than 140 
and less than 200 for unnumbered objects; M means 
those found by Mayall and discussed by him and 
Eggen (1953), except for M V, which is a later identi- 
fication; the Anon. object at the end of the table of 
M31 objects is one of two that Hubble (1932, footnote 
p. 54) considered as “probable.” For those in M33, we 
used an unpublished print marked by Hubble to show 
the objects we list as Ha, Hb, and Hd, while Mc is an 
unpublished identification by Mayall. In NGC 185, we 
observed two of the three reported by Baade (1944), 
and Ba is our designation for the object on his chart 
(1944, Plate IV) with coordinates 33.0 mm £ and 
37.4 mm N of the center of the system, and with Bc 
the one 12.4 mm E and 3.5 mm S. These objects in 
M31, some not previously charted in the literature, are 
marked in Figs. 4 and 5. In the Fornax dwarf system, 
the Anon. cluster is the one 6/7 E and 1/5 S of CD 
—35°919 (Baade and Hubble 1939). The three-color 
measurements are not so complete as we should like, 
mainly because they represent results so close to the 
useful limit of our infrared photometer that they could 
be obtained only under the very best observing con- 
ditions. 

An additional group of clusters of considerable 
interest is that in M33. With the Crossley, these objects 
are faint and difficult to see, and for this reason we were 
able to measure only four (Fig. 6) of more than 100 
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discovered by Sandage (Bowen 1956). We shall therefore 
mention here only the fact that these four M33 clusters 
are unusually blue. The two bluest objects may be rich 
Population I clusters like some of those in the Large 
Magellanic Cloud. The nature of the other two cannot 
be determined from the present data, but if they are 
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Fic. 4. Outlying star clusters observed in M31. 
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Population II objects, their stellar content must k 
considerably different than that of globular cluster 
found in the other major systems. Of the two clustet 
measured in NGC 205, H III seems normal, whil 
H VIII seems blue and faint like the clusters in M3. 
The clusters Ba and Bc in NGC 185 are of norm: 
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Fic. 5. Star clusters observed over or in main body of M31 and around NGC 205 (upper right). 


observed color, but they are of rather low luminosity, 
unless NGC 185 is somewhat more distant than M31. 
Their lower luminosity, however, is consistent with the 
concept that the brightest objects in a low-luminosity 
parent system may be systematically fainter than those 
in a system of normal or high luminosity. 


VIII. DISCUSSION OF THE THREE-COLOR 
OBSERVATIONS 


Stars 


The three-color data for the relatively unreddened 
stars in Table III are plotted in Fig. 7. The main 
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sequence, and the giant and dwarf separation for late- Fig. 8 are plotted the three-color individual observa- 
type stars, at about spectral type KO, may readily be _ tions for the 12 F- and G-type supergiants that show 
seen. This plot therefore establishes rough thermal loci a large range in reddening (Kron 1958). The latter 
for relatively unreddened, probably single stars. In define a reddening line for objects whose intrinsic colors 


TasLe VII. Photometric data for extragalactic star clusters. 
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HVIII 16.72 2 +0.44 Sh eo. pice yous Nets M33 
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« Outlying objects (Fig. 4). 
b Object No. 81 in Harvard Circ. No. 276. 
© Cluster 6’.7E and 1’.5.S of CD-35°919; Baade and Hubble (1939). 
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Fic. 6. Star clusters observed in M33. 


may be comparable with those of the globular clusters. It is noteworthy that the adopted color system has 


The line is straight within the accuracy of the ob- good resolution within the color range represented, for 


servations. A least-squares solution for an assumed €Ven among these few supergiants the effects of space 
linear relationship gave reddening are well separated from those of thermal 


i, reddening. As the figure shows, the later spectral types, 
(P-V)=+0.20+0.05+ (0.77+0.047) (V-Z). G2, lie above the least-squares line, while the earlier 
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Fic. 7. Three-color (two color-index) diagram for unreddened stars observed to define main, 
dwarf, and giant sequences in P, V, J system. 


types, F6 to F9, lie below it; the unlabeled points all 
correspond to type GO. 


Reddening/Absorption Ratios 


With the six-color [V-] data in Table IIIvfor the 
F- and G-type supergiants, least-squares solutions of 
the form (P-V) and (V-J)=a+0[V-I] were made, 
since the plots of (P-V) and (V-J) against [V-Z ], Fig. 9, 
show satisfactorily close linear correlations. Solutions 
for the 12 stars tabulated, and for only the 7 GO Ib 
stars, agreed so closely that only the result for all the 
stars is given here, with probable errors: 


(P-V)=+0.642+0.017+ (0.330+0.016)[V-I)]. 


Hence 

ALV-I]= (1/0.330)A(P-V), 
or 

Ew_n= (3.0340. 15)Ewp_y). 


According to Morgan, Harris, and Johnson (1953), an 
investigation (unpublished) by Harris and Johnson gave 
the following fundamental relation between total visual 
absorption, Ay, and the six-color Ejy_n_ excess: 
Ay=0.96+0.05)E,y_n. On this basis we find: Ay 
= (2.9140.22)E(p_y), which agrees within the limits 
of error with the result, Ay=(3.040.2)E(z_y), found 
by Morgan, Harris, and Johnson (1953) from much 
more extensive material involving (B-V) colors. By an 


exactly similar procedure, we find from the (V-J) 

colors of the same 12 supergiants that A4y= (2.33+0.14) 

XL w_n. j 
Clusters 


In Fig. 10 are plotted all the three-color data that 
have resulted from the present study of star clusters 
in the Galaxy. The data are taken from Table VI. The 
main sequence, dwarf and giant branches from Fig. 7, 
and the reddening locus from Fig. 8 for the supergiants, 
are drawn in Fig. 10 as light lines. To obtain a reddening 
line for the galactic globular clusters, linear least- 
squares solutions for the regression lines were made for 
the globular cluster three-color data given in Table VI 
and plotted in Fig. 10, with the average result: 


(P-V) = +0.022-£0.025+ (0.706+0.030) (V-J). 


Although probably not significant, the slope happens to 
agree closely with that found for the F- and G-type 
supergiant stars. Several features are apparent. 

(1) The open clusters in the Galaxy separate fairly 
well from the globular clusters in this three-color 
photometry, even in the presence of appreciable 
reddening. 

(2) The distribution of the points for all objects is 
parallel to the reddening lines; this result suggests that 
the most important factor contributing to the elongated 
distribution of points in the diagram is space reddening. 

(3) There is some indication that open clusters may 
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Fic. 8. Three-color (two 
color-index) plot and red- 1.0 
dening line for 12 F- and 
G-type supergiants. 
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be classified according to some parameter, which may 


be age, from three-color photometry of their integrated 


_ light. Some older open clusters like M67 and NGC 7789, 


: 


for example, plot among the globular cluster group, 
whereas younger open clusters like NGC 869 (h Per) 


plot below it. 
(4) The use of an infrared color may be helpful in 


(V-I) 


PHOTOMETRY 


597 


<— (P-V)=+0.20+0.77 (V-I) 
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distinguishing clusters in which there are some bright 
red stars. Thus, among the galactic clusters whose 
(V-I) colors place them systematically redward of the 
others, all but one, NGC 559, are known to contain 
high-luminosity red stars. Unfortunately, the presence 
of bright red stars, or the effect of high space-reddening 
are not easily separable. 
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Fic. 9. Relationships between (V-J) and (P-V)"vs six-color [V-J] indices for 12 F- and G-type supergiants used to 
evaluate the ratio, between total and selective absorption. 
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Fic. 10. Three-color (two color-index) diagram for all observed galactic open and globular clusters. 


These results for galactic clusters are only of a pre- 
liminary and exploratory nature, but they suggest that 
further photometry of their integrated light ‘may be 
worthwhile, especially if comparisons are to be made 
with similar objects in other galaxies. 


Spectral Types 


As a preliminary step to study any possible corre- 
lation between the three-color data and spectral types 
of globular clusters, we next examine a group probably 
having little or no reddening. This group was formed 
by selection of the 22 observed clusters at high galactic 
latitudes, || 220°, and for which spectral types were 
available. The three-color data for these objects are 
plotted in Fig. 11. In this figure it is clear that there is 
no strong correlation between these color data and 
spectral type, a result that is not surprising in view of 
the composite nature of the stellar content of globular 
clusters. Nevertheless, two A-type clusters are grouped 
together in the lower left-hand corner, whereas F-type 
clusters tend to lie above and to the right of them. The 
only high-latitude G-type cluster observed is NGC 
6981, which in the figure falls among the F-type clusters. 
NGC 6681 and 6723 are other G-type clusters, not at 
high latitude but apparently little reddened, that would 


plot near NGC 6981. Four clusters in the upper right- 
hand part of the plot appear to be well separated from 
the main group. For one of these, NGC 6254, however, 
Arp (1955) found a color excess of Hc1=0.4 mag., more 
than enough to account for its position in the diagram. 
This case suggests that the other three clusters in this 
part of the diagram may also be there because of red- 
dening. This possibility is strongly supported by the 
fact that three of these clusters, NGC 6171, 6218, and 
6254, all lie rather close to each other near /=340° 
and 6=-+22°, which is a region of low galaxy counts 
(Shane and Wirtanen 1954). Furthermore, by means 
of a method (Kron 1958) using the six-color photometry 
of Stebbins and Whitford (1945), we find that the K3 
giant, HD 152601, which lies within 2° of NGC 6254, 
has a (P-V) color excess of between 0.2 and 0.3 mag., 
even though its distance is doubtless much less than 
those of the globular clusters. The fourth cluster, NGC 
6864, may be a similar case of appreciable reddening. 
It lies in the same region as the G3 Ib supergiant star 
HD 192876 (a! Capricorni), which has a (P-V) color 
excess of +0.18 mag. (Kron 1958). 

We conclude that globular clusters, even if allowance 
is made for space reddening, have intrinsic properties 
that cause dispersion both perpendicular and parallel 
to a reddening line in a plot like Fig. 10. Since certain 
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TABLE VIII. Comparison of spectral type estimates. 
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of these intrinsic properties are doubtless indicated by 
spectral type, we decided to try to correct for the loose 
correlation between color and type in estimating color 
excesses from these three-color data, as described in the 
next section. 


Ix. COLOR EXCESS AND ABSORPTION 


As the first step in using spectral types to estimate 
color excesses and absorptions for the galactic globular 
clusters, we compare the previously available, most 
extensive set of types (Mayall 1946), based on Mount 
Wilson standards, with those more recently determined 
by Morgan (1956), referred to the Yerkes MK system. 
For some clusters, Morgan’s estimates depend on the 
old, low-dispersion Crossley data, while for others they 
come from new, higher dispersion material obtained 
with the Cassegrain spectrograph of the 82-inch 
McDonald reflector, as indicated in Table VIII. The 
comparison shows no appreciable systematic difference 
between the estimates depending on the Crossley 
plates, and only a tenth of a spectral class for the 
McDonald types. Since the latter are undoubtedly 
preferable, all the 1946 types were increased by one- 
tenth spectral class, and, where a 1956 McDonald type 
was available, it was averaged in with double weight. 
These revised types are given in the last column of 
Table VIII, and in the second column of Table X. 
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Fic. 11. Three-color, or two color-index plot for globular clusters of known spectral type and for which |}| 2 20°. 
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Fic. 12. Plots of (P-V) and (V-Z) against spectral types of globular clusters; those in less obscured regions are indicated by 
filled circles with associated numbers denoting field galaxies counted, those in obscured areas by open circles; the curves marked 
T and II are the two assumed alternative relationships between intrinsic color and spectral type, used to estimate color excesses. 


In the second step, the (P-V) and (V-J) colors were 
plotted separately against the revised spectral types, 
as shown in Fig. 12. Those clusters probably least 
affected by galactic absorption are represented. by filled 
circles and associated numbers that indicate the field 
galaxies counted on one-hour exposures with the 


TaBLE IX. Assumed intrinsic colors for spectral types. 


Spectral (P-V) - Spectral (V-D) 
Type I II Type JHE 
A6-F0 0=48 A6-A8 0™67 
F1-F2 0.49 A9-F2 0.68 
F3-F4 0.50 F3-F4 0.69 

FS 0™m51 0"50 F5 0.70 

F6 0.53 0.51 F6 0.71 

F7 0.55 0.52 F7 0773 = 0™72 
F8 0.57 OSS F8 0.75 0.73 
F9 0.59 0.53 F9 0.78 0.75 
GO 0.62 0.54 GO 0.81 0.76 
G1 0.65 0.55 G1 0.86 0.78 
G2 0.69 0.56 G2 0.92 0.79 
G3 0.72 0.57 G3 0.97 0.81 
G4 0.75 0.58 G4 1.02 0.83 
G5 0.79 0.59 G5 1.07 0.84 
G6 0.82 0.60 G6 1.14 0.86 


Crossley reflector (Mayall 1946, Table 5); the ope 
circles denote clusters in fields where no galaxies wei 
counted. Several clusters are identified by their cat 
logue numbers: NGC 6440, the reddest one measure 
by us and by Stebbins and Whitford (1936) falls fe 
outside the diagram; NGC 6522, whose color was use 
by Baade (1951) to estimate the absorption an 
distance toward the galactic center; NGC 6171, tk 
reddest cluster in a field in which a few galaxies coul 
be counted; NGC 6723 and 6681 apparently litt 
reddened despite their presence in fields showing r 
galaxies; NGC 6981, whose estimated spectral tyy 
seems somewhat too late for its color; and NGC 665 
and 6864, which were used as limiting cases to draw tl 
“bluest-color” bounding curve I, described in tt 
following. 

As might be expected for composite stellar system 
Fig. 12 shows no close dependence of color upon spectr: 
type. Nevertheless, there appears to be a gener. 
tendency, among the least obscured clusters, for tk 
later types to be redder. Unfortunately, however, the: 
are so few clusters of latest type in unobscured fielc 
that the picture is badly confused in the range fro1 
GO to G5. Thus any attempt to obtain a curve the 
would relate the ‘‘unreddened,” or “normal,” « 


TABLE X. Revised spectral types, color excesses, and absorptions. 
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0.53 0.58 0.9 sleal 12 ies 1.0 lee, 0.34 0.41 
Ona 0.22 0.6 0.8 0.4 0.5 0.5 0.6 Oily 0.21 
0.37 0.53 0.7 de 0.9 iL 0.8 LD 0.28 0.41 
0.00 0.1 0.2 0.0 0.0 0.1 0.00 0.03 
0.27 0.29 0.6 0.7 0.6 0.7 0.6 0.7 0.21 0.24 
0.13 0.15 0.2 0.3 0.3 0.3 0.2 0.3 0.07 0.10 
0.28 0.36 0.4 0.7 0.6 0.8 0.5 0.8 0.17 0.28 
0.63 0.81 0.8 il. J 1.6 1.9 i7 15 0.41 0.52 
0.53 ie? ey 12 12 0.41 0.41 
0.28 0.41 0.6 1.0 0.6 0.9 0.6 1.0 0.21 0.34 
0.66 0.89 1.0 1.6 1.5 2.0 ee 1.8 0.41 0.62 
0.48 ileal ical iloil igi! 0.38 0.38 
0.66 0.85 1.4 1.9 a5) 2.0 1.4 2.0 0.48 0.69 
0.74 0.97 aul iL sy Dp 1.4 2.0 0.48 0.69 
1.07 1.30 1 Beth 2.5 3.0 2.3 2.8 0.79 0.96 
0.59 ileal i ig 12 0.41 0.41 
0.05 0.1 0.1 0.1 0.1 0.03 0.03 
0.35 0.58 12 iL Ss} 0.8 1.3 1.0 1.6 0.34 0.55 
0.69 0.92 1.8 2.4 1.6 2.1 ag De 0.59 0.76 
0.36 0.55 0.9 il 4 0.8 iL.@) 0.8 1.4 0.28 0.48 
On8We TACs ei A varls7e AO Teele 4s isie 230s 07525 05692 
0.70 0.78 1.4 il? 1.6 1.8 15: 1.8 0.52 0.62 
0.37 0.42 0.7 0.9 0.9 1.0 0.8 1.0 0.28 0.34 
1.44 1.67 3.0 305) 3.3 3.8 One 3.6 1.10 1.24 
0.83: 1.06: awe  alyer 1 Oe 2/3453 TSS ES 0.455) 2102622 
ibe 1.56 3.0 S50) S)58) 3.6 ae, 3.4 1.10 ibe aly 
0.63 0.65 ils 1.4 1.4 125 1.4 1.4 0.48 0.48 
0.49 0.72 1.4 2.0 heal lef ie? 1.8 0.41 0.62 
1.42 1.47 3.0 Bigs) S186) 3.4 S52 3.4 1.10 1517 
0.93 1.06 ile Dreal Ail 2.4 1.9 Di 0.65 0.76 
1.10 iS) eed DAS DS} Sal Qed Deol 0.72 0.93 
0.49 0.72 0.9 15 ibs) LAME 1.0 KS 0.34 0.55 
0.27 0.50 0.6 1.2 0.6 1a | 0.6 il % 0.21 0.41 
0.42 0.50 0.9 ee 1.0 ileal 1.0 2 0.34 0.41 
0.11 0.34 0.2 0.8 0.3 0.8 0.2 0.8 0.07 0.28 
0.44 0.63 0.6 lea! 1.0 1.4 0.8 2) 0.28 0.41 
0.00 0.19 0.1 0.6 0.0 0.4 0.0 0.5 0.00 0.17 
0.69 0.70 0.9 1.0 1.6 1.6 a i538 0.41 0.45 
(0) 0.06 (0) 0.1 (0) 0.1 (0) 0.1 (0) 0.03 
0.52 0.75 0.7 153 2 a7 1.0 1.5 0.34 OR52, 
0.25 0.27 0.4 0.5 0.6 0.6 0.5 0.6 0.17 0.21 
(0) 0.09 (0) 0.1 (0) 0.2 (0) OZ, (0) 0.07 
1 ehil 1.36 Asi) 2.9 3.0 Sad 2.8 3.0 0.96 1.03 
0.33 0.6 0.7 0.8 0.7 0.8 0.24 0.28 
OEtS Ord 0.1 0.3 OZ 0.2 0.07 0.07 
0.20 0.48 0.4 1.0 0. i 0.4 1.0 0.14 0.34 
0.17 0.30 0.1 0.4 0.4 0. 0.2 0.6 0.07 0.21 
0.03 0.08 0.0 0.2 0.1 0. 0.0 0.2 0.00 0.07 
(0) 0.00 (0) 0.0 (0) 0. (0) 0.0 (0) 0.00 
0.16 i023 0.4 0.4 0.4 0.14 0.14 
0.02 0.3 0.4 0.2 0.2 0.07 0.07 
0.04 0.2 0.1 02 0.2 0.07 0.07 
0.05 0.0 0.1 0.0 0.0 0.00 0.00 


* Types in parentheses are assumed. 
b New estimate by Morgan (1956). 
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“intrinsic” color of a globular cluster with its estimated 
spectral type is bound to involve considerable uncer- 
tainty and judgment—at least with the present material. 
This uncertainty, of course, enters directly into the 
deduced color excesses, which are obtained as the dif- 
ferences between the measured colors and those given 
by the assumed or inferred “normal color”-spectrum 
relationship. Compared to this source of uncertainty, 
the errors in the measured colors probably are of minor 
importance in color excess estimates. 

The third and last step in deriving color excesses was 
to draw in Fig. 12 curves representing such “normal 
color’’-spectrum relationships for both the (P-V) and 
(V-I) colors. Obviously, such curves could be drawn 
in several ways, such as an average through the pre- 
sumably least reddened clusters (filled circles) which 
would yield a number of negative excesses, or as 
“‘bluest-color” bounding envelopes that give either zero 
or a small color excess for the highest latitude clusters. 
While for the latter some color-magnitude diagrams of 
individual stars have been interpreted as indicating a 
small amount of reddening, we have preferred to risk 
under- rather than over-correction. We have, therefore, 
drawn limiting curves through the bluest clusters, as 
shown in Fig. 12 by the single lines for types A6 to F5 
for (P-V), or to F6 for (V-Z). Since for later types it 
was difficult to know whether or not to place full faith 
in the estimates of spectral types for NGC 6681 and 
6723, we decided to use the two alternatives repre- 
sented by the branching curves labeled I and II. The 
first ignores the three clusters in question, while the 
second gives them full weight. NGC numbers identify 
those clusters used to define the curves, which are 
dashed near their extremities to show that their course 
is conjectural because of insufficient data. The colors 
corresponding to these curves are listed in Table IX, 
and they are given to 0.01 mag. only to facilitate further 
computations. Their probable errors are doubtless much 
larger, because of the subjectivity of the procedure. 
We mention as reasonably conservative, estimated 
probable errors of +0.05 mag. for types A5 to about 
GO, and of +0.1 mag. for GO to G5. We also note that 
the ranges in the tabulated colors are: for alternative I, 
0.34 mag. in (P-V) and 0.47 mag. in (V-J); for II, 
0.12 mag. in (P-V) and 0.19 mag. in (V-Z). Thus alter- 
native II is not far from the assumption of the same 
intrinsic color for clusters of all spectral types—an 
assumption whose usefulness was demonstrated nearly 
25 years ago by Stebbins and Whitford (1936). Alter- 
native I, on the other hand, allows for a fairly large 
variation in color. Depending on which alternative is the 
more tenable, color excesses estimated with I may be 
too small, those with II too large; perhaps the truth 
lies somewhere between. But because we cannot be 
sure, we shall in the subsequent discussion carry through 
parallel calculations for each alternative, in order to 
show the range of values in such quantities as corrected 
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colors, total magnitudes, distance moduli, and diam- 
eters. 4 
The foregoing procedure was used to compute color 
excesses for all clusters of known spectral type, and the 
results are given in Table X under the column headings: 
Ew_y) and Eq_p. For clusters of unknown spectral 
type, F5 was assumed for those with |5| 220°, GO for 
those with —20°<b<+20°, and G5 for those i in the 
group around the galactic center direction; these 
assumed types are in parentheses in the table. 

Table X also contains the visual absorptions, Ay, 
obtained from 2.9 E.p_y) and 2.3 Eqw—r) (Sec. VIII). 
Since Ay should be the same regardless of the spectral 
regions used for its estimation, we have a useful internal 
check on the material and procedure. If the two poor-. 
data cases of NGC 6401 and 6453 are disregarded, the 
frequency distribution of the differences without regard 
to sign, |AAy| =a AE: Ew-1n—-2.9 Ew_y), is as follows, 
for 65 clusters: 


|AAy| 0.0 0,1 0.2 0.3 0.4 (0.5) (OGRU Sven Ocm 
I 14 Ae 11 11 5 3 1 1 2 
Il (RYE FA 10 ii 8) B 1 0 2 


Thus there is a satisfactorily high concentration of 
small differences, which amounts to 94% (I) and 95% 
(IL) for |AAy| <0.5 mag. The two largest differences 
of 0.8 mag. are for NGC 6266 and 6553; for the former 
the (V-I) measurements have the unusually large AD 
of 0.11 mag.; for the latter, the spectral type is assumed 
as G5, but a still later type would give more consistent 
color excesses by this procedure. If averages of all 65 
differences are taken, the mean difference without 
regard to sign is 0.21 mag. (I) and 0.18 mag. (II); with 
sign, +0.15 mag. (I) and +0.09 mag. (II). The size 
of the latter systematic differences suggests, as one pos- 
sibility for better agreement, the use of slightly different 
total-to-selective absorption ratios. For example, if 3.1 
for (P-V) and 2.1 for (V-Z) were used, the systematic 
difference between the V absorptions derived from the 
two color indices would be close to zero. Another possi- 
bility is, of course, to use somewhat different curves in 
Figure 12. Such refinements seem neither justifiable or 
desirable to us, however, so we have taken averages for 
alternatives I and II of Ay obtained from (P-V) and 
(V-I); these are listed as (A)y in the next to the last 
two columns of Table X. Finally, we computed from 
(A)y/2.9 the values of (£)(p_y) listed in the last two 
columns. 

Since the procedure used here to derive color excesses 
for the galactic globular clusters is admittedly arbitrary, 
it is of interest to compare the present results with 
others recently obtained by different means. Usually 
the latter have involved color-magnitude diagrams 
determined from individual cluster stars, with the 
space-reddening estimated from the position of the 
giant branch, of the main sequence, or of the RR Lyrae 
variable-star “‘gap.’”’ Most of the modern photoelectric 
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Fic. 13. Cosecant-law repre- 
sentations of globular cluster 
color excesses estimated from 
alternatives I and IT. 
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material refers to the larger, brighter and little 
reddened high-latitude clusters, so that we cannot get a 
good check on the larger color excesses. If we limit the 
comparison to the best-observed clusters, we find the 
results listed in Table XI. 

We consider this comparison to be generally satis- 
factory, probably to a precision of about 0.05 mag., 
except possibly for the case of NGC 6254 (M10). This 
cluster has no known RR Lyrae variables, and the color- 
magnitude diagram depends on only one 60-inch 
reflector plate in each color. Under these circumstances, 
the photometric estimate of color excess may be as 
uncertain as the one based on its spectral type, by 
+0.1 mag. 

It is also of interest to examine the correlation 
between galactic latitude and the color excesses derived 
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here. We assume that the actual patchy distribution of 
absorbing material may be treated statistically as a 
thin, uniform layer centered on the galactic plane and 
sun. If the clusters are outside such a layer, then the 
reddening and absorption in magnitudes vary linearly 
with cosecant b. Figure 13 shows a plot of E,p_y) for 
alternatives I and II, with the lines resulting from the 
least-squares solutions, which with probable errors are 
(with 6440 omitted because of its dominating weight): 


I, Ew_y)= —0.018-+0.022+ (0.052-£0.003) csc] oI, 
Il, Evp_y)= +0.001+0.022+ (0,063+0.003) esc|b}. 


Similar results were obtained many years ago by 
Stebbins and Whitford (1933, 1936, 1937) and, to make 
a comparison, we reduce their data to the P, V system, 
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as follows. Stebbins’ original (1933) solution for 42 
clusters gave E.= (0.0237+0.0018) csc|b|, which was 
later (1936) discussed in terms of the cluster color 
scale E. Since this scale is equivalent to £; for stars, 
and £,=1.68 E2, then £,= (0.040+-0.003) csc|b|. Also, 
in their first report (1937) on the colors of extragalactic 
nebulae, they obtained E,= (0.015-+40.005) csc|b| as 
the average from two solutions of 66 nebulae and of 19 
globular clusters outside Hubble’s zone of avoidance. 
Moreover, Stebbins and Whitford (1936) found that 
1.6 E,= Ecr, we assume that Ecor= E:p_y), and Morgan, 
Harris, and Johnson (1953) deduced that Ay=6.1 F). 
These relationships provide the comparison in Table 
XII of these photoelectric determinations of the redden- 
ing and absorption toward the galactic poles. 

Considering all the assumptions and uncertainties 
involved, especially the difficult decision of what to 
use for a ‘normal” color of a globular cluster, we con- 
sider that both the older and newer data agree well 
enough to confirm Stebbins and Whitford’s earlier 
conclusion: on the average, there is very little space- 
reddening and absorption toward the galactic poles. 

We wish to emphasize that this result based on the 
assumed cosecant law is mainly useful in comparisons 
of photometries of the globular cluster system as a 
whole; computed cosecant law color excesses for in- 
dividual objects have little significance because of the 
very patchy, irregular distribution of interstellar 
matter. 


xX. LUMINOSITIES, COLORS, DISTANCES 
AND DIAMETERS 


With the color excesses and absorptions deduced as 
described in the preceding section, we may obtain 
corrected apparent total magnitudes, integrated'colors, 
distances and linear diameters. But first we shall 
derive a set of absolute total magnitudes, or luminosities, 
since they may be obtained independently of absorption 
whenever the distance moduli depend only upon ap- 
parent magnitudes of brightest stars or variables 
measured in the same spectral region. Such a list of 
apparent distance moduli, based on the assumption 
that M,,=0.0 for RR Lyrae variables, has recently 
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Evp_v) 

NGC M I II Ecx Reference 

4147 0700 0700 <0™15 Sandage and Walker (1955) 
4147 0.00 0.00 0.00 Newburn (1957) 

52/2 oa O00 Ons 0.00 Arp (1955) 

5904. 5 0.07 0.07 0.16 Arp (1955) 

6205 13 0.00 0.03 0.16 Arp (1955) 

6254 10 0.17 0.28 0.4 Arp (1955) 

6341 92 0.03 0.03 0.00 Arp (1955) 

6656 22 0.41 0.45 0.4-0.5 Arp and Melbourne (1959) 
LOTSa pop 0207) 0107 0.00 Arp (1955) 

7089 2 0.07 0.07 0.05 Arp (1955) 
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TABLE XII. Polar reddening and absorption. 


No. Cl Sol. Eve_v) Ay 
64 I 0"052-++-0™003 0™15+0701 
64 II 0.063+0.003 0.18+0.01 
42 (1936) 0.064-+0.005 0.24+0.02 
19\ (1937) 0.024+0.008 0.09-40.03 


given in the third column of Table XIII, euenpe that | 
the moduli for NGC 2419, 6426 and 6638 represent | 


revisions or additions based on later work. 


Luminosities 


Since most of the moduli depend on photographic } 
observations, we give in the sixth column of Table XIII | 
the absorption-free, total absolute photographic mag- } 
nitudes, Mp;, obtained from V; in Table IV, from } 
(P-V) in Table VI, and from (m—M) p in Table XII; | 
Vit (P-V)— (m—M ) po: Also, because mucll { 
modern photometry of clusters is being done in visual } 


thus Me,= 


light, we have derived a set of apparent, visual-light 


moduli, (m—M)y in the fourth and fifth columns of | 
Table XIII, from (m—M)9—(E)p_v). To correspond | 
with these moduli, we also computed total absolute | 
visual magnitudes, Mv; in the seventh column, from | 
Mr,—(P-V).; for this purpose the corrected colors, | 
(P-V), in the eighth and ninth columns, were averaged } 
for alternatives I and II, since the maximum difference } 


is only 0.21 mag. 


The values of Mv, which involve measurements in | 
both photographic and visual light, are obviously | 
absorption-free only if the color excesses have been | 
correctly estimated. We have a partial check, however, ‘| 
from the few globular clusters having modern, visual- 
light photometry (Table XI). Thus with M’y, obtained | 
from V; and the directly observed, visual-light apparent 
moduli, (m—M)’y, we obtain the comparison with Mv; | 
shown in Table XIV. Since the average differences with | 
and without regard to sign are only —0.01 and 0.14 | 
mag., with none greater than 0.3 mag., the comparison © 
seems satisfactory for these clusters. As mentioned | 
previously, however, this small group includes only one © 


appreciably reddened cluster (NGC 5264), so this check 


on the estimated color excesses, while necessary, is — 


hardly sufficient. 


The frequencies of the total absolute photographic | 


and visual magnitudes are given in Table XV and are 
shown as histograms in Fig. 14; clusters omitted are 


NGC 2419, 5634, 5897, 5401, 6416, 6453, and 6544 _ 


because of uncertain V;, and 6838 because of uncertain 
(m—M) »,. This material has also been subdivided into 


three groups in order to look for any systematic varia- — 


tion of cluster luminosities within the Galaxy. These 
groups are (1) outlying clusters with |b| > 10°, (2) those 
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TaBLE XIII. Apparent moduli, luminosities, corrected colors,* and magnitudes.® 


(m—M)yv (P-V)c 
M (m—M) pg 1 II Mp, Mv; i II 
19.4 19.3 <-7.9 <-—8.4 0™50 
ileal if sal —6.4 —6.9 0.48 
68 16.1 1.7 —7.2 —7.7 0.49 
53 16.7 16.7 —8.4 —8.9 0.50 
16.4 16.4 —6.0 =6-5 0.52 
3 1584 Si) 15.6 —8.9 —9.4 0.56 0.49 
16.4 16.3 —6.6 —7.1 0.50 
seal 17.1 17.0 <-—6.7 <-—7.2 0.53 0.49 
Viel 17.6 —6.8 —7.3 0.47 
: ° oo 0.49 0.46 
16.2 16.1 ELE NEMD 0.50 
5 io ilsypal —8.6 =9.2 0.56 
80 16.0 15.9 15.8 —7.9 — 8.9 0.62 0.55 
4 14.0 13.7 13.6 —7.2 —7.8 0.58 0.51 
16.0 15.8 —O-2 == (ats) 0.64 0.61 
107 16.2 15.9 15.8 —7.0 —7.6 0.69 0.56 
13 14.8 14.8 —8.4 —9.0 0.57 0.54 
12 15.2 15.0 —7.7 —8.2 0.56 0.53 
Not 17.6 —7.6 —8.2 0.57 0.54 
10 15e2 15.0 14.9 —7.8 —8.4 0.62 0.51 
62 16.7 16.3 16.2 Ot —9:6 0.54 0.43 
19 16.0 15.6 —8.2 —8./ 0.49 
tidied 17.1 17.0 —7.4 —8.0 0.70 0.57 
Wes 16.9 16.7 —6.7 —7.2 0.72 0.51 
16.7 16.3 —7.4 —7.9 0.48 
s . 0.69 0.48 
16.4 16.0 —7.7 —8.2 0.46 
97 Te 52, nil —8.6 0.50 
0.82 0.65 
Lin? 16.9 16.7 —7.7 —8.4 0.78 0.58 
0.65: 0.48: 
14 16.7 16.2 16.1 —8.0 =8..6 0.60 0.50 
17.8 17/55 <-4.2 <—4.8 0.57 0.51 
.- tee 0.71 0.57 
07595 | O42: 
17.8 NESE) —8.2 —8.7 0.54 
0.57 0.50 
0.62 0.51 
0.67 0.46 
0.77 0.56 
28 15.6 T5539 15a —7.7 —8.3 0.63 0.56 
WES) Wl? INF! —7.4 —8.0 0.68 0.55 
Pere nag Sci age 0.78 0.61 
22 14.2 13.8 —8.2 —8.6 0.45 0.41 
16.5 16.2 16.0 —7.3 —7.9 0.70 0.52 
54 17.4 722 —9.0 —9.5 0.53 0.49 
1573 “2° 15e2 —7.5 —8.0 (0.61) 0.54 
17.4 16.4 —6.7 —7.2 0.59 0.52 
56 16.3 16.1 16.0 —7.3 —7.8 0.50 0.46 
55 14.5 14.4 —7.7 —8.2 0.48 
71 14.6: 14.5: 14.3; —5.4: —5.6: 0.81 0.61 
75 18.1 18.0 17.9 —8.8 —9.4 0.64 0.50 
16.9 16.9 16.8 =—7.2 —7.8 0.62 0.55 
72 17.0 see 17.0 —7.0 —7.8 (0257) Ons 
18.8 18.7 —7.5 —8§,0 0.47 
15 15.9 15.8 —8.9 —9.4 0.52 
2 16.2 16.1 —9.1 —9.6 0.50 
30 1525 15.5 —7.4 —7.2 0.48 
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a Numbers in parentheses are uncorrected (colors) or uncertain (magnitudes), 
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No. No. 
5 Outlying 
|b] > 20° 
2208 
{e) 
On, 85 —9.0 -8:5 -80 -75 —-7:0 =—65 =60 My, -!0.0 -9.5 9.0. -85 ~80 ‘-7.5  =fORe=s.0 
5 5 
Intermediate te 
— 20%b<+20° 4 
10 ® 
| 
ie) ie) Fa 
Mi=—95 —9.0 —85 —80 —7.5 ~—/.0 “=co 60 My,—!0.0 -95 -90 -—85 -80 -7.5 -7.0 -65 
5 
Center 
group 
10 ® 
oO 
Mp.-9.5 -90 -85 -80 -7.5 -7.0 -65 -6.0 C WGee 10.0 —95 -9.00 —85 -80 =-75°>-703—65 


Fic. 14. Histograms for three groupings of globular cluster total absolute magnitudes, in photographic light (left) and visual 
(right), computed with the assumption that Mp=0.0 for cluster-type (RR Lyr) variables. 


of intermediate concentration with —20°<b<+20°, 
and (3) those in the direction around the galactic 
center. The histograms for the frequencies for the three 
groups combined are given in the lower part of Fig. 20, 
for comparison with similar data for M31 (Sec. XII). 

The small number of clusters in the intermediate and 
center groups, 10 each, is not encouraging for refined 
analysis. Nevertheless, there is little suggestion that 
cluster luminosities are significantly different in different 
parts of the Galaxy. So far as they go, the frequency dis- 
tributions look rather similar for the three groups, and 
the average absolute magnitudes in the penultimate 


TaBLE XIV. Comparison of total absolute visual magnitudes. 


column of Table XV are remarkably the same. Thus the 
over-all means in photographic and visual light, —7.7 
and —8.2, respectively, may be fairly stable, statistical 
quantities as distance indicators, for example, of the 
galactic center (Sec. XIII). 


Colors 


The frequencies of the corrected, integrated-light 
colors, (P-V),, are given in Table XVI, and are shown 
as histograms in Fig. 15; the poor-data cases of 6401 
and 6453 are omitted as for the luminosities, and this 
material also has been subdivided into three groups, 
according to concentration of the clusters toward the 
galactic plane and center. 

Comparison of the histograms obtained from alter- 


, if 

peal MBS NS si aia natives I and II shows, as expected, that alternative IT 
tre Ki ips A Ae aoe Boe a produces frequency distributions having pronounced 
5272 3 Gk AS Er 2293) VOL ae Ons peaks, especially for the outlying and intermediate 
Sons 2 : Mf ~ ES —9.2  —9.2 0.0 cluster groups. The small assumed intrinsic color 

j ? ee ee a: variation also tends to suppress any systematic dif- 
6254 10 6.6 14.8 —8.2 —8.4 0.2 ferences between the groups. For example, the histo- 
on - ee eh e mee noe a i grams resulting from alternative I suggest rather 
7078 «15 6.4 15.9 29.5 ong 0.1 definitely that the clusters become redder, both with 
7089 2 6,5.) Gre —9.7 —9.6 0.1 concentration toward the galactic plane and with 


location around the center. This effect, however, is 
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TaBLE XV. Distribution of luminosities of galactic globular clusters. 


pen eee O69 «8-84 7:9) | 7.4. ° =6.9 —6.4. ‘Total (M) o 
Interval —9.5 SU) 8358) —8.0 —7.5 —7.0 —6.5 —6.0 No. (mag.) (mag.) 
Outlying ie il 4 3 5 5 2 22 —7.7 0.8 
|b] >20° \Mv, 1 5) 2 5 5 2 2 22 —8.2 0.8 
Intermediate Mp; 1 0 2 2 3 1 1 10 —7.5 0.75 
—20°<b<+20°\My; 1 0 3 1 3 1 1 10 —8.1 0.7 
Center { Mp; 1 0 2 2 4 1 10 —7.75 0.6 
Group \Mv; 1 0 2 5 1 1 10 —8.3 0.6 
All i 3 4 7 9 12 4 3 42 7.7 0.8 
Mv: 3 5 df 11 9 4 3 42 —8.2 0.8 


closely related to the way in which color excesses were 
estimated, particularly for clusters of unknown spectral 
type: by the assumption of types: F5, GO, and G5 for 
such clusters in the three groups. Moreover, nearly half 
of the center group clusters are unobserved for spectral 
type. Thus until further photometric or spectroscopic 
work is done for them, we cannot conclude with con- 
fidence from our data that the closer globular clusters 
are to the galactic plane and center, the redder is their 
intrinsic color. But this is an interesting possibility 
pertaining to galactic evolution. 


(I) 


Outlying 
b> 20° 
266 


(P-V), 45 50 55 60 65 70 


Intermediate 
—20°<b < +20° 


16@ 


The frequency distributions of the corrected colors 
for 65 clusters regardless of group are given in the last 
two lines of Table XVI; they are shown as histograms 
in the lower part of Fig. 19, for comparison with the 
colors of the M31 clusters. Here, however, we wish to 
discuss only these combined data for the galactic 
globular clusters. For alternative II, the frequency dis- 
tribution has a fairly symmetrical, prominent maximum 
about a mean (P-V)-=+0.52 mag., with a small dis- 
persion c=0.05 mag. While we do not rule out such a 
result on the basis of the present material, there are 
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Fic. 15. Histograms for three groupings of globular cluster colors corrected for color excess by alternatives I and II. 
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TaBLE XVI. Distribution of corrected colors of galactic globular clusters. 
ne 
(P-V) -+-0.40 0.45 0.50 0.55 0.60 0.65 0.70 0375 (0:80 ©0785.) Dotaledeare. o 
Interval +0.44 0.49 0.54 0.59 0.64 0.69 0.74 0.79 0.84 89 No. (mag.) (mag., 
Outlying f I 7 9 6 3 1 26 «©+0.53 0.06 
[b| > 20° {IL 9 {4 3 26 «640.51 0.03 
Intermediate | (1 ee eee en bei ls 16 +0.60 0.08 
—20°<b<+20° \II 1 2 9 4g 2 16 +0.52 0.05 
Center {I 3 2 0 2 3 5 7 ee cia 2 23 «40.69 0.12 
Group \II 1 5 4 9 1 3 23 +0.55 0.06 
All I 11 13 12 10 4 6 4 3 2 65 +0.60 0. 4a 
II 2 16 27 14 3 3 65 +0.52 0.05 


some reasons to view it skeptically. If the three clusters 
(NGC 6681, 6723, 6981) used to define alternative IT 
are found to have earlier spectral types by future and 
better spectroscopic work, or fail to have their relatively 
blue colors confirmed by further photometry, then this 
particular path of color excess estimation is invalidated. 
If so, there would then be stronger support for alter- 
native I, which conforms much more to the expectation 
that appreciably redder intrinsic colors accompany 


later spectral types. For this reason, we regard mor 
favorably the color-frequency histogram resulting fron 
alternative I. It does not have nearly so high a maxi 
mum, and it is unsymmetrical with a long tail towarc 
redder colors. Moreover, unless the unobserved cluster 
materially modify the frequencies, it seems unlikely t 
us that the most frequent color, about (P-V)-=+0.5: 
mag., or the range from +0.45 to +0.90 mag., wil 
change drastically. The corresponding dispersion 


TABLE XVII. Corrected moduli, distances, linear diameters, and mean surface brightnesses. 
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TABLE XVII—Continued 


(m—M). r (kpc) Do.o(pc) SB)» 
NGC M I Il I II I at ocr 
6402 14 14.7. 14.3 8.7. 7.2 7 23 Sd ain 
6426 167) 16.5 22 20 ut ay, 
e440 Ble Gane af 
6517 
6522 15.9 15 18 ora 
6539 
6544 
6553 

6569 Aa 

6624 , ¢ hy 
6626 28 14.3 14.0 72 163 2 17 Ee Wee Hels bil 
6637 69 ; Se a 
6638 AgtAnL 15¢9 19 15 24 19 Si) eG 
6652 t & 
6656 2 19 6<- 12.5 3.3. 3.2 5. 24 O5% 
6681 70 f a ia 
6712 1S 14.5 11 7.9 39-28 LOT wceonh 
6715 54 16.7 16.6 22 21 31-29 Sais Neon 
6723 Ber ast 10 34 =0:4 
6760 13,6. on 184 ae) Lares 13. 12 hon eens 
6779 56 is4 15.2 12 it 35. 32 Bory, <0 
6809 55 14.2 6.9 42 hy, 
6838 71 n- 048.3 be 8 AG 20: 14: ta.) EGS: 
6864 15 17.8 17.3 36 29 BL bat Pocst Ae sind 
6934 16.9 16.6 24 21 23 20 i oie ited 
6981 72 af 17.0 25 47 40.6 
7006 18.3 46 40 0.0 
7078 15 15.6 13 36 are 
7089 2 15.9 15 30 —2.2 
7099 30 15.5 13 26 0.8 


® Unit of (SB) is My/pc?. 


o=0.11 mag., with such a skew distribution, may not 
be very meaningful. 


Distances and Diameters 


With the values of Ay listed in Table X and of 
(m—M)y in Table XIII, corrected distance moduli 
were computed. These are given as (m—M), in 
Table XVII, with corresponding distances in kilo- 
parsecs, r(kpc). On the basis of these distances, linear 
diameters containing 0.9 of the total light, Do.», were 
computed from the angular diameters do.. in Table IV, 
and, finally, the mean absolute surface brightness, (SB), 
from Mv.+5 logDo.. Results from these calculations 
are plotted in Fig. 16 for alternatives I and II, with 
Mv; as abscissa. 

Figure 16 shows, first, that there is little over-all 
difference in the character of the plots for the two dif- 
ferent ways used to estimate color excesses; except in 
a few individual cases, averages for I and II would be 
satisfactory. Considering, next, the linear diameters 
containing 0.9 the total light, we note they are nearly 
independent of luminosity in the range from Mv:= — 6.8 


to —9.6 and in the diameter range from about 20 to 
50 pe. There is some indication of systematic difference 
between the three groups of clusters: for those with 
|b| > 20°, —20°<b<+20°, and the 10 clusters near the 
galactic center direction, the mean diameters are 
(D)o.9=34, 28, and 22 pc, respectively. While this 
trend may be real, there exists the possibility that it is 
due to increasing field star density, which makes more 
uncertain those diameter determinations for clusters in 
rich star fields. The mean absolute surface brightnesses, 
on the other hand, show a marked correlation with the 
absolute magnitudes, which is to be expected from the 
small range in linear diameter compared to the much 
larger range in luminosity. The intrinsically brightest 
clusters appear to have mean surface brightnesses more 
than 3 mag. brighter than the least luminous clusters. 


XI. CLUSTERS IN THE M31 GROUP 


Magnitudes and Colors 


We consider here the magnitudes V, in this case 
equivalent to V; and the color indices (P-V) and (V-J) 
from Table VII, which are plotted in Figs. 17 and 18. 
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Fic. 16. Linear diameters for 0.9 total light (above) and corresponding mean absolute surface brightnesses (below), for three 
groups of globular clusters, plotted against total absolute visual magnitudes; results are shown separately for alternative I 


(left) and II (right). 


Although this sample is incomplete, especially for the 
fainter objects, it has the advantage of including 79 
clusters all presumably at nearly the same distance. 
Looking, first, at the lower left part of Fig..17, we 
find that there are eight objects with (P-V)<+0.50 
mag. Four of these are in M33 and have already been 
mentioned in Sec. VII. Here we wish to note that the 
brightest one in M33, Mc with V=15.97, is apparently 
more than two magnitudes fainter than the brightest 
cluster in M31, M II with V=13.75. Even with allow- 
ance for as much as 0.5 mag. local reddening in M33, 
it seems probable that the upper limit of luminosity for 
clusters in M33 is significantly lower than that in M31. 
Of the remaining four clusters with (P-V) <+0.50 mag., 
one is H VIII in NGC 205 and the other three are HS, 
H7 and H21 in M31. Although the latter appear in the 
spiral structure, their location in the principal plane is 
necessarily uncertain because of possible projection. 
Object H5, however, was observed in three colors, and 
if plotted in Fig. 10 would be associated with the 
galactic open clusters. For this reason it was observed 
with the Crossley nebular spectrograph. A five-hour 
exposure recorded only a narrow trace without con- 


spicuous absorption lines, and a continuum somewhat 
stronger shortward of H and K than found previously 
for objects M II, III, IV, and VI (Mayall and Egger 
1953). If these three relatively bluer objects in M3! 
are, in fact, like open clusters of higher luminosity, they 
probably are not greatly reddened, for the colors of the 
galactic open clusters listed in Table VI range up tc 
(P-V)=+0.90 mag. 

Considering, next, those objects with (P-V) 2 +-0.5( 
mag. in Fig. 17, we note a fairly high concentration ir 
the approximate range +0.5< (P-V) <+1.0 mag., anc 
a few redder ones that tend to be fainter with increasing 
color index, the most extreme case being B327 witt 
V=16.97 and (P-V)=+1.94 mag. There appears tc 
be, in fact, a fairly well-defined sloping right-hanc 
boundary, as would be found for the highest luminosity 
clusters increasingly dimmed and reddened by inter. 
stellar matter in M31. The effect is indicated in th 
figure by the dashed lines connecting those points cor. 
responding to clusters that are brightest for their colors 
the dotted line shows, with less certainty, the simila1 
effect expected for clusters of high to normal luminosity 

Except for the reddest clusters, it is, of course, im- 
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TABLE XVIII. Distribution of colors of M31 clusters. 
(P-V) { +0.50 0.55 0.60 0.65 0.70 0.75 0.80 0.85 Total (P-V) o 
Interval +0.54 0.59 0.64 0.69 0.74 0.79 0.84 0.89 No. (mag.) (mag.) 
Outlying only 3 5 1 5 1 4 24 +0.68 0.12 
+0.50< (P-V) <+0-90 3 6 5 9 11 4 11 7 56 +0.72 0.10 


possible without further information to identify indi- 
vidually those clusters of small or moderate reddening, 
because of the ranges in true color and brightness. 
Nevertheless, their influence as a group may be avoided 
by considering only those clusters which appear in Fig. 
4 to be sufficiently outlying that little or no space- 
reddening from the spiral structure seems probable. 
There are 24 of the observed objects that may reason- 
ably be placed in this category, and they are indicated 
in Table VII by the superscript a, and in Figs. 17 and 
18 by encircled dots. Their range in color is from 
(P-V)=+0.50 to +0.89 mag. (inclusive), which may 
correspond closely to the true range because the sample 
is of fair size. The frequencies in 0.05 mag. intervals, 
means and dispersions for both the outlying clusters and 
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for all those with +0.50<(P-V)<-+0.90 mag., are 
given in Table XVIII, and their distribution is shown 
as the histogram in the upper part of Fig. 19. 

This distribution in color for these M31 clusters does 
not indicate a well-defined maximum within the (P-V) 
color range of 0.40 mag. Although there may be some 
slight suggestion of subsidiary maxima at about 
(P-V)=0.72 and 0.82 mag., these data are not relevant 
regarding the existence of “blue’- and “ted”-color 
groups of clusters, as suggested by Kron and Gascoigne 
(1952), nor do they indicate a clear separation into 
“halo” and ‘‘disk-type” globular clusters, as identified 
by Morgan (1956). Our data, however, are observa- 
tionally biased toward the apparently brightest objects, 
and the possibility is by no means excluded that obser- 
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Fic. 17. Color-magnitude diagram for star clusters observed in the M31 group, with some probable (dashed) 
and possible (dotted) reddening runs indicated. 
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vations of the many fainter clusters may reveal groups 
like those mentioned. Also, there is little tendency for 
the outlying clusters to be systematically bluer than 
those seen projected over the spiral; in fact, six outlying 
ones have (P-V)>-+0.80 mag. Hence the use of these 
objects, most of which probably are bonafide globular 
star clusters, as indicators of small space reddening by 
interstellar matter either from the Galaxy or in M31 
seems likely to be unsatisfactory, or, at best, attended 
with results of low precision. 


Color Excess and Absorption 


If it is assumed that the eight clusters corresponding 
to those points in Fig. 17 connected by the dashed lines 
are the highest luminosity clusters of nearly the same 
color, (P-V)~0.7 mag., then it is possible to estimate 
the ratio of V absorption to (P-V) color excess in M31. 
A least-squares solution for the relationship between V 
and (P-V) gave the following result with probable 
errors: V=12.11+0.17+ (2.50+0.14)(P-V), which is 
represented by the solid straight line drawn in Fig. 17. 

Thus these M31 data give Ay/E.p_y)=2.50+0.14, a 
value somewhat smaller than that derived by Stebbins 
(1950) from the integrated light of small areas around 
the nuclear region: 3.0 (when reduced to visual). It is 
also smaller than the values usually obtained from 
space-reddened stars in the Galaxy. Hiltner and 
Johnson (1956), for example, found a ratio of 3.040.2 
on their (B-V) color system that has nearly the same 
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color base-line as (P-V). Although a slope of 3.0 would ~ 
not represent so well these few reddest clusters observed — 
in M31, the difference in the ratios, 0.50-L0.24, is too 
poorly determined to suggest seriously that the absorb- 


ing material in M31 differs significantly from that in © 


the Galaxy. Similar conclusions have been reached by © 


Holmberg (1950) for NGC 5195, and by Sérsic (1958) , 


for NGC 5128. 


Color excesses and total absorptions may also be f 


estimated, with low precision, for the few individually — 


reddest clusters represented in Figs. 17 and 18, if a_ 


f 


“normal” or intrinsic color index is assumed for them. — 


Moreover, since some of these were observed in both — 


the (P-V) and (V-Z) color systems, comparison of the 


results provides a useful internal check on the pro-— 


cedure. For this purpose we use, first, the (P-V) data 


for the 12 objects in M31 with (P-V) 2 +0.90 mag., and 
assume that their mean intrinsic color is (P-V)o 


=+0.70 mag., which is the average of the two mean 


values in Table XVIII. For each of these clusters, the 


visual absorption, Ay, was computed from 2.9[(P-V) 
—0.70 ]=2.9 Ew_y)= Ay, where the factor 2.9 from 
Sec. VIII was used because it is probably less subject 
to observational selection error than the factor 2.5 
found from M31. Next, from the three-color plot of 
Fig. 18, it seemed reasonable to obtain an average 
intrinsic (V-Z) color from the 12 clusters represented 
in the rectangle, since among these points are five 
referring to outlying clusters. Thus, for these 12 clusters 


(V-I)o=+0.86 mag., and the corresponding (P-V)o © 
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Fic. 18. Three-color (two color-index) diagram for star clusters in M31. 
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Fic. 19. Comparison of 
_ histograms of colors of star 
_ clusters in M31 and of 
globular clusters in the 
Galaxy; for the latter the 
colors are corrected for 
color excess by alternatives 
Tand IJ. 


50 


=+0.71 mag., which agrees satisfactorily with that 
from the preceding grouping. For the remaining nine 
clusters with (V-Z)>-+1.00 mag., the visual absorption 
was found from 


2.3 (V-1) —0.86 ]=2.3 Ew_1=2.3 Ev_n=Av, 


where the factor 2.3 was obtained as described in Sec. 
VIII. Table XIX gives the results from this procedure. 
The agreement between Ay from the two different 
color systems seems generally satisfactory, except for 
B227. This object is located on the edge of a dark lane, 
and evidently the (V-I) observations (Fig. 18) were 
more affected than the (P-V) ones (Fig. 17). If this 
discrepant (V-Z) case is omitted, the last three columns 
of the table give, respectively, the mean, (A )y, rounded 
to the nearest 0.05 mag., the corresponding (Z)p_y) by 
division by 2.9, and the corrected colors, (P-V),. The 
mean corrected color for these 13 objects is (P-V). 
=0.68 mag., which is 0.03 mag. smaller than (P-V)o 
because of the use of a factor of 2.9 instead of 2.5. 
This treatment of the M31 photometric data indi- 
cates that interstellar matter in the spiral may have 
reddened a few of the observed clusters by about a 
magnitude and dimmed them by several magnitudes. 
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On the whole, however, most of the clusters in this 
sample do not appear to be strongly reddened by local 
absorption, as Baade (1950) had noted previously from 
his color-filter photographs. 


Luminosities 


It isevident from Fig. 17 that the apparently brightest 
clusters in M31 attain high luminosities. For V<14.5, 
an apparent distance modulus of 24.6 (Sandage 1958), 
and without any correction for possible local absorption 
in the spiral, there are four clusters having My<—10; 
namely, M II, B282, H12, and H42. Of these, only M II 
is outlying and, if unreddened, it has My = —10.8, the 
position of H12 in Fig. 17 does not suggest that it is 
reddened, while the locations of B282 and H42 in Fig. 
17 suggest that they are. If the V magnitudes for the 
two latter clusters and for the others in Table XIX are 
corrected by the tabulated values of (A)y, then the 
following 10 clusters have Mv.< —10.0: 


Cluster My. Cluster My. 
H1i2 —10.3 H100 —10.3 
H32 —10.0 B227 —10.8 
H39 —11.2 B282 —11.0 
H42 —11.2 B327 —11.5: 
H93 —11.0 MII —10.8 
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TABLE XIX. Color excess and absorption for M31 clusters. 


Ew-v) Avy Eqn Av 
Cluster (P-V)o= +0770 (V-I)o= +0786 (Ay) (E)p-v) (P-V)e 
H 16 0m22 064 0720 0746 0255) (Ont »e0e7s 
B 282 sie Be Ueera 8 0.25 0.58 0,005, “0.215 S070 
H 38 0523: .067, ses see 0.70 0.24 0.69 
H 100 0.23 0.67 0.43 0.99 0.85 0.29 0.64 
B 259 0.24 0.70 cee see 0.70 0.24 0.70 
42, 0.27 0.78 0.47 1.08 0.95 0.33 0.64 
H 97 0.31 0.90 ao tee 0.90 0.31. 0.70 
B 258 0.32 0.93 ee oes 0.95 0.33. 0.69 
H 93 0.48 1.39 0.47 1.08 1.25." {07435 ROWS 
H 32 0.62 1.80 ste see 1280, ) 105625 COR 70 
H 39 OB hee CARAS) 1.04 2.39 2.30 0.79 0.69 
B 227 0°81. - 2:35 -S(0l61)(1740) 22735 OPS Ore) 
B 327 1.24 3.60 1.85 4.25 3.90 1.34 0.60 


8 Not included because (P-V) = +0.81 <0.90 mag, 


The upper limit of luminosity for clusters in M31 
appears from these results to be of the order of My 
=-—11. The highest luminosity of —11.5 for B327, 
however, is to be regarded as very uncertain, because 
this object was the faintest observed, with P=18.9. 
Also, if the intrinsic colors of any of these clusters are 
redder than the assumed mean value of (P-V),.=0.70 
mag., their listed luminosities will become fainter by 
A(P-V)oX2.9, which amounts to 0.6 mag. if, for ex- 
ample, (P-V)o>=0.90 mag. In order to show the influ- 
ence of these uncertain absorption effects on the 
luminosity frequencies, we give the latter for: (1) out- 
lying clusters only, (2) clusters with +0.50< (P-V) 
<+0.90 mag. uncorrected for absorption, and (3) 
clusters with (P-V)>+0.50 mag. including those 
corrected for absorption as described in the foregoing. 
Table XX contains these three sets of frequencies in 
0.5 mag. intervals, and the distribution for (1) and (2) 
is shown graphically as a histogram in the upper right 
part of Fig. 20; the upper left part gives the histogram 
of frequencies for P-magnitudes. 

These frequencies suggest that the luminosities of 
these M31 clusters have a well-defined maximum around 
My=—9, with a moderate dispersion of c=0.6 mag. 
if the highly space-reddened ones are disregarded. This 
result, however, could be considerably affected by ob- 
servational selection, because the observing program 
included mainly those objects of photographic mag- 
nitude 16.5 and brighter. Thus the frequencies for 
My=-—8 and fainter are hardly represented in the 
sample. On the other hand, the 24 outlying clusters, 
whose census probably is almost complete, give nearly 
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the same distribution. Obviously, observations of more 
and fainter clusters are needed to establish more reliably 
the fainter part of the luminosity function. 


XII. COMPARISON OF GLOBULAR CLUSTERS IN 
THE GALAXY AND IN M31 


Colors 


In the left-hand side of Fig. 19, the scales of abscissae, ( 
or (P-V) colors, have been aligned vertically in order to 
compare the color frequency distributions of globular ~ 
clusters in M31 and in the Galaxy. In the histogram for — 
M31, only those clusters are included that have - 
0.50< (P-V)<0.90 mag., and in that for the Galaxy 
the (P-V). colors are those corrected on the basis of 
alternate I. These two histograms indicate that in M31 
a substantial majority of the clusters have colors in the 
range from (P-V)=0.65 to 0.90 mag., whereas in the 
Galaxy a large proportion of the globular clusters have 
colors in the range from (P-V)=0.45 to 0.65 mag. 
Whether this systematic difference is real, and not due 
to observational selection of clusters in M31 and to 
incorrect color excesses for clusters in the Galaxy, prob- 
ably remains to be determined by future work. It 
does not seem probable to us, however, that reasonably 
different procedures for deriving color excesses for the 
galactic globular clusters would produce a distribution 
of colors comparable to that in M31. Many of the 
galactic globulars with |b|>20° are in the range of 
(P-V)-=0:45 to 0.65, and they are so little reddened 
that it does not make too much difference how their 
color excesses are estimated. Also, although most of the 
clusters with (P-V),=0.65 to 0.90 mag. are located 
around the center, the observing program included a 
large fraction of the known number, and observations of 
additional ones could hardly produce a frequency maxi- 
mum comparable or higher than that in the bluer range. 
If alternate II is used to estimate color excesses, there is 
even less possibility of getting a maximum in the redder 
range. As for the clusters in M31, since uncorrected 
colors were used in the histogram, we cannot rule out the 
possibility that many of those seen in projection over 
the spiral may be reddened by amounts averaging 0.25 
mag. in (P-V). Nevertheless, a substantial percentage 
would have to be reddened by this amount in order to 
produce a frequency maximum in the range from 
(P-V),=0.45 to 0.65 mag. This condition seems im- 
probable if half are located in front of the principal 
plane containing the absorbing material, and half in 
back. Possibly a more promising way to get a higher 


TABLE XX. Luminosity distribution of M31 clusters. 
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Fic. 20. Comparison of the histograms of absolute photographic (left) and visual (right) magnitudes of star clusters in M31 
(above) and of globular clusters in the Galaxy (below), computed with the assumption that Mp=0.0 for cluster-type (RR Lyr) 


variables. 


maximum in the bluer range would be to observe the 
remaining two-thirds of the clusters known in M31. 
Since nearly all of these are fainter than photographic 
magnitude 16.5, however, it is also probable that many 
will be considerably reddened, and some will be open 
clusters, with the result that inference of their intrinsic 
colors will be a difficult and uncertain procedure. Until 
such additional data are available, we conclude that the 
present results give some support for a systematically 
different color frequency distribution, in the sense that 
the clusters in M31 may be intrinsically redder, by 0.2 
to 0.3 mag. on the average, than the globular clusters 
in the Galaxy. 


Luminosities 


In Fig. 20, the left-hand side shows the frequency 
distributions in photographic light, and the right-hand 
side in visual light, of the apparent magnitudes of the 
clusters in M31 (above) and of the absolute magnitudes 
of globular clusters in the Galaxy (below). In this 
figure the maximum frequencies are vertically aligned, 
with the result that the scale of abscissae have a 
constant difference of 23.6 mag. in photographic light 


and 24.1 in visual light. The bearing of these figures on 
the distance modulus of M31 is considered in the next 
section, which is a discussion of globular clusters as 
distance indicators. 

Considering, first, the brighter side of the luminosity 
distributions represented in Fig. 20, we note that the 
histograms for the M31 clusters appear to fall off faster 
and to extend to higher luminosities, than the ones for 
the galactic globular clusters. This characteristic, if 
real, probably is due to the fact that for M31 this 
program included all but the most heavily absorbed 
high luminosity clusters, whereas for the Galaxy the 
intrinsically brightest ones in the southern hemisphere 
are not included. [Note added in proof. We regret having 
overlooked, when writing this section, the important 
paper by S. C. B. Gascoigne and E. J. Burr (Monthly 
Notices Roy. Astron. Soc. 116, 570, 1956) on 47 Tucanae 
and w Centauri. From photoelectric photometry they 
found Mp,=—10.0 for 47 Tuc and —10.4 for w Cen. 
If these values were included in the lower left histogram 
of Fig. 17, the agreement between the luminosity dis- 
tributions for the Galaxy and M31 would be much 
improved. | This kind of observational selection argues 
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against making a too-detailed comparison of the high- 
luminosity “tails”; about all we are willing to con- 
clude is that the present data do not indicate any 
substantial difference in this respect between M31 and 
the Galaxy. Second, we note that all the histograms 
show fairly well-defined maxima, but, as stated in Sec. 
XI, those in M31 probably are misleading, if not 
spurious, because of the lack of observations for the 
great bulk of clusters fainter than photographic mag- 
nitude 16.5. For the Galaxy, on the other hand, the 
frequency distributions are reasonably complete (about 
60%), and as significant as permitted by our present 
knowledge of distance moduli, and by the assumption 
of Mp=0.0 for cluster-type (RR Lyrae) variables. 


XIII. GLOBULAR CLUSTERS AS DISTANCE INDICATORS 


In this section we consider how the photometric data 
obtained in the present work may be used to estimate 
the distance moduli of the Andromeda nebula (M31), 
Magellanic Clouds, and the galactic center. This was 
one of the original aims of the program, which, if suc- 
cessful, could be extended with larger telescopes to 
other, more distant galaxies. A few of these have been 
resolved into outlying objects whose colors, distribution, 
order of luminosity, and possible nonstellar character 
suggest they probably are globular clusters (Bowen 
1952, Baum 1955). Since they are among the highest- 
lacunosity. objects in the systems, their importance as 
distance indicators is obvious. 


Andromeda Nebula (M31) 


For this spiral we use the results shown in Fig. 20, 
in which the scales of abscissae have been adjusted to 
make the maximum frequencies coincide: In photo- 
graphic light (left), the correspondence gives P—Mp; 
= 23.6 mag. when the filled circles are used. If only the 
open circles were used, which means only outlying 
objects in M31 and high-latitude clusters in the 
Galaxy, then a lining up of maxima would give 
P—Mr,= 24.1 mag. This latter figure is undoubtedly 
preferable, because in M31 the observations are com- 
plete only to P~16.5 mag., while in the Galaxy the 
distances of the high-latitude clusters are best deter- 
mined. In visual light (right), the maxima given by the 
filled circles yield V—Mv,= 24.1 mag., while the open- 
circle maxima could give either 24.1 or 24.6 mag., 
depending on how the sharper maximum for M31 is 
placed with respect to the broader one for the Galaxy. 

Thus the foregoing comparison of apparent magni- 
tudes of star clusters in M31 with absolute magnitudes 
of globular clusters in the Galaxy gives a difference 
ranging from 23.6 to 24.6 mag. If these objects are 
comparable in luminosity, then a reasonably conser- 
vative estimate of the M31 distance modulus is 
(m—M)=24.10.5 mag., based on the premise that 
Mp=0.0 for cluster-type (RR Lyrae) variables, and 
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with no correction for galactic absorption in the direc- 
tion of M31. Although in the past estimates of ee | 
correction have ranged up to 0.7 mag., there now seems | 
to be little justification for them, because of the work \ 
of Code and Houck (1956). Prout three-color observa- i 
tions, they found that the colors of .several of the } 
brightest stars in an outer spiral arm, distant 75’ from‘ 
the nucleus, are bluer than the galactic OB stars. On, i 
the other hand, there may be reason to revise the, } 
deduced modulus downward by 0.5 to 0.6 mag., } 
because of recent work on the luminosities of cluster- | 
type variables and of RR Lyrae (Pavlovskaya 1953, | 
1954; Arp 1958; Sandage and Eggen 1959; Baum ef al. | 
1959); if so, then our result would become (m—M) |} 
~23.50.5 mag. However, since the luminosities of | 
cluster-type variables in the high-latitude or halo-type | 
globular clusters may be closer to Mp=-+0.2 mag. than jf 
to an uncertain mean in the range of Mp from about | 
+0.2 to +1.0 mag., we conclude that the M31 distance 

modulus, as estimated from its probable globular } 
clusters, falls within the range (m—M)=23.5 to 24.0: } 
mag., with an uncertainty of 0.5 mag. 


Magellanic Clouds 


Since only five probable globular clusters were | 
observed in the Small Cloud and six (including the | 
poor-data case’ of NGC 2121) in the Large Cloud, we | 
cannot construct any significant luminosity distribu- 
tions for estimation of separate distance moduli for the | 
two systems. This procedure, however, has been used _ 
for 97 SMC clusters by de Vaucouleurs (1959), who | 
employed small-scale photographic photometry, with — 
Gascoigne and Kron’s (1952) photoelectric cluster data 
as standards. Instead, we note that the mean, uncor- 
rected apparent total magnitudes are (P;)=11.6 or — 
10.9, and (V,)=11.0 or 10.2 mag., for the Small or / 
Large Cloud, respectively. Thus there is an average 
systematic difference of 0.7 or 0.8 mag. between these 
Cloud clusters, with those in the SMC being fainter. 
This difference has already been noted before (Kron 
1956), and also seems to be confirmed by de Vaucouleurs’ 
photographic results. Nevertheless, the photoelectric — 
data are presently too few to place such a difference 
beyond question, so in this discussion we shall use the © 
over-all means of (P),;=11.2 and (V),=10.6 mag. We — 
next review some of the data bearing on correction of 
these values for galactic and local Cloud absorption. 

Until a few years ago, it was customary to use galactic 
absorption corrections of A,,=0.3 or 0.4 mag., com- ~ 
puted either (Shapley 1951, de Vaucouleurs 1954) from 
Harvard galaxy counts fitted to a cosecant law, or 
Apy=0.4 to 0.6 mag. calculated (Thackeray and 
Wesselink 1953) from Oort’s (1938) discussion of 
Hubble’s (1934) data. Lack of more suitable photometric 
results left no alternative to this procedure, which was 
known to be unsatisfactory because of the considerable 
unevenness in the background galaxy population: poor 
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TABLE XXI. Comparison of results for galactic globular clusters. 


NGC Apert. Vs Vim J-KM (P-V); J-KM NGC Apert Vs Vem J-KM (P-V); J-KM 
2419 138” 11705 10797 +0708 +0"57  —0=07 6553. 106’" = (9™37) =9™44 (—0"™07) +1750 +0™11 
5904 = 250 6.67 6.45 +0.22 0.66 +0.03 6569 250 (9.10) 8.90 (+0.20) 1.24 +0.13 
6205 49 a8 oe +s: 0.56 —0.01 6624 250 8.51 8.45 -+0.06 0:99" 7 =0-01 
6205 82 7.85 8.08 —0.13 0.58 +0.01 6626 250 7.49 7.29 --0,20 P05 55-0208 
6205 250 6.54 6.46 +0.08 0.60 -+0.03 6638 106 OS) D560) OE 1.08 --0.12 
6218 250 7.95 7.70 +0.25 0.79 +0.02 6656 250 6.38 67595 9-0; 08 0.99, --0:13 
6254 250 deol 7.39  +0.12 0.88 -+0.09 6809 82 9.78 ONO07 ate 0.60 +0.05 
6273 250 i393 7.26 +0.07 0°98 +-0:.08 6809 250 7.62 7.43 -+0.19 0566") =-0. 11 
6293 250 8.77 8.65 +0.12 0.92 -+0.06 7006 82 11.08 11.08 0.00 0.66 +0.05 
6304 250 (8.77) 8.82 (—0.05) 1.22 0.00 7006 250 LOKOS O60 Oot O77 420216 
6316 250 (9.25) 9.36 (—0.11) 1.10 —0.07 1078 82 Looe 7.38 +0.14 0.62 +0.03 
6341 250 6.97 6:92, -£0.05 0.53 0.00 7078 250 6.88 6-15) 2 --On15 0.62 +0.03 
, 6402 106 (8.86) 9.09 (—0.23) 1.27 0.15 7089 82 7.63 Leon = OU 05620 0205 
6402 250 (8.20) 8.06 (+0.14) 1.20 +0.08 7089 106 7.41 7.48  —0.07 O502 000 
6517 49 (11.80) 11.75 (40:05) (1.71) (+0.04) 7089 250 6.88 6.7/6 --0.12 OF 02.5 0705 
6522 49 9.99 9762 -+0.37 1.08 -+0.06 7099 82 8.65 S019) --0,02 0.51 -F0.03 
6522 250 8.75 8.71 +0.04 +1.04 +0.02 7089 =. 250 7.95 7h, 80\nera|- O09 e EOS SS ane OuUS 


| around the Small Cloud and rich around the Large 


(Shapley 1957). The question has more recently been 
considered by Feast, Thackeray and Wesselink (1958) 
and by Arp (1958), on the basis of photoelectric obser- 
vations of colors of stars in the galactic foreground and 
in the Clouds, but they obtained somewhat different 


answers. The Pretoria workers found A,,=0.2 to 0.8 


mag., depending on assumptions made for the intrinsic 
colors of the stars observed. They concluded that ‘‘the 
early-type supergiants studied in the two clouds suffer 
a small but definite amount of absorption, partly due 
to galactic dust and partly due to dust within the 


_ Clouds themselves,” but they refrained from quoting 


any average values. Arp, on the other hand, in his 
extensive program of southern hemisphere photometry 
at the Cape, found that some faint, early-type galactic 
stars yielded a reddening of +0.01+0.01 mag. in the 
SMC field, and +0.04+0.02 mag. in the LMC field. 
He obtained significant reddening, among the stars 
observed, for only two of six supergiants identified in 
the SMC. Thus he confirmed local absorption within 
the Clouds, but not through the Galaxy in the direc- 
tions to them. 

The foregoing review suggests that there may be 
little, if any, galactic absorption towards the Clouds, 


but that some objects within them are locally reddened 


and obscured. Of the Cloud globular clusters observed 


_ photoelectrically in this program, nearly all are outlying 


or fringe objects. Only one cluster, NGC 1835 at one 
end of the central bar in the LMC, is an interior object, 
but its color, (P-V)=+0.61 mag., seems normal. Thus 
these 11 clusters probably are negligibly locally reddened 
in the Clouds. We shall, therefore, apply no corrections 
for galactic or Cloud absorption to the mean apparent 
magnitudes, (P),=11.2 and (V),=10.6 mag., of these 
Cloud clusters. If they are of the same luminosity as 
those in the Galaxy, for which we found (Sec. X, Table 


~ XV) (M)p,=—7.7 and (M)v.=—8.2 mag., then the 


estimated distance moduli from the two colors are 
(m—M)p=18.9 and (m—M)y=18.8 mag. Since first- 
decimal accuracy can hardly be expected in these values, 
we prefer to round off the mean to 19, and to estimate 
its uncertainty at 0.5 mag. 


Galactic Center 


We estimate this distance in two ways: (1) from the 
group of globular clusters seen in the general direction 
of the center, and (2) from the central bulge RR Lyrae 
variables studied by Baade (1951, 1953) in the field of 
NGC 6522. 

1. With the omission of the two poor-data cases of 
NGC 6401 and 6453, there are 23 clusters, identified by 
the superscript letter b in Table VI, which appear to 
form a fairly compact group around the direction to the 
galactic center. For this group the mean corrected total 
apparent visual magnitude, obtained from the indi- 
vidual values of V;, in Table XIII, is (V);.=7.5 or 
7.1 mag., for color-excess alternates I or II, respec- 
tively. From Table XV and the discussion in Sec. X, 
it seems reasonable to assume for this center group 
that (M)v,=—8.2 mag., which gives distance moduli 
of m—M=15.7 (1) or 15.3 (II) mag., and corresponding 
distances of 13.8 or 11.5 kpc. How uncertain these 
figures may be is difficult to determine. Their formal 
probable errors, if computed from those of (V):. and 
(M)v1, are of the order of +0.1 to +0.2 mag., which 
are unrealistic because they do not include the more 
serious uncertainties in zero-point and sampling. Some 
indication of the latter may be obtained from the 10 
clusters in the group for which individual corrected 
moduli are available (Table XVII). The corresponding 
mean distances (and ranges) are (r)=13.0(7.2 to 20) 
or 10.9(6.3 to 16) kpc, for color excess alternates I or 
II, respectively. With such large ranges, the group 
means are probably not very trustworthy, although 
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their formal probable errors are only +1.0 I and 
-+0.7 II kpc. We conclude from these results that this 
group of clusters indicates a galactic center distance of 
the order of 12.5 kpc, with an estimated uncertainty 
of 1.5 kpc. 

2. From photographic photometry based on photo- 
electric standards (Stebbins, Whitford, and Johnson 
1950), Baade (1951) found that the most frequent 
apparent magnitude of a large number of RR Lyrae 
variables in a field at 1=328°, b=—4° is my,=17.5 
mag. From the fact that the frequency maximum is 
steep and of small spread, he concluded that these 
variables are located in the galactic nuclear bulge. 
Thus they could be used to estimate the distance to the 
galactic center, with certain assumptions. The first is 
that M,,=0.0 for these variables, the second that their 
space absorption is the same as that of the globular 
cluster NGC 6522 in the same field. Using Stebbins and 
Whitford’s unpublished photoelectric color measure- 
ments of this cluster. Baade derived an absorption of 
A p= 2.8 mag. A few years later, following further work, 
Baade (1953) presented revised values of 17.3 and 2.75 
mag., or m—M=14.6 mag. and a galactic center 
distance of 8.16 kpc, which has subsequently been widely 
used. Since the crucial point is how the absorption 
correction is obtained from the observed color of NGC 
6522, we shall make an estimate based on our data for 
this cluster. From Table X the color excess for NGC 
6522 is (E)~p_y)=0.48 mag., which multiplied by 4 
gives A »,=1.92 mag. Thus for M,,=0.0 for RR Lyrae 
variables, the corrected distance modulus is m—M 
=17.3—1.9=15.4 mag., or a distance of 12.0 kpc. 
Again, it is difficult to assess the accuracy of this result 
by a formal probable error calculation, but if the color 
excess error is +0.1 mag., then the modulus i is uncertain 
by at least 0.4 mag. and the distance by +1. 2 kpe. 

Although these two ways of estimating the galactic 
center distance give results accordant to better than 
1 kpc, such close agreement is undoubtedly accidental 
because of the several shaky assumptions and sys- 
tematic errors involved. Nevertheless, our results lead 
to a galactic center distance nearly 50% greater than 
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that in current use, unless the zero-point for luminosities” 
of RR Lyrae variables is considerably revised. Thus, 
to obtain a galactic center distance as small as 8.2 kpe 
from our data and procedure, it would be necessary to. 
use Mp~+0.9 or +0.8 mag. for RR Lyrae variables. 
This revision has already been suggested, as mentioned 
in the above discussion of the M31 modulus, and future 
work may establish it. If so, then the different 
photographic space-absorption correction for NGG 
6522 found by Baade, 2.8 mag., and by us, 1.9 mag., 
will serve as an illustration of how difficult i is to get. 
reliable color excesses and absorptions for obscured 
galactic center objects, especially when they are of 
composite character and of uncertain intrinsic color 
(Whitford 1960). 


ADDENDUM 


During preparation of this paper for publication, 
some new photoelectric photometry of star clusters 
was reported: by Johnson (1959) for 27 galactic 
globular clusters, and by Hiltner (1960), for 46 star 
clusters in M31 and M33—23 in each system. In both 
cases the observations were made on the U, B, V 
system. After conversion of their results to our P, V 
system according to the relationships given in Sec. III, 
we have the comparison for objects in common as shown 
in Tables XXI and XXII. 

In Table XXI the columns headed Vy; and (P-V)5 
contain magnitudes and colors by Johnson reduced to 
our P, V system; those in parentheses represent con- 
versions in the color range where the transfer formulae 
are not well determined. The column headed V xm gives 
our magnitudes for the apertures used by Johnson, 
listed in the second column. For this purpose we used 
the same luminosity distribution curves that gave the 
total magnitudes and diameters in Table IV. The two 
columns of differences, in the sense J-KM, yield average 
systematic differences and AD’s, respectively, of +0.06 
and 0.10 mag. in the magnitudes, and +0.05 and 0.04 
mag. in the colors. The largest difference in the inte- 
grated magnitudes is +0.37 mag. for NGC 6522 meas- 
ured with the 49” aperture; however, this difference is 


TABLE XXII. Comparison of results for extragalactic star clusters. 


H-KM 


Cluster Via (P-V)x  H-KM Cluster Va H-KM (P-V)a H-KM 

ve ye, 

Tele, 16™27 —0=10 +031 +0™20 H 94 15™88 —0™02 +089 +0704 
26 15.54 —0.03 0.89 +0.18 99 15.74 +0.08 0.69 +0.04 
Pai 15.68 0.00 0.76 +0.05 106 dS AZ, +0.12 0.75 —0.06 
42 14.33 +0.01 1.08 +0.11 11D 5 43 —0.04 0.66 +0.03 
44 15.10 +0.14 0.88 +0.55 115 15e57 —0.01 0.93 +0.11 
55 15.60 —0.04 0.90 +0.13 116 15.46 +0.05 0.91 +0.05 
57 16.05 +0.08 0.90 +0.05 B 231 15.30 +0.04 0.82 +0.07 
64 14.77 +0.04 0.84 +0.03 258 15.80 —0.02 Lea? +0.15 
73 15.50 +0.09 0.90 +0.10 282 14.25 +0.04 +0.88 +0.07 
76 15.68 +0.03 +0.56 +0.06 

M33 
Ha=e 16.21 +0.02 +0.70 +0.33 Hb=d 16.67 -+0.13 +0.11 +0.09 
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only +0.04 mag. for the 250” aperture. All the other 
integrated magnitude differences fall in the satisfac- 
torily small range from —0.23 to +0.25 mag. On the 
other hand, the differences in integrated color range 
from —0.07 to +0.16 mag., which, although not so 
good, may serve as an indication of the difficulties, and 
hence of the lack of high precision, of integrated-light 
and large-aperture photometry in rich star fields of 
nonuniform density. 

In the discussion of his data on the basis of the same 
intrinsic color, (B-V)=-+0.60 or (P-V)=+0.52 mag., 
for all globular star clusters, Johnson found for NGC 
6522 a maximum color excess of E,z_y)=+0.53 mag. 
and a total blue absorption of A4g=+2.1 mag.; this 
value agrees satisfactorily with ours of A,,=+1.9 mag., 
as deduced in Sec. XIII. Johnson also compared his 
estimate with Baade’s (1958), but the comparison was 
made with the misprinted value of +2.15 on page 329, 
instead of with the correct number +2.75 on page 61 
of the cited reference. 

In Table XXII the columns headed Vy and (P-V) x 
contain magnitudes and colors determined by Hiltner 
and reduced to our P, V system. For M31 the average 
systematic differences and AD’s are, respectively, 
+0.02 and 0.05 mag. in the magnitudes, and +0.08 
and 0.08 mag. in the colors. The magnitudes thus agree 
very well, but not the colors, which seem to have an 
appreciable systematic error. If 0.08 mag. is subtracted 
from the differences, the AD becomes 0.05 mag., which 
would be satisfactory. Since we cannot account for 
such a systematic difference, however, we do not wish 
to consider it as a correction to either series of obser- 
vation; instead, we should be inclined to take averages 
from the two sets. For M33 only two clusters were 
observed in common, and the difference of +-0.33 mag. 
between the measurements of color for Ha= Hiltner’s e, 
is especially large. Nevertheless, Hiltner’s much more 
extensive material for M33 leaves little doubt that its 
star clusters are systematically bluer and fainter than 
most of those observed in M31, as we had concluded 
previously (Secs. VII and XI) from fewer data. This 
result thus emphasizes again that there are significantly 
different types of star clusters in different types of 
of galaxies. 

Lastly, Sandage and Wallerstein (1960) have pub- 
lished their detailed study of NGC 6356, a cluster in 
Morgan’s (1959) spectroscopic group having moderately 
strong metallic lines. From three-color photoelectric 
observations and from spectrograms of 12 early-type 
stars in the same field, they found_color excesses ranging 
from E,z_vy)=0.0 to 0.9 mag. Using these, they adopted 
a value of E,z_y)=0.5-+40.1 mag. for the cluster. Their 
estimate may be compared with our alternate ones of 
Evp_vy)=0.28 (I) or 0.48 (II); the agreement seems 
about as good as may be expected, considering the dif- 
ferent assumptions involved. For RR Lyrae My=0.0 
they also derived an apparent modulus of (m—M)y 
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=17.7, compared to our mean of 16.8 (Table XIII). 
Since the latter modulus depends upon Dr. Helen S. 
Hogg’s (1953) preliminary photographic measurements 
of brightest stars in NGC 6356, the photoelectric data 
indicate a correction of +0.9 mag. to the photographic. 
Accordingly, our corrected moduli and distances (Table 
XVII) for this cluster should be revised to: (m—M), 
=17.0 (I) or 16.2 (II), and r(kpc)=25 (I) or 17 (ID), 
compared to Sandage and Wallerstein’s estimates of 


.(m—M),=16.2 and r(kpc)=17.4. On the ad hoc as- 


sumption that NGC 6356 is close to the galactic 
nucleus, they suggested that their values obtained from 
RR Lyrae My=-+0.9, (m—M),=15.3 andr (kpc) =11.5, 
are to be preferred. 


ACKNOWLEDGMENT 


We are indebted to S. C. B. Gascoigne and to Joel 
Stebbins for permission to make use of the photo- 
metric and colorimetric data of the open clusters. These 
data originated in a joint project undertaken by 
Gascoigne and Kron at the Lick Observatory in 1956; 
much of the observing, however, was done by Stebbins 
and Kron. 


REFERENCES 


Arp, H. C. 1955, Astron. J. 60, 317. 

. 1958, ibid. 63, 118. 

Fear. ae Handbuch der Physik (Springer-Verlag, Berlin), Vol. 

»P. fd. 

Baade, W. 1944, Astrophys. J. 100, 137. 

—. W. 1951, Publs. Univ. Michigan Obs. 10, 16. 

. 1953, Symposium on Astrophysics (University of Michigan, 

Ann Arbor, Michigan), p. 25. 

. 1958, Le Probleme de Populations Stellaires (Pontificiae 
Academiae Scientiarum Scripta Varia, Vatican City), Vol. 16. 

Baade, W. and Hubble, E. 1939, Publs. Astron. Soc. Pacific 51, 40. 

Baade, W. and Mayall, N. U. 1950, Problems of Cosmical Aero- 
dynamics (Central Air Documents Office, Dayton, Ohio), p. 166. 

Baum, W. A. 1955, Publs. Astron. Soc. Pacific 67, 328. 

. 1959, Astrophys. J. 130, 749. 

Becker, W. 1938, Z. Astrophys. 15, 225. 

Bowen, I. S. 1952, Carnegie Inst. Washington Yearbook 51, 19. 

. 1956, ibid. 55, 48. 

Christie, W. H. 1940, Astrophys. J. 91, 8. 

Code, A. D. and Houck, T. E. Astron. J. 61, 173. 

de Vaucouleurs, G. 1954, Australian J. Sci. Suppl. 17, 7. 

——. 1959, Publs. Astron. Soc. Pacific 71, 202. 

Eggen, O. J. 1951, Astrophys. J. 114, 141. 

——. 1955, Astron. J. 60, 65. 

Feast, M. W., Thackeray, A. D., and Wesselink, A. J. 1958, Ob- 
servatory 78, 156. 

Gascoigne, S. C. B. and Kron, G. E. 1952, Publs. Astron. Soc. 
Pacific 64, 196. 

Hiltner, W. A. 1960, Astrophys. J. 131, 163. 

Hogg, H. S. 1953, David Dunlap Obs., Communication No. 34. 

. 1959, Handbuch der Physik (Springer-Verlag, Berlin), Vol. 
53, p. 166. 

acibeee E. 1950, Lund Obs. Meddelanden, Ser. II, No. 170. 

Hubble, E. 1932, Astrophys. J. 76, 44. 

. 1934, ibid. 79, 8. 

Johnson, H. L. 1952, Astrophys. J. 116, 272. 

. 1953, ibid. 117, 313. 

Johnson, H. L. and Morgan, W. W. 1951, Astrophys. J. 114, 522. 

Keenan, P. C. 1940, Astrophys. J. 91, 113. 

Kron, G. E. 1956, ibid. 68, 230. 

. 1958, ibid. 70, 285. 

Kron, G. E., Greeby, R. W. and Willson, J. R. 1956, Publs. 
Astron. Soc. Pacific 68, 544. 


620 Gs. MES. 


Kron, G. E. and Smith, J. L. 1951, Astrophys. J. 113, 324. 

Kron, G. E., White, H. Si and Gascoigne, S.C. B. 1953, Astrophys. 
J. 118, 502. 

Mayall, N. U. 1946, Astrophys. J. 104, 290. 

Mayall, N. U. and Eggen, O. J. 1953, ” Publs. Astron. Soc. Pacific 
65, 24. 

Morgan, W. W. 1956, Publs. Astron. Soc. Pacific 68, 509. 

——. 1959, Astron. Je 64, 432. 

Morgan, W. W.., Harris, D: L., and Johnson, H. L. 1953, Astro- 
phys. J. 118, 92, 

Mowbray, A. G. 1946, Astrophys. J. 104, 47. 

Nassau, J. J. and Mac Rae, D. A. 1955, ’ Astrophys. Js A2ZhS2: 

Nassau, J. J. and Seyfert, C. K. 1945, Astrophys. J. 102, Bit, 

Oort, J. 1938, Bull. Astron. Inst. Neth. 8, 233. 

Pavlovskaya, P. 1953, Variable Stars (U.S.S.R.) 9, 233; 1954, 9, 
349. 

Sandage, A. R. 1958, Astrophys. J. 127, 513. 

Sandage, A. R. and Eggen, O. J. 1959, Monthly Notices Roy. 
Astron. Soc. 119, 255. 


KRON AND WN. UleMAYARL 


Sandage, A. R. and Wallerstein, G. 1960, Astrophys. J. 131, 598. 

Sawyer, H. B. 1947, David Dunlap Obs. Publ. 1, 383. d 

Seares, F. H. and Joyner, M. C. 1945, Astrophys. J. 101, 15. 

Sersic, J. L. 1958, Observatory 78, 24. 

Shane, C. D. and Wirtanen, C. A. 1954, Astron. J. 59, 285. 

asa H. 1951, Publs. Univ. Michigan Obs. 10, 79. 

—. 1953, Proc. "Natl. Acad. Sci. U.S. 39, 358. 4 

oS : 1957, The Inner Metagalaxy (Yale University Press, New 
Haven, Connecticut) Pp. 97-100. q 

ete H. and Sayer, A. R. 1935, Proc. Natl. Acad. Sct. U. S.4 


Stebbins, J. 1933, Proc. Natl. Acad. Sci. U. S. 19, 222. 
Stebbins, J. and Whitford,-A. E. 1936, Astrophys. J. 84, 132. 
==: 1937, ibid. 86, 247. 
—. - 1945; ibid. 102, 318. y 
Stebbins, Ea Whitford, A. E. and Johnson, H. L. 1950, Astrophys. — 
J. 112, 469, 7 
Thackeray, A. D. and Wesselink, A. J. 1953, Nature 171, 693. 
Whitford, A. E. 1960, Pubdls. Astron. Soc. Pacific (to be published). 


THE ASTRONOMICAL JOURNAL 


VOLUME 65, 


NUMBER 10 DECEMBER, 1960 


Effect of Precession and Nutation on the Orbital Elements of a 
Close Earth Satellite 


YOSHIHIDE Kozat 
Smithsonian Astrophysical Observaiory and Harvard College Observatory, Cambridge, Massachusetts 


(Received August 26, 1960) 


Perturbations due to the motion of the equatorial plane of the earth are derived for the orbital elements 
of a close earth satellite. It is suggested that, for precise studies of satellite motion, a system be adopted in 
which the inclination and the argument of perigee are referred to the equator of date, and the longitude of 
the node is measured from a fixed point along a fixed plane and then along the equator of date. 


INTRODUCTION 


a a previous paper (Kozai 1959b) we studied the 
_# motion of an artificial satellite by assuming that 
the equator of the earth is fixed to an inertial system 
of reference. 
However, the earth is moving around the sun, and 
because of precession and nutation the equatorial 
plane is not fixed in space. The effect of the orbital 
motion of the earth can be taken into account as 
direct perturbations due to the sun and the moon 
(Kozai 1959a). The effects of precession and nutation, 
which are studied in the present paper, may be regarded 
as indirect luni-solar perturbations. 


DISTURBING FUNCTION 


The declination 6 appearing in the expression for the 
gravitational potential of the earth, 


GM As 
y=—| 1+ sints) + : | (1) 


Tt 7; 


is measured from the instantaneous equator of the 
earth. However, the equator of the earth undergoes 
precessional and nutational motions in inertial space. 
Therefore, in solving the satellite equations of motion 
in inertial space, several time-dependent terms are 
introduced in the expression for the earth’s potential. 
It is the purpose of this paper to examine the dynamical 
terms introduced by these time-dependent terms. 
- The equator at a certain initial epoch is adopted as 
the fixed plane. The argument of latitude, the argu- 
ment of perigee, the longitude of the ascending node, 
and the inclination with respect to this fixed plane are 
denoted by L, ¥, N, and J, and those with respect to 
the equator of date are denoted respectively by u, , 
Q, and z (Fig. 1). The longitudes of the node NV and Q 
are measured from the fixed equinox £ at the initial 
epoch. In Fig. 1, two angles 9 and Vo represent the 
inclination and the longitude of the ascending node 
of the equator of date referred to the fixed plane. 

The mutual inclination of the two equators @ is 
assumed to be a very small angle, so that its square &” 
can be neglected in the following discussion. Then, 

_ between the two sets of elements the following relations 


hold: 


J—1i=8 cos(N—WN), 
N—Q=—6 cotJ sin(V—N>), (2) 
Yy—w=86 cscJ sin(V—N>). 


The declination 6 is expressed by the elements referred 
to the fixed plane as follows: 


sinéd=sinL| sinJ—6 cosJ cos(N—N) | 
—6cosL sin(N—WN). 


Introducing this relation into Eq. (1) leads to a 
potential function which is not symmetric about the 
equator of the fixed coordinate system. 

In relation to the short-periodic perturbations due 
to the oblateness of the earth, the effect of the motion 
of the equator is very small and can be neglected. The 
most important effect comes from the secular part of 
the disturbing function. This is given by 


A,GM 
R=— (1-2-4 sin’) 
ie (02 


+26 sin2J cos(V—WNo) |. (3) 


The terms of the first part of R are of the same form 
as those given for the secular part in the previous 
paper (Kozai 1959b) except for a difference in the 
definitions of the inclination. The term in the second 
part of R arises from the fact that the axis of symmetry 
of the potential field deviates from the pole of the 
fixed coordinate system by the angle 0. 


Orbital plane 


Equator of ddte 


Fixed equator 


Fic. 1. Equators and orbit. 
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PERTURBATIONS 


Equations expressing the variations of the orbital 
elements with the disturbing function R show that the 
semi-major axis a and the eccentricity e are constant 
and, of course, they have the same values in the two 
systems of coordinates. 

The variations of the other four elements produced 
by the term proportional to 6 in R are then given by 
the differential equations 
d:-6J Az 


dt P 
d-6N Ao» fF cos2/ 
nN [ 


né cosJ sin(N—No), 


cos(V—No)-tsin/- ar 


he ie sind 


dbp As cosJ cos2J 
= nfo 


—$sin2J ) cos(V— No) 
dp 


(4) 


sinJ 
+5 sinJ cons | 

d:6M 3A, 
= n(1 


dt 2p 


e)} sin2J[@ cos(N—No)—6J ], 


where 8J in the right-hand sides is derived from the 
first of these equations, and 


p=a(l—e). 


Now let us denote obliquities of the fixed and 
moving equators to the ecliptic by ¢ and ¢’, respectively, 
and the angular distance between the two equators 
measured along the ecliptic by o. The ecliptic may be 
assumed to be fixed over several years. As angles (i) we 
and e—é’ are very small, the following relations hold: 


6 sinNo=c sine’, _ 
(5) 


6 cosNo= €’—e. 


The principal terms whose coefficients are larger than 
one tenth of a second of arc in the expressions of 
o sine’ and e—é’ are 


o sine’ = 686 sinQ-+051 sin2Lo+20"04T-+A, 
(6) 
e’— €=9"21 cosQ+0"55 cos2Lo+B, 


where 2 and Lo are, respectively, the longitude of the 
ascending node of the lunar orbit referred to the 
ecliptic, and the mean longitude of the sun. In the 
expressions (6) T denotes years from the initial epoch, 
and A and B are additional constant terms chosen so 
that the two quantities o sine’ and e/—e may vanish 
at the initial epoch. 

Then 6 cos(V—WNo) and 6 sin(N—Np»), appearing in 
the right-hand sides of the differential equations (4), 
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have the following form: 


6 cos(N— No) =8"0 cos(N—Q)-+1"2 cos(W+0.) 
+0"5 cos(N—2Lo0)+20"04T sinV 
+C cos(W¥—Q), 
()a 
6 sin(N—N»)=8"0 sin(N—Q)-+1"2 sin(N-+0) 
+0°S sin(W—2Lo)— 20"04T cosN 
+C sin(V¥—Qp), 


where C and Q) are constants related to A and B by 


A=C sin, 
B=C cosQo. 


On the assumption that V is a known function of 
time and that N= —A.n cos//p’, the expression of 6J 
is derived by the usual method of approximation as 
follows: 


N N 
6J =8"0—— cos(N—-Q) +17 2- cos(N+2) 


N—-% N+% 


+075 cos(N—2Lo)+20704T sinN 


N—2no 
% ” 


+C cos(WV—Q))+ is 
N 


cosV, (8) 


where N, %, and no are, respectively, the daily mean 
motions of V, Q:, and Lo expressed in degrees, and 


= —0°053, 
no=0°986. 


For a typical close earth satellite whose inclination is 
not near 90°, all terms except for cos(W—2Lo) have 
nearly the same period, and both N/(N-20¢ and 
N/(N+QQ are near unity. 

From Eggs. (8) and (2), inequalities appearing in the 
inclination to the equator of date are derived as follows: 


i=8" kts cos(N—Q—) — 172 cos(V¥+Q,) 
Nc Q NAM 
nO § 
+1”%0 ———_ coslV. 


cos(N—2Lo)+ 
N-—2no N 


In this expression (9), terms having large amplitudes 
in (8) have disappeared. And, if N=—3°, the largest 
term in (9) is the last one, with an amplitude of 1”. 

If the mean value of 7 is adopted as the inclination 
appearing in the expression of N, ¥, and m as in the 
previous paper, 6/ in the right-hand sides of (4) can 
be written as 6’+6cos(V—No), since the relation 
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between J and iy is 
J=ipt6i+6 cos(V—WNo), 


where 7 is the mean value of 7. 
Then dN, dy, and 6M are easily derived from Eqs. (4) 
-as follows: 


6N = —F coti—G tant, 
6~=F csci--5G sinz, 
6M = 3G(1—e?)? sini, 


(10) 


where 


N 
f=8:0——- sin(N—-0)+172— 
N-M% N+AND¢ 


sin(V¥+00 


N 
+0. >—— 


—2no© 


sin(WV—2Lo)—20"04T cosN 
uv 


eg sinV+C sin(N—Q), (11) 


G=8"0————- sin(N—Q,) — 1"2 _—— 


eID ey ei Ly) 
(N-Q)? (N+0,)? 
Nno Buta 
+1”0——_—__—- sin(V— 2Lo0) +—— sinJV. 
(N—2no)? N 


By combining Eqs. (2) and (10) expressions of 
inequalities in the elements referred to the moving 
equator are derived as 


6Q=—A coti—G tant, 


(12) 
6w= H csci-+-5G sint, 
where 
Q¢ 
H=8'0—— sin(N—{)—1*2- sin(VN+2¢ 
N-% N+Q 
— sinN. 


HO) 
+1"0—— sin(WV— 2Lo0)+— 
N—2no N 


_ The amplitudes in the expressions of 6Q and dw are 
clearly much smaller than those in 6N and 6y for a 
typical close satellite. 


PRP ECGrs 623 


DISCUSSION 


As an example, numerical expressions are given here 
for a case in which N= —3°0 and i=45°. 


d1=071 cos(V—Q,)—0"2 cos(N—2Lo)—1"0 cos, 
6Q= —0°2 sin(V¥—Q()+0"3 sin(N—2Loe)+2"0 sin, 
dw=075 sin(WV—Q,)—0"6 sin(V—2Lo)—4"9 sin, 
6M = (1—e)*(0"2 sin(V—Q.)—0"2 sin(NW—2L) 
—2"1 sinN]. 


The coefficients in these expressions are very small, but 
not negligible for some satellites. Furthermore, if a 
fixed equator is adopted as the reference plane, in- 
equalities appearing in the orbital elements become 
much larger. 

The difference between the perturbations appearing 
in the two coordinate systems is, in essence, the 
following. In the fixed coordinate system, the precession 
and nutation lead to an asymmetric potential field and 
produce secular perturbations which typically amount 
to something more than 20 seconds of arc during the 
course of a year. By a shift to a coordinate system 
attached to the moving equator of date, the major 
perturbations are eliminated and the typical values 
are reduced to several seconds of arc. 

The writer therefore suggests that, for the analysis 
of a close earth satellite, the inclination and the 
argument of perigee can most conveniently be referred 
to the equator of date, and the longitude of the node 
be measured from a fixed point. Explicitly, the node 
should be measured from the equinox of an initial 
epoch along the fixed equator to the ascending node of 
the moving equator, and then along this equator. 
Note the assumption that the initial epoch, to which 
the fixed equator is referred, is not far from the date. 

This system of coordinates has already been adopted 
by G. Veis and C. Moore of the Smithsonian Astro- 
physical Observatory in their orbital improvement 
program for artificial satellites, and the orbital elements 
given in the Smithsonian Special Reports are referred 
to this system. 
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V=—1/r+J2P2(siné)/+J4P4(sind) /r®, 


with J. and J4 assumed to be small quantities of the first and the second orders, respectively, and with the 
value of the orbital inclination 7 lying in a neighborhood of tan2~63°4. The method of attack is based on 
the removal of the short-periodic terms from the Hamiltonian by the-von Zeipel method, followed by a 
Taylor series expansion of the energy integral up to quantities of the second order. As far as the Delaunay 
variables G’, g’ are concerned, the motion then becomes formally identical with that of a simple pendulum, ° 
and the solution is reduced to elliptic functions. In this form all the essential features of the motion are 


clearly revealed. 


1. INTRODUCTION 


HE published solutions (Brouwer 1959, Garfinkel 

1959, Kozai 1959) of the satellite problem 
generally contain long-periodic terms with a critical 
divisor 5 cos’%’—1, where z is the “mean” orbital 
inclination. Such a solution therefore fails to represent 
the motion with sufficient accuracy in some 
neighborhood 


(1) 


of the singularity at i=7*. A nonsingular solution, as 
given by the author in his 1959 lectures at the Summer 
Institute in Dynamical Astronomy, is reproduced here 
with a few simplifications. The method involves the 
removal of the short-periodic terms from the 
Hamiltonian, followed by a Taylor series expansion 
of the energy integral. 
Let the problem be defined by the potential 


V=—1/r+J2P2(sind)/7+J4P4(sind)/r°, 


|i—t|<w, im@=tan12, 


2) 


with the constants Jz and J, satisfying the relations 
O<Jo1, Js=O(J2?), and with the orbital inclination 
7 lying in a sufficiently small neighborhood (1). 


2. THE NEW HAMILTONIAN 


A canonical transformation based on the von Zeipel 
(1916) method is used to remove the short-periodic 
terms from the Hamiltonian. The new Hamiltonian 
may be readily obtained from the author’s earlier 
paper (1959) with two modifications: 

(1) Since the long-periodic terms in that paper had 
been removed along with those of short period, the 
former must now be restored. This is achieved by 
setting S,*=0 on page 357. 

(2) Since the long-periodic terms have been shown 
to be independent of the parameters ¢1, cs, c3 of the 
intermediary orbit that was used there as a first 
approximation, these parameters may now be set equal 
to zero. Useful checks are then provided by a com- 
parison with the new Hamiltonian F’ of Brouwer (1959) 


who had used as a first approximation the Kepler 


ellipse, with c1=¢2=c;=0. 
If c;=0, there follows goi=0, w=\1=A2=1; further- 


more, if S;*=0, the new Hamiltonian F’(L’,G’,g’) can © 


be written 


F’=F’+F* =11°+Ri4 R,+E+4R*+6*, (3) 


where R, and R, are the portions of the disturbing — 


function of orders Jz and J,, respectively, ® is of order 
J,?, and the bar and the star denote the secular and the 
long-period parts of the variable, respectively. It is 
understood that the Delaunay variables L, G, g occur- 
ring as the arguments in the right-hand member of (3) 
are to be replaced by the new variables L’, G’, g’. The 
quantities Ri, Re, &, Ro*, &* are given respectively by 


the equations (28), (91), (68), (95), (80) of the author’s 


earlier paper (1959), with c;=0, 2k=Jo, k’=Ja: 


Ri= (1/4) JonG*(3y’— 1), (4) 

Ry= (—3/128)JanG—" (5— 322) (3—30y?+-35y4), Gya 
B= (3/128) Jn —5+-4e-+ Sa? | 

+y?(10—240— 18x2)+4(35-+36a+5a2) ], (6) 

R,*= (15/64) nG-e sin’s(1—7y*) cos2g, (7) 

&*= (3/64) J2nG—é sini (1—15y*) cos2g, (8) 

with the abbreviations 
a= (1—@)!=G/L, y=cot=H/G, n=l aa) 


The expressions (3)—(8) agree, as they should, with 
Brouwer’s (1959), who removed /’ from the Hamiltonian 
in the first of his two von Zeipel transformations. 
Since J’ is absent from F’, we have L’=const, so that 
our problem is. reduced to one degree of freedom, with 
G’, g’ appearing as conjugate variables. Z 

The quantity F*’ and the derivatives of R, with 
respect to G will be needed in the next section. By the 
addition of (7) and (8), 


F*' = R,*¥+6*= —3,J2nGe 


X sin’ cos2glb—2bo(5y2—1) ], (10) 
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CREMMCAL INCLINATION SATELLITE 


with the constants 6 and bo defined by 


j= 1S J 2, 
(11) 

bo=1—70/4, 
where 0 is the Vinti index (Garfinkel 1959, Vinti 1959). 
From (1) it follows that the quantity A defined by 
A=|5 cos%i—1| satisfies the inequality A<4w; accord- 
ingly A is a small quantity, say of O(J2"), with the 
number k to be specified later. Then sin’s=4+O(J,"), 
and (10) can be written 


F¥'= — (3/40) J:%bneG* cosdg-+O(Js*). (12) 
Next, by the differentiation of (4) with respect to G, 
Rig= — 3 J nG— (Sy’—1)=O(Jo!*), 


: (13) 
Rige= 35 nG-*>+O (J+). 
3. THE ENERGY INTEGRAL 
Since # is absent from F’, 
F’=const (14) 
is the integral of energy. Let 6 be defined by 
6=G’—G", (15) 


where G”’ is a constant to be chosen later, and let F’ 
be expanded about G” in a Taylor series in powers of 6. 
In view of (3), Eq. (14) becomes 


Fl qb + AF enqn 8+ F*' ++ +-+c=0, (16) 


with the constant term F’(G’’) absorbed in the con- 
stant c. The dominant terms in (16) are given by 


We = Rig: --= Bs: - , 
FP’ eng =Rigngu + etene = B/+ vee, 
F*'(G")= A cos2g’, 


(17) 


where A, B, B’ are constants corresponding to the 
dominant terms in (12) and (13): 


A=— (3/40)J2bne’G~, 
B=— (3/4) JonG~(5 cos’i— 1) =O(J""**), 
B'= (3/2)JnG-*>0. 


(18) 


It is understood here that L, G are to be replaced by 
L’, G” wherever they occur in (18) and (9). Now (16) 
can be approximated by 


A cos2g’+ Bé+3B'e+c=0, (19) 


Note that if i=z, then B=0 from (18), and 
5=0|A/B’|*=O(J-}); on the other hand, if i¥7«, then 
5=0(B/B’)=O(J2*). It is therefore natural to choose 
k=%, so that the quantities |5cos%’—1| and 6 are 
both of order (J). Then all the terms of (19) are of 
the second order in J2, and the terms rejected are 
of O(J23). 
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The canonical equation in G’, 
G! =F '= —2A sin2g’+0(J.28), (20) 
becomes, in view of (15), to O(J.2), 
6=—2A sin2g’, (21) 


which, together with (19), determines 6 and g’ as 
functions of the time. 

In the classical solution A is proportional to the 
coefficient of the singular long-periodic term. When 
either e=0 or 6=0 in (18), this term vanishes, and the 
problem of singularity disappears. Indeed, with A=0, 
(21) is satisfied by 6=const, while the canonical 
equation in g’, 


j= Fe = —B+0(2)), (22) 


has the solution g’= — Bt, with —B in (18) furnishing 
the classical expression for the secular rate of the 
argument of the pericenter, to the order J,*. Without 
any loss of generality, it can therefore be assumed here 
that A¥0. 


4. CHANGE OF VARIABLES 


It will appear later that the location of the point of 
stable equilibrium depends on the sign of the Vinti 
index 6 (Garfinkel 1959, Vinti 1959). To make the 
equations of motion independent of 0, define a new 
variable é: 


&=g'—in(1+sgn A). (23) 


If b>0, then 4<O and €=g’; on the other hand, if 

b<0, then A>O and £=g'—}. Therefore cos2g’ 

=—sgn A cos2é, so that (19) becomes 
—|A|cos2é+ Bo+4B'e+c=0. (24) 


Now let G” and the origin of time be so chosen that 


60) =05(0)=0. (25) 
Then c=|A|, and (24) can be written 
2|A|sin2é-+ Bd-+3B'2=0. (26) 


It has been shown that 6 is of order (J2)3. Noting 
that A+0 and that B’>0, define new constants: 


C= |5/J2b|}, 


a=—B/|4AB’|?= (5 cos’*—1)(p/2e)C, (27) 
wo= |44 B’|?=3Joimep|b/5|}, 
where p=a(1—e?)=G’, and a new variable 7: 
n= | B’/4A |*6—a. (28) 


Note that the parameter a is a measure of the difference 
|i—ix|, and that 7 is a linear function of G’; the scale 
has been magnified so that both @ and 7 are of order 
unity. In terms of the new variables the equations of 
motion (26), (21), and the initial conditions (25) 
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become 
sin’é+7?=a? 
n= wo sin2é, (29) 
£(0)=0, 1(0)=—a. 


5. SOLUTION OF THE EQUATIONS OF MOTION 


Equations (29) are identical in form with the 
equations of motion of a simple pendulum. If @ is the 
central angle, wo the angular frequency of zero oscilla- 
tions, and a? is the ratio of the actual energy to that 
required to raise the pendulum to the top, perfect 
analogy is established by 


6=2, 6/wo= —2n, 


The phase-plane orbits are furnished by the contour 
lines of the function 


f(&n) = sin?é+ 1. 


Since f is a periodic function of 2£, the domain may be 
confined to —3r<é<}z. The stationary points of v 


btw? sing=0. (30) 


(31) 


are furnished by the condition f:=f/,=0, which is 
satisfied by the two points M and S: 
M: £=0)) H=0%) Fat =0 
(32) 
S: £=$2; y=0; [al=1, 


That M is a minimum point of f and that S is a saddle 
point follow from the values of the second partial 
derivatives of f. With the substitution 2=0 the phase- 
plane diagram becomes identical with that of a simple 
pendulum, as given in the standard textbooks of 
mechanics (e.g., Goldstein 1950). The contour lines 
passing through S, 


n=-kcosé, (33) 


separate the domain into two regions: (1) The libration 
region |a|<1, where £&(t) is periodic, and the -orbits 
are closed curves enclosing M, the point of stable 
equilibrium, and (2) the circulation region |a|>1, 
where &(é) is monotonic, and the orbits are open 
curves; the advance or the regression of the argument 
of the pericenter occurs as a> 1 or a<—1, respectively. 
The formal solution of (29) is 


sinE=a sn (wot, |a|)=sn (awot,1/|a|), 
(34) 
n= —a Cn (wot, |a|)= —a dn(awi,1/|a|). 
The first form is convenient when |a| <1, the second 
when |a|>1; in either case the indicated modulus x of 


the elliptic functions i is less than unity. Note that, in 
view (27) and (17) 


—B=—Rier, (35) 


so that the product aw» in (34) is the classical rate of 
advance of the argument of the pericenter. 


Ady) = 
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In the limiting case |a|=1, (34) becomes 


siné= tanh (awot), 
(36) 
n= —a sech (awyt). 


This is an asymptotic orbit described in infinite time; — 
as t—> », it follows that £— 7/2, n—>0, so that the 
configuration approaches the saddle point S, which is _ 
a point of unstable equilibrium. It is noteworthy that 
the stable and the unstable points are situated on the © 
line £=0 and on the line £=4z, respectively. These are — 
the line of nodes g’=0, m and the line g/=7/2, 3/2 — 
if 6>0, and vice versa if b<0. The latter inequality — 
holds for the earth with the commonly accepted values 
of Jz and Js. From (27), the condition |a| =1 implies 

|i—i«| ~|5 cos**—1| =5 0/5)! (37) 
The ‘“‘width” wo of the Ree. region is thus only a | 
small fraction of the width w of the “resonance” region 
defined in our earlier paper by 


Geile 
=—|— (38) 
2p|120 
Indeed, 
Wy/w=2|6J2|#~0.16, 
(39) 


max w=1°, max wo=560”, 


the numerical values corresponding to p=1+e, e=1, 
and Jo=1.08X10-, J4= —1.71X 10-6, with 6=—0.47 
(O’Keefe ef al. 1959). Somewhat smaller values of w 
and wo correspond to J4= (—1.3-40.2)10~* (King-Hele 
1959), with the resulting value >= —0.1+0.2. 

The general behavior of the phase orbits can be © 
summarized as follows. As 7 increases from 0 to 7/2, 
a decreases from +0 to —o. The value i=i#~ 6324 
corresponds to a=0, with the orbit degenerating into 
a point M. Let & and é denote the linear and the 
periodic portions of £, and w the common angular 
frequency of £* and 7. (1) In the libration region 
|a| <1, 

k=|a|, &=0, &=sin“(@sn aol), 
w= (1/2K)w, 
0<|é|<sin“]a|, 
0<|n|S|al, 


where KK is the complete elliptic integral of the first 
kind. As |a| increases from 0 to 1, w decreases from wo 
to zero; the secular rate remains zero. The amplitudes 
of £ and 7 increase from zero to 7/2 and 1, respectively. 
(2) In the circulation region |a|>1: 


(40) 


Kk=1/la|, &=(a/2K own, 
&*=sin1(sn oe (1/2K Jaw, 
w= (1/K)awo, (41) 


0<| | Ssinfalt— (n/2K"P) 
— sna (n/2K)* 3}, (@—1)!< [n] <a. 


CAR IE IC 


As |a| increases from 1 to ©, w and é increase from 
zero, asymptotically approaching 2aw) and awo, re- 
spectively. Since aw»=—Rig, these results agree with 
the classical theory as |a| — ©. The amplitudes of 
£* and 7 increase from zero, asymptotically approaching 
(8a?) and (4a), respectively. The latter values are 
derived from (41) with the aid of a Fourier series for 
sn (Whittaker and Watson 1952), and the relations 
(42). 


*=§x sin(r/K)awoe+---, «=1/lal, 
(42) 


Ha— (= 1)!]= (4a) +++, 


(3) The asymptotic orbit, |a|/=1, separating the 
' libration and the circulation regions, is described with 


an infinite period. On this orbit.the angular frequency 
® iS minimized at value zero, the amplitudes of the 
oscillation are maximized, and the secular rate £& 
vanishes, as |a| ranges from ~ to 0. 


6. COMPARISONS AND CHECKS 


In a recent paper Hori (1960) treated the problem 
by removing the long-periodic terms from the 
Hamiltonian F’ by the method of von Zeipel, but with 
the generating function S expanded in powers of (J2)?. 
With this modification of the classical solution the 
singularities are avoided. Although he carries out the 
expansion S=So+S:+S1+ ::- explicitly only as far as 
3, higher-order terms can be obtained if desired. His 
derivative S:,, is equal to our 6, and the results are 
equivalent to order (J2)?. 

Further checks of the theory are furnished by a 
comparison of (42) with the dominant long-period 
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terms in Eqs. (84) and (99) of our earlier paper (1959). 
For the combined effect of the second and the fourth 
spherical harmonics these equations yield: 


de! = J obEG4(5 cos’i—1)-? sin2g’+- 5 oy 
(43) 
5* logG’ = zoJ obP’G—* (5 cos’i—1)— cos2g’/+---, 


with L, G, replaced by L’, G’’ wherever they occur as 
arguments in the right-hand member. An agreement 
between (43) and (42) is established with the aid of 
(27) and (28), and the relation G=/p. Thus the 
“resonance” and the classical solutions asymptotically 
merge into each other. 

The form of the solution adopted here clearly reveals 
all the essential features of the motion; i.e., the 
libration, the asymptotic orbit, and the interchange of 
the stable and the unstable points with change of sign 
of the Vinti index. The latter phenomenon has a 
perfect analogy in a simple pendulum, where such an 
interchange is governed by the sign of the acceleration 
of gravity. 


REFERENCES 


Brouwer, D. 1959, Astron. J. 64, 378. 

Garfinkel, B. 1959, Astron. J. 64, 353. 

Goldstein, H. 1950, Classical Mechanics (Addison-Wesley 
Publishing Company, Inc., Reading, Massachusetts), p. 290. 

Hori, G. 1960, Astron. J. 65, 291. 

King-Hele, D. G. 1959, Xth International Astronautical Congress, 
London. 

Kozai, Y. 1959, Astron. J. 64, 367. 

O’Keefe, J. A., Eckles, Ann, and Squires, R. K. 1959, Science 129, 
565-566. 

Vinti, J. P. 1959, J. Research Natl. Bur. Standards 63B, 105. 

von Zeipel, H. 1916, Ark. Mat. Astr. Fys. 11, No. 1. 

Whittaker, E. T. and Watson, G. N. 1952, A Course in Modern 
Analysis (Cambridge University Press, New York), p. 510. 


THE ASTRONOMICAL JOURNAL 


The System of VV Cephei 


VOLUME 65, 


NUMBER 10 DECEMBER, is 
. 


LAURENCE W. FREDRICK* 
Flower and Cook Observatory, University of Pennsylvania, Philadelphia and Sproul Observatory, Swarthmore College, 
Swarthmore, Pennsylvania 


(Received September 9, 1960) 


The system of VV Cephei is studied by combining photometric, spectroscopic and astrometric information. 


The absolute parallax is shown to be approximately 0005. The inclination of the orbit is very nearly 90° 
and the size of the giant M component is on the order of 600 solar radii. Discrepancies between the spectro- Ps 


scopic data and the astrometric and photometric data are discussed and suggestions are made to explain or 


resolve the discrepancies. 


INTRODUCTION 


HE system of VV Cephei presents a rare oppor- 
tunity in that it may be studied as an astrometric, 
spectroscopic, and eclipsing binary. Gaposchkin (1937) 
and Goedicke (1938) have presented data on the 
spectroscopic orbit and van de Kamp (1951) has given 
results of a provisional study of the astrometric 
material. The present study brings the Sproul astro- 
metric material up to date and combines it with an 


analysis of the photoelectric light curves obtained ° 


during the recent eclipse. 


ASTROMETRY 


The system has been observed at the Sproul Observa- 
tory in every season starting in 1938. A total of 1126 
measurable plates have been taken on 316 nights until 
the end of 1958. The exposure times run between one 
and two minutes as VV Cephei is reduced by 5%76 by 
means of a rotating sector. Data for the reference 
stars used in the astrometric reductions are given in 
the following. The dependence values Dep are given 
to show the balance of the reference frame, Dia ‘is the 
diameter of the image in millimeters on the average 
plate and Am is the magnitude difference from refer- 
ence star 2. The spectra are by.. Vyssotsky (inter- 
observatory communication) who notes that star 1 is 
extremely reddened. 


Diam 
Star Dep (mm) Am My Sp 
1 0.334 0.068 +0.08 10.68 G 
2 0.213 0.076 +e 10.60 KO 
3} 0.210 0.098 —0.38 10.22 B8 
4 0.243 0.081 —0.03 10.57 F5 
VV Cep tee 0.073 +0.05 10.65 M2 


The observations, measurements and reductions 
were carried out in the usual Sproul fashion (van de 
Kamp 1942). Because of the enormous amount of 
material, the nights have been combined to form 
normal points where the night weights have a sum of 
approximately twelve. Thus each normal point on the 
average contains four night’s observations. The normal 
points are given in Table I, where the columns are (1) 
the year of observation, (2) the interval of observation 


* Now at Lowell Observatory. 


for the normal point, (3) the sum of the weights, (4) 
the mean epoch reckoned from 1950.00, (5) and (6) 
the mean parallax factors in a and 6, (7) and (8) the 
mean elliptical rectangular coordinates in X and Y 
(for periastron passage 1951.2), (9) the mean displace- 
ment £ in a, and (10) the mean displacement 7 in 6. 
The values of £ and 7 are given in units of 10 mm 
and a color effect in & (Lippincott 1957) has already 
been allowed for. Figure 1 presents the astrometric 
data plotted as yearly normal points in order to 
minimize the parallax effect. The orbital motion is 
clearly superimposed on the proper motion. q 

In order to reduce the computational time to a 
minimum, ten supernormal points were formed. This 
virtually eliminates any parallactic effects and allows 
the material to be represented by the equations of 
condition where the bars indicate that supernormal 
points are being used and the symbols have their 
usual meaning (van de Kamp 1951). 


E=cotuatt (B)xt+ (Gy; 1= Gti (A)et (F)y. 


Least-squares solutions for the unknowns were made 
using P=20%34, e=0.20, and T=1942.56, these 
elements being adopted from the papers of Gaposchkin 
and Goedicke cited earlier. The reliability of the 
results were then tested by varying the eccentricity. 
The sum of the squares of the residuals from these 
solutions indicate that either the eccentricity is zero 
or that the astrometric material can better be fitted 
by assuming a different value of P or T. Since P is 
probably known with an accuracy far greater than T, 
it was assumed that T should be investigated. 

Using the yearly normal points shown in Fig. 1, 
the proper motion in « and in y was removed graphically 
and T and e were determined by means of displace- 
ment curves (van de Kamp 1947), thus yielding 
T=1951 and e=0.50. The eccentricity was then varied 
through a large range, and least-squares solutions 
were made using the supernormal points. The sum of 
the squares of the residuals suggest that T=1951.2 
and e=0.50 should be used. 

Adopting P=20934, e=0.5, and T=1951.2, a 
solution was then made using the normal points in 
Table I and the equations of condition. 
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Year 


1938 


1939 


1940 


1941 


1942 


1943 


1944 


1945 


1946 


1947 


1948 


1949 


1950 


1951 


1952 


1953 


Taste I. Astrometric normal points. 


ViVi 1G Pie Eat 


Interval w t ips Ps a y é n 
Jul. 7-Aug. 15 16 E1141 23 +.96 eS iiny +.42 ey mate 
‘Aug. 26-Sep. 4 13 —11.34 148 +.94 Ao {38 +41 sais —49 
Oct. 15-Nov. 9 6 24118 85 Wy —1.40 + ,38 abit —24 
Jul. 2-Aug. 6 13 —10.46 “150 + .86 as +.26 ap —20 
Aug. 12-Sep. 3 12 BOSS — .09 + .96 —1.46 4 24 +18 ping 
Sep. 5-Sep. 21 11 —10.30 40 +84 —1.46 4.24 0 =19 
Nov. 16-Dec. 7 8 —10.10 94 =e17 mea +.20 ai 20 
Jun. 22-Jul. 9 11 — 9.50 +.68 +,75 —1.49 + .09 23 =19 
Jul. 14-Jul. 21 12 — 9.46 +.50 + .88 —1.49 +.09 +10 =z 
Aug. 1-Aug. 5 11 = i0t Al 4295 AX07 —1.50 +.08 +15 = 15 
Aug. 9-Aug. 20 11 = 0138 +.07 HL .08 —1.50 B07 0 27 
Sep. 2-Sep. 15 11 20.32 —.29 +.90 ==1-50 + .06 +20 ag 
Sep. 18-Oct. 10 11 — 9.26 — .58 aEe 70 —1.50 +.05 aid my 
Oct. 12-Nov. 28 12 = O16 —.85 +.16 51350 +.03 oe = 9 
Jun. 27-Aug. 3 11 218,46 S248 +.86 —1.49 —.09 stab, + 6 
‘Aug. 8-Aug. 17 14 — 8.39 aah + ,98 —1.49 = 10 +15 6D 
Aug. 29-Sep. 9 15 RR, IG oer 01 —1.49 a ty 16 
Oct. 11-Nov. 17 10 = 819 = 81 HE 33 Peo =) i4 1 eh 
Jul. 16-Aug. 8 ~—12 ee +.40 e202 —1.45 ei ~16 455 
Aug. 27-Oct. 8 15 7220 AO + .80 —1.43 30 ei ANG 
Jul. 1-Sep. 6 12 6.42 Ey +.88 —1.36 eye ry) 0 
Sep. 13-Nov. 18 11 — 6.20 =i +39 ean 3d = 48 eS +12 
Jul. 24-Sep. 2 15 — 5.39 +.10 +.96 =1,99 — .60 ai9 sae 
Oct. 22-Dec. 4 9 ists = 68 =.02 = 1.18 — .63 = td 417 
Jun. 23-Jul. 12 11 — 4.49 +.65 B78 —1.06 eit a) “etd 
‘Aug. 12-Sep. 28 11 54134 2,10 88 =1..03 SE 6 ae 
Oct. 16-Nov. 18 17 24,10 £783 +.34 = 1:00 ==, 75 —16 eel 
Jun. 27-Jul. 18 8 273) 47 +.56 + .83 ae ey: ao = 8 
Jul. 31-Aug. 12 12 — 3.40 224 +.98 E82 282 = 15 +17 
Aug. 24-Aug. 31 12 pte .34 ants +.96 = 5180 2183 = 10 as 
Oct. 3-Nov. 27 9 =e, 18 278 a0 eet = 84 a4 aaa 
Jul. 13-Aug. 11 12 =), 43 teen + .94 = {55 — .86 ee) =a 
Oct. 20-Oct. 26 9 at) 164 + .36 aT 2H =12 me 
Jul. 8-Sep. 16 (ee 130, 12 +.88 Sat — 83 oe 0 
Oct. 21-Nov. 14 11 ze a 2.88 ag 14 ey ~13 eal 
Aug. 9-Aug. 10 6 = 0,40 +.16 +.98 etd. =) —30 416 
Sep. 20-Sep. 25 12 0.27 =),53 +.76 in AB os =~ 13 
Oct. 1-Oct. 13 12 = 0723 = 70 +.57 op 2,62 —10 = 
Oct. 20-Nov. 20 12 aoa 200 +.09 an) 22 — .60 6 arf, 
Jul. 18-Aug. 9 12 PO: 57 +.34 +.94 fey?) iil +14 6 
Sep. 17-Oct. 6 12 No a eT + .46 =) 24 and +9 
Oct. 13-Oct. 30 15 + 0.80 ERT +.41 es E47 20 16 on 
Nov. 1-Nov. 18 14 + 0.86 03 +.07 ed >is a6 + 6 
Jul. 31-Sep. 15 11 + 1.66 17 +.90 + .46 24 ee x 9 
Sep. 23-Oct. 5 11 + 1.74 2560 +.69 ends + .28 10%) — 6 
Oct. 13-Nov. 20 13 eaves aed + .29 + .43 EP, = 4 0 
Jun. 25-Jul. 13 2) "429-50 +.68 +.74 4 iuiGR +18 ag 
Jul. 15-Aug. 19 6, + 2.60 +.20 +.95 ot +.61 +11 +10 
‘Aug. 23-Sep. 1 13 + 2.56 1,12 +.96 + .19 +.63 417 +56 
Sep. 5-Oct. 5 11 oh 39 = .45 +.80 a. 17 +.65 +14 1 
Oct. 6-Oct. 23 13 + 2.80 eee +.44 ee 14 +.67 a9 +15 
Oct. 24-Oct. 30 11 9 83 ay) ae 98 + .13 + .68 +7 ae 
Nov. 12-Nov. 24 9 288 eek — .06 tlt +.69 Eas =f 
Jul. 5-Oct. 3 12 + 3.68 hr Te Aaa TS Beano +10 —10 
Oct. 7-Oct. 12 12 Boe ey =) +.56 aN; +.83 +9 —2 
Oct. 16-Nov. 1 12 SBE eG. aed 1,32 ey. 583 ARS, =e 
Nov. 2-Nov. 19 12 +) 3.87 ay cot = 94 + .84 se) ag 
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Taste I. (Continued) 


Year Interval w t Ps Ps © y iB n 
1954 Jul. 27-Aug. 13 12 + 4.59 + .26 + .96 2S Sy, eae 0 
Aug. 14-Aug. 31 14 + 4.64 — .06 LO =" 49 S87 avd, + 6 
Sep. 2-Sep. 29 12 oe aD ei oe a7 = Ae Pee ey, ey | as 
Oct. 23-Nov. 15 15 + 4.84 — .90 +.15 — .55 +.86 0 +4 
1955 Jul. 13-Aug. 6 9 + 5.56 + .40 TEN) = 76 +.84 3 +10 
‘Aug. 27-Sep. 21 11 + 5.69 eee + .86 ey, ERS aed atti) 
Sep. 28-Oct. 17 14 + 5.77 —.69 +.59 — .8i +.82 —2 + 8 
Oct. 20-Nov. 18 12 + 5.84 — .88 +.21 — ,83-- + .82 0 +11 
1956 Jul. 28-Aug. 22 10 + 6.57 + .32 +.88 —1.00 +.75 — 1 +17 
Sep. 21-Nov. 5 12 + 6.80 eT 42,39 22405 oe) = +13 
Noy. 6-Nov. 27 11 + 6.87 — .92 0 —1.06 +.72 —4 + 5 
1957 Sep. 11-Oct. 11 10 + 7.76 —.64 +.63 —1.22 + .60 —16 +21 
Oct. 20-Dec. 2 8 + 7.84 — .88 — .02 —1.23 +.59 —13 +20 
1958 Jun. 26-Jul. 1 12 + 8.49 EY +.70 a Ue +.49 +4 +24 
Aug. 9-Aug. 21 14 + 8.62 +.05 + .98 —1.34 +.47 —22 +20 
Noy. 7-Nov. 12 12 +. 8,86 = 192 + .06 — 1.86 + .43 =%6 +21 
mm I 
+0.002- 
+0.001 - O 
Oo ri) 
As 
fe) fe) O 
-0.001- O St ae 
fe) 
pee Sum of weights Scale 
< 30 ° 
” mean probable error 
30-45" « O P 


1940 


> 45 Ge) 


1945 1950 1955 
Fic. 1. Mean yearly displacements £ (R.A.) and y (Dec.). 


1960 


WV Ci Part Bal 631 


~ 0.003 0:004 0:005 a 


Fic. 2. The parallax mw as a function of mass ratio for the extreme values of 6 and for a:=13.3 a.u. 
Dashed line on the left is the limiting value for 7. 


£=c,+yut+7,P.+ (B)x+ (G)y, (G), (A), and (F) yield the elements 
N= Cyt Myt+ayPst (A)x+ (F)y¥. a=0.00178mm =0"0336-+0"0030 
i= 90°39’ 
The coefficients determined by the method of least ait 


squares are 
The values of the relative parallax in « and y are in 


G2= +0.00027 ve good agreement. Assuming the present parallax deter- 
cy= —0.00017 : ‘ minations from 1951 on have no systematic errors, a 
oe Rion = Foyer ean pia study by Miss Lippincott (1957) would correct the 
t= +0.00046=+0.0086 -+0.0026 value of the relative parallax in « to 070041, which is 
Gino 00071= 0 0134 aD ose in excellent agreement with the value in y. Adopting 
(G)=+0.00116=+0.0219 +0.0063 this correction the combined relative parallax is 
(A) =—0.00063=—0.0119 +0.0037 
(F)=—0.00098=—0.0185 +0.0063. r= +0"0042+0"0020. 


The conversion to seconds of arc is made by using the The reduction to absolute parallax may be made by 
Sproul scale factor, 18”87/mm. The coefficients (B), using the precepts of the Vyssotskys (1948) yielding 


0.5 


0:004 0005 o!006 TT 


Fic. 3. Same as Fig. 2 except here a;=10 a.u. 
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Fic. 4. Twenty-day means of visual observations by McLaughlin. The long-term M’ wave is indicated by the dashed line. 


FRE DR VC K : 
a correction of +0"0010 or by using the precepts of 
Binnendijk (1943) yielding a correction of --0"0015. 
The value +070012 is adopted here, giving for the 
combined absolute parallax, 


= +070054+070023. 


This is'a direct trigonometric parallax. Another ang 
more reliable value for the parallax may be deriv 
using the amplitudes-of the photocentric orbit given 
in the foregoing. 

The fractional distance 6 of the M star to the photo- 
center is given by (van de Kamp 1951) 


B=ls/ (m+n), 


where the subscripts M and B refer to the M- and 
B-type stars involved. Normalizing Im+/p to unit 
light gives @=/p. This can then be obtained from the 
depth of eclipse curve given in Fig. 12 and discussed 
later. At the effective wavelength of the Sproul re- 
fractor we find 0.13<@<0.21, the uncertainty arising 
from the intrinsic light variations, which will be 
discussed later. 

For a given mass ratio Sp/Im, there is the mass 
function B,=9Np/(Itm-+Ms). The spectroscopic orbit 
gives a. Since a=4,/Bz, the substitution can be made 
in a= (B,—£)a which gives 


ee (1—6/Bs)ai, 


where a is the semi-major axis of the photocentric 
orbit with respect to the center of mass. The parallax, 
which is here called a dynamical parallax mg, then 
follows directly from the relation 


na=a’/a. 
Adopting Goedicke’s value for a,sini, since the 


inclination is very nearly 90°, there results a:=13.3 a.u. 
If the relation for zg is written in the form 


ma=al"'/(1—B/Be)an, 


it is possible to investigate ag as a function of the mass 
ratio. Figure 2 exhibits this for the extreme values of 6. 
It is easily seen that for this value of a; the lower limit 
for the parallax is 070026. 


SPECTROSCOPY 


A new analysis for the spectroscopic elements will 
soon be available from the Michigan material. From a 
small-scale reproduction of the radial velocities the 
author has derived the following elements using the 
method of Lehmann-Filhes as presented by Smart 
(1951): 

w= 200° 
e=0.2 
Vo=—18.2 km/sec 
T=JD 2430500 
a, sint= 15.2108 km. 
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TaBLe II. Variable and comparison stars. 


Star BD No. R.A. (1950) Dec. My Spectra 
VV Cephei Variable +62°2007 21555™14s +63°23/13” 5.25 M2 
20 Cephei Comparison +62°2029 22503™29s +62°32'28” 5.14 K5 
tee Check +62°1994 21551™10s +62°28/35”” 6.25 B8 


[t follows that a;=10 a.u. A plot of za as a function of 
mass ratio is given in Fig. 3 for a;=10 a.u. Here the 
ower limit for the parallax is 070039. 

For a range in the mass ratio of 0.7 to 1.0 it is 
evident that the parallax lies between 07004 and 0”006. 
The new spectroscopic material will probably favor 
the smaller value of a, and hence the larger value for 
the parallax. In what follows r=0%005 is used. 


PHOTOMETRY 


McLaughlin made visual estimates of the magnitude 
of VV Cephei on 1539 nights from June 1932 to May 
1957. These observations have been averaged over 
twenty-day intervals and are plotted in Fig. 4. 
Examination of this figure shows very pronounced 
peaks. There are 25 such peaks within the 8380-day 
interval giving an average time between peaks of 349 
days. The peaks have an average amplitude of 073. 
There is also a much longer term variation (here called 
the M’ wave) of the order of 13.7 years with an ampli- 
tude of something like 0715. It is to be noted that the 


% Transmission 


seouve® 
eoesee 


recent eclipse occurred near the top, but on the down- 
ward slope, of this long-term variation. 

The system of VV Cephei has been observed photo- 
electrically at the Flower and Cook Observatory. The 
Pierce photometer (Biltzstein 1958) coupled with the 
15-inch siderostat was used for most of the observations. 
The photometer was not used as a dual photometer 
for these observations. This photometry is defined by 
the following filters: 


Filter Type Wavelength Halfwidth 
Yellow Interference 5295 A 55A 
Green Interference 4890 50 
Blue Interference 4290 65 
Ultraviolet Glass 3800 290 


The interference filters were fabricated by Baird- 
Atomic, Inc. and purchased through a grant from the 
Swarthmore College Faculty Research Fund. The 
transmission intensities were traced on a Beckman 
spectrometer and are shown in Fig. 5. 

Observations were made on 131 nights. The observa- 
tions were made in the order: sky, comparison star, 
variable, comparison star, variable, etc., in all wave- 


Objective —— 
Ultraviolet filter —.— 
Objective plus UV filter -—- 
Blue filter 

Green filter ---- 

Yellow filter ——-- 


Fic. 5. Transmission characteristics of the filters and objective defining the author’s photometry. 
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TaBLe III. Photometric normal points. 


Yellow Green Blue Ultraviolet 
JD Epoch Am n Epoch Am n Epoch Am n Epoch Am 
2435438 .5328 —0.117 3 BSRBS +0.013 3 5484 —0.164 2 
5441 .4838 —0.114°3 .4854 +0.025 3 -4868 —0.144 3 .5359 —1.222 
5442 .5249 —0.129 2 patil +0.003 2 . 5289 —0.158 2 -4868 —1.369 
5448. 5134 —0.175 2 .5148 —0.065 2 25159 —0.255 2 .4680 —1.414 
5452 .4847 —0.149 2 4859 —0.016 2 -4873 —0.255 2 
2435679 .8075 +0.024 3 -8562 +0.446 1 . 8289 —0.317 
5681.7797 —0.019 6 . 7574 +0.162 4 . 7788 +0.447 3 7384 —0.237 
5682. . 7505 —0.254 
5683 .7799 0.000 5 8018 +0.162 3 .8268 +0.421 3 7874 —0.233 
5684. 7156 —0.022 3 .6973 —0.125 
5685. 7652 +0.002 10 7553 +0.216 4 .7778 +0.488 3 . 7936 —0.171 
5700 .6702 +0.005 4 
5701.7498 +0.018 8 7596 +0.004 
5703 .7695 +0.030 5 
5708.7133 +0.023 7 . 7494 +0.155 3 7134 +0.521 3 6947 —0.035 
2435709 .7143 +0.022 8 . 7108 +0.043 
5710.7274 +0.023 5 .7011 +0.052 
SyAlehs 7454 +0.066 
5713 .6864 +0.006 1 - 6480 +0.157 
5715.6799 +0.010 1 
5721. .6610 +0.108 
5725 .6626 +0.024 5 6186 +0.168 2 .6846 +0.533 3 6323 +0.067 
5726.6515 +0.033 6 .6782 +0.185 3 .6910 +0.545 2 .7134 +0.107 
5729 .6638 +0.018 14 
5735.6367 +0.062 6 .6760 +0.205 1 6380 +0.584 6 .6605 +0.142 
5737 .6453 +0.072 4 6454 +0.598 4 .6457 +0.140 
5738.6251 +0.054 1 6239 +0.568 1 .6227 +0.081 
5754. 5871 +0.168 1 .5883 -+0.703 1 .5913 +0. 203 
5755.5784 +0.178 1 6036 +0.702 1 .6038 +0.232 
5756.5819 +0.183 8 - 5830 +0.704 9 .6079 +0.301 
5758 .6203 +0.198 1 .6212 +0.675 1 6221 +0.327 
5762.5500 +0.185 1 5749 +0.711 1 
6019 - 7807 +0.575 
6020 . 1667 +0.608 
6021 7860 +0.613 
6022 . 7942 +0.594 
6025 ‘ 7871 +0.576 
6027 . 1336 +0.525 
6030. . 7618 +0.488 
6034 7415 +0.475 
6035 . 7858 +0. 433 
6037 7181 +0.470 
- 2436038. . 7408 +0.400 
6039. . 7470 +0.432 
6040. 7574 +0.451 
6051. BChey il +0.682 2 7232 +0.378 
6052.7218 —0.019 5 7245 +0.673 4 1675 +0.524 
6053. 7557 —0.026 2 ' . 7676 +0.659 2 . 7833 +0. 490 
6056. 7093 —0.031 5 Sth +0.658 6 . 7538 +0.501 
6061.7185 —0.034 5 . 7246 +0.172 4 7203 +0.670 5 .6716 +0.511 
6062. 7312 —0.054 2 aleVAl +0.162 2 7330 +0.647 2 .7021 +0.498 
6063. 7135 —0.055 5 7144 +0.178 5 7135 +0.666 4 -6635 +0.527 
6074.6792 —0.053 5 .6818 +0.191 4 6783 +0.663 5 7312 +0.541 
6078. 6368 —0.047 7 .6379 +0.167 7 6428 +0.654 7 6641 +0.545 
6087 .6279 +0.014 7 6288 +0.224 7 -6296 +0.691 7 6631 +0.571 
6101.5842 +0.106 5 5850 +0.320 5 . 5860 +0.765 5 - 6082 +0.585 
6106. 5894 +0.145 2 5908 +0.372 2 -5923 +0.781 2 -6255 +0.673 
6107 .5438 +0.143 3 .5447 +0.355 3 -5456 +0.781 3 - 5666 +0. 603 
6108. 6038 +0.139 2 6047 +0.368 2 .6056 +0.793 2 -6285 -++0.650 
6109.5715 +0.166 6 5724 +0.374 6 .5733 +0.796 6 -6085 +0.648 
6115.5919 +0.204 2 .5944 +0.403 3 5954 0.825 3 6230 +0.648 
6121.5402 SRA 5) .5412 +0.436 5 5421 +0.834 5 5937 +0. 639 
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TaBLeE III. (Continued) 


Yellow Green Blue Ultraviolet 
JD Epoch Am nN Epoch Am n Epoch Am Nn Epoch Am n 
6122.5956 +0.229 3 .5966 +0.438 3 5975 +0.842 3 .6278 0.613 3 
6123.5360 SO.220° 3 . 5369 +0.439 3 5379 +0.833 3 .5657 t0.635 3 
6124.5507 +-0.237 5 .5516 +0.454 5 5525 +0.848 5 5835 +0.643 5 
6125. 5484 +0.229 3 . 5485 +0.446 3 5449 +0.844 3 OS) +0.649 1 
6126.5339 +0.242 6 .5348 +0.441 6 5357 +0.846 6 .5729 +0.589 6 
6131.5188 +0.245 3 .5159 +0.437 5 5168 +0.798 4 .5631 +0.573 5 
6132.5502 +0.219° 5 .5415 +0.440 5 5423 +0.786 6 .5923 +0.577 4 
6134.5171 +0.228 2 .5180 +0.420 2 5189 +0.771 2 .6201 +0.490 2 
6140. 5583 +0.190 8 .5670 +0.405 6 5520 +0.665 6 SVE) +0.313 5 
6149 4981 +0.195 4 .5294 +0.412 5 5218 +0.781 4 .5144 +0.586 3 
6152.5718 +0.192 4 Oma +0.392 4 5667 +0.772 2 .5728 +0.474 3 
6153.5442 +0.192 6 .5451 +0.392 6 5460 +0.747 6 .4948 +0.483 5 
6154.5506 +0.218 5 .5506 +0.407 6 5451 +0.776 5 4871 +0.500 7 
6155.5389 +0.189 7 . 5398 +0.385 7 5407 +0.773 7 4848 +0.478 5 
6158.5014 +0.175 1 .5023 +0.385 1 5032 +0.700 1 -4781 +0.370 4 
6160. -4925 +0.329 1 
6163.5251 +0.166 3 .5259 +0.385 3 5267 +0.627 3 4783 +0.234 6 
6164. 4635 +0.156 3 
} 6167 .5063 +0.290 4 .5072 +0.399 4 5082 +0.576 4 -4695 +0.054 4 
6168 .5530 +0.170 1 .5599 +0.392 2 5608 +0.631 2 5434 +0.167 1 
6169.5507 +0.174 4 .5511 +0.397 4 5471 +0.688 3 5045 +0.171 5 
6172.5207 +0.188 7 .5208 +0.405 7 5217 +0.707 7 -4675 +0.271 4 
2436174.5211 +0.194 8 .5186 +0.401 8 $231 +0.684 7 4685 +0.248 4 
6175.5170 +0.171 6 .5209 +0.400 7 5185 +0.677 6 4715 +0.243 5 
6178.5353 +0.195 3 BO2a2 +0.410 3 5371 +0.687 3 4564 (+0.314) (1) 
6185.5149 +0.189 2 5111 +0.363 1 5166 +0.610 2 .4926 +0.046 3 
61884954 +0.154 2 .4963 +0.360 2 4972  (+0.587) (2) .4675 +0.021 4 
6190. .5060 —0.035 4 
6194. .4604 —0.122 2 
6195. 4875 +0.141 2 4884 +0.314 2 4893 +0.514 2 4683 —0.140 2 
6200.4791 +0.127 1 -4800 +0.296 1 4809 +0.537 1 .4628 —0.104 3 
6202.4709 +-0.127 1 .4718 +0.293 1 4727 +0.481 1 .4700 —0.224 1 
6205. 4627 —0.323 2 
6207.4719 +0.150 1 .4731 +0.261 1 4771 +0.384 2 -4657 —0.390 2 
6208. 4734 +0.075 1 .4743 +0.262 1 4755 +0.346 1 4674 —0.428 2 
6209. 4704 +0.071 1 .4713 +0.192 1 4722 +0.305 1 4692 —0.569 1 
6210. 4696 +0.046 1 4705 +0.202 1 4714 +0.252 1 -4687 —0.613 1 
6212.4677 +0.052 1 .4668 —0.662 1 
6370. 7980 —0.123 3 . 7988 —0.016 3 7996 +0.043 3 7972 —1.251 3 
6371.7936 —0.128 2 . 8005 —0.017 3 8013 +0.053 3 .8032 (—1.189) 1 
6372.7831 —0.120 2 . 7743 —0.017 1 7906 +0.060 3 7823 —1.234 2 
6382.7583 —0.162 5 7592 —0.037 5 7600 —0.030 5 7954 —1.368 6 
6387 .7558 —0.129 2 .7425 —0.034 3 7473 +0.032 4 
6394. 7368 —0.112 6 . 7380 +0.024 7 7388 +0.123 7 7938 —1.027 7 
6401.7345 —0.127 3 .7276 —0.002 4 7361 +0.112 3 - 7693 —1.142 4 
6414. 6491 —0.052 1 
6420.7227 —0.129 5 . 7099 —0.053 3 7203 +0.036 4 . 7690 —1.207 5 
6424.6978 —0.119 4 . 6986 —0.012 4 7004 +0.046 4 . 7304 —1.193 3 
6425.6925 —0.133 5 6885 —0.015 4 6893 —0.011 4 7254 —1.320 3 
6430. 6766 —0.078 6 6854 +0.011 3 6678 +0.064 3 .7144 —1.163 3 
6433 ..6856 —0.040 3 6824 +0.045 7 6778 +0.079 3 .7168 —1.145 4 
6434. 6520 —0.042 2 6437 +0.042 3 6520 +0.071 3 .6701 —1.156 1 
6435. 6884 —0.041 6 .6902 +0.056 7 6913 +0.083 6 . 7300 —1.140 2 
6438. 7441 —0.049 5 6904 +0.073 4 7271 +0.104 4 .7106 (—1.133) 1 
6445.6231 —0.002 5 .6187 +0.111 6 6215 +0.120 7 6637 —1.150 5 
6455 .6542 +0.032 2 6527 +0.177 3 6527 +0.181 3 .6777 —1.156 2 
6457.6148 +0.076 4 .6170 +0.179 7 6135 +0.171 5 .6570 —1.178 7 
6459 .6556 +0.059 4 6504 +0.177 4 6514 +0.158 5 6870 —1.189 1 
6469. 5900 +0.074 4 -5682 +0.174 4 5680 +0.079 5 .6171 —1.354 3 
6470.5712 +0.095 3 .5835 +0.169 3 5837 +0.097 3 .6028 —1.328 3 
6475. 5589 +0.110 5 .5813 +0.184 4 5813 +0.120 4 oe tee e 
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Taste III. (Continued) 


Wier Es MOMRET  CAke 


Yellow Green Blue ‘Ultraviolet 
JD Epoch Am n Epoch Am n Epoch Am n Epoch Am a 

6482.6220 +0.110 5 .6005 +0.186 4 6028 +0.139 3 .5819 —1.201 4 % 
6489 . 6283 +0.045* 7 6036 +0.169 3 6038 +0.146 3 5849 —1.205 3 

*Flare } 

2436493. 5054 +0.031 3 .5568 —1.183 2 : 

6497 .4855 —0.037 3 5004 +0.103 3 .5007 +0.086 3 ES2a2 —1.311 6 1 
6507 .5331 —0.045 1 .5323 +0.071 “4 .5339 —1.201 1 " 
6508 . 5083 —0.047 2 .5467 +0.058 1 -5091 +0.075 2 6225 —1.218 5 
6514. 5364 —0.095 2 5348 +0.028 2 .5356 +0.038 2 .5492 —1.310 2 


lengths. The sky was monitored as often as necessary. 
The integration time was 60X10-> days. On almost 
every night the check star was observed in order to 
check the constancy of the comparison star. Data for 
the stars are given in Table II. The observed normal 
points in the sense VV Cephei minus comparison star, 
corrected for relative extinction, are given in Table ITI. 
The epoch is heliocentric from Greenwich mean noon 
and m indicates the number of observations making 
up each normal point. These are plotted as dots in 
Figs. 6, 7, and 8. The check star minus comparison 
star is shown for each wavelength. 

In order to extend the observing season during 
egress, observations were made by Binnendijk on the 
28-inch . reflector on eight nights at the same time 
observations were being made on the siderostat. Five 
additional nights were then obtained. These observa- 
tions, reduced to the system of the Pierce photometer, 
are given in Table IV, and are plotted in Figs. 6, 7, 
and 8 as crosses. 

In addition to the foregoing observations, G. 
Larsson-Leander (1957, 1959) has observed the.system 
extensively in blue and yellow light. His observations 
are not restricted by the seasons due to his location 
(60°). 

During eclipse the light of the M star varies rather 
smoothly. The long slope at JD 2435760 is similiar to 
one at JD 2436100 and occurs again at JD 2436450. 
The average time between these’ slopes is 345 days. 
This is so close to the 349-day value between visual 


4 
peaks on McLaughlin’s light curve that they no 
doubt represent the same phenomenon. The 345-day 
cycle is made up of two waves, one of small amplitude 
(S wave) followed by one of large amplitude (LZ wave). 
The S wave plus L wave may be represented as the 
resultant sum of a sine and cosine function, but there 
is not sufficient coverage in time to confirm this. 

A curious and abrupt pulselike phenomenon occurred 
at JD’s 2435729, 2436141, 2436166, and 2436469. 
These are interesting because they have larger ampli- 
tudes in the shorter wavelengths whereas the waves 
already referred to have larger amplitudes in the 
longer wavelengths. Table V lists the amplitudes in 
magnitudes of the various variations. In the table the 
depth of eclipse d in the various wavelengths is included 
and plotted in Fig. 12. The wavelengths have been 
rounded off, the amplitudes of the wave were inter- 
polated from Larsson-Leander’s material, and the epoch 
is the central date of the phenomenon. The L wave 
appears to increase linearly with wavelength. The 
amplitudes of the pulses increase with decreasing wave- 
length just as the depth of eclipse does. Because of this, 
it is tempting to associate these with the B star. How- 
ever, A. Deutsch could not find any trace of the B star 
on a spectrogram taken on JD 2436142 and Wright 
(1958) places the system only a week past third contact 
on JD 2436187. These pulses may be due to either a 
shell about the B star or perhaps a fluorescence in the 
atmosphere of the M star due to a brief increase in 
radiation from the B star. The depth of eclipse was 


TasLe IV. Photometric normal points (28-inch reflector). 


Yellow 

JD Epoch Am n Epoch 
2436185. 5494 +0.162 2 .5502 
6188.5197 +0.153 1 .4780 
6194, 5247 +0.138 3 -4922 
6197 .5326 +0137" 2 .5054 
6200.5158 +0.143 2 -4863 
6202. 5087 +0.131 4 4837 
6207 .5820 +0.087 2 6077 
6208. 5274 +0.085 1 .5150 
6209. 4953 +0.066 2 4833 
62135420 +0.075. 1 .5670 
6215.5145 -+-0.055 , 2 .4914 
6223. .5107 
6227.4932 +0.038 1 .4947 


Blue Ultraviolet 

Am n Epoch Am nN 
+0.625 2 .5495 +0.062 2 
+0.604 1 4701 +0.003 3 
+0.520 3. | -4753 —0.143 3 
+0.525 2 .4957 —0.116 2 
+0.512 2 .4762 —0.133 2 
+0.475 2 .4978 —0.238 2 
+0.373 2 
+0.387 1 .5079 —0.433 1 
+0.303 2 .4700 —0.571 3 
+0.323 1 .5767 —0.630 1 
+0.313 2 .4737 —0.584 3 
+0.240 4 -4955 —0.723 4 
+0.276 1 .5100 —0.723 1 
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TaBLe V. Amplitudes. 


Color SY 16 val P» P3 d 
Yellow 5300 (0718) 0=30 0704 0705 0700 020 
Green 4900 vee 0.25 0.04 0.05 0.01 0.28 
Blue 4300 (0.12) 0.20 0.14 0.15 0.11 0.64 
Ultraviolet 3800 tee Ostoa 0223.7 028) 0520' 1.76 


determined by a rectification procedure which will be 
sxplained. 

On JD 2436489 a rather sudden change in the yellow 
light occurred during a 36-minute interval. The light 
became brighter by 0703+ and subsided during this 
time. This is shown in Fig. 9 where the ordinate is in 
instrumental intensity units. A rapid change such as 
this must be quite local in origin,.but it is impossible 
to say which star is involved. Observations in the 
other three wavelengths are lacking because telescope 
time on that night was limited and the observations 
in the other wavelengths were already complete. The 
flare is real as can be seen from the constancy of the 
sky and comparison star. Abrupt changes of 0™1 or 
more have been reported to occur in a short period to 
time, but this is the only evidence for such changes 
found during this investigation. 


SOLUTION 


To find a solution, Larsson-Leander’s light curves 
were used because they are continuous. A straight- 
forward solution is impossible, due to the intrinsic 
variations. Following a suggestion by J. E. Merrill, 
the light curve was rectified by subtracting the S and 
L waves, making allowance for the descending slope 
of the long-term wave M’. The rectified light curves 
were then converted into intensities and folded back 
upon the egress branch and a y solution was made 
following the method of Russell and Merrill (1951). 
The resulting solution in blue is shown in Fig. 10. 
The elements obtained are 


Blue Yellow 
Ee 0.6 Assumed 1.0 
k 0.1 0.05 
06 14°81 14°35 
0; 11°24 12°87 
t 92° (88) 91°(89) 
ia 0.23 0.24 
Te 0.03 0.01 
1st contact 2435639 2435678 
2nd contact 5717 5708 
Mid-eclipse 5947 5947 
3rd contact 6177 6240 
4th contact 6251 6270 


These solutions would have reasonable validity provided 
the assumptions about the constancy of the M star’s 
intrinsic variations hold and provided the B star is 
not variable. It is obvious from a comparison between 
the times of the contacts in the blue and yellow given 
in the foregoing that something is wrong. Treating the 
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branches separately it is possible to fit ingress very 
well with a small value, k=0.1, and egress with a 
somewhat larger value, k=0.3. This would suggest 
that the B star changed its effective size drastically 
during the eclipse, and from the shape of the egress 
branch that it changed in brightness also. However, 
this is rather startling and we may look for another 
explanation. 

An alternate and consistent solution may be derived 
from Larsson-Leander’s color curve shown in Fig. 11. 
Assuming that the B star’s color does not change 
drastically as its light varies then any abrupt change 
in the color of the system will be due to eclipse. The 
color changes sharply more positive between JD 
2435655 and 2435680, and sharply less positive between 
JD 2436180 and 2436210. The slopes are identical but 
of opposite sign. Rounding of the color curve before 
and after these dates can be caused by a shell about 
the B star. The asymmetry in the bottom of the 
eclipse coincides with the M star’s L wave and may 
indicate an abnormal brightening in the longer 
wavelengths. 

It is easily seen on inspection of the light curves 
that third contact in the ultraviolet light occurred 
long before third contact in the other colors. This 
adds support to the notion that the shell causes the 
rounding off of the color curve. If this is so we can 
estimate the times of contact for the B star and its shell. 
We find 


B Star Shell 
1st contact 2435656 2435626 
2nd contact 5683 5713 
Mid-eclipse 5931 5931 
3rd contact 6179 6149 
4th contact 6206 6236 


Assuming i= 90°, we have for the B star 
(D—d)/(D+d)=r,/r,= k=0.05, 


and for the shell (subscript, y=ring, or in this case 
= shell), 
1,/%g=k=0.17. 


a 
JD 2436489.6300 


ry 
eo 
° 
° 
| 


Fic. 9. Flare of VV Cephei in instrumental intensity units. 
Open circles are comparison star readings, half-filled circles are 
VV Cephei readings, and filled circles are sky plus dark current 
readings. 
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t then follows that the ratio of the size of the shell to 
hat of the B star is three. McKellar, Wright, and 
‘rancis (1957) have estimated the size of the shell at 
agress to be twice the size of the B star. 

The time of mid-eclipse given by the color curve 
grees within a day with the time of mid-eclipse given 
yy Gaposchkin’s (1937) elements. 


DISCUSSION 


The apparent visual magnitude of the M_ star 
obtained from the light at mid-eclipse) is approxi- 
nately 5725. For a parallax of +07005 and an 
ibsorption of +010 (Risley 1943) the absolute visual 
nagnitude of the M star is —1.36. If the parallax were 
)"004 the absolute magnitude would be —1.84. In 
iny case, this disagrees with the absolute visual 
magnitude associated with the classification M2+Ia 
—Tab assigned to the M star by Keenan and Wright 
(1957). 

The absolute visual magnitude of the B star follows 
rom the depth of eclipse and the relation 


Lp/Lm=0.17/0.83= 2.512136) . 


nence Mg=+0.36. For a parallax of 07004, we get 
Mp=—0.12. Thus the B star in any case is an under- 
luminous shell star. Assuming a bolometric correction 
of 2.2 magnitudes, the radius of the M star can be 
estimated by means of Stefan’s law. The M star’s 
bolometric magnitude is about —4, which leads to a 
radius of 200 Ro, or 300 Ro when Goedicke’s value for 
the temperature is used. Based on the circular model, 
this would indicate that the radius of the orbit is 
about 6 a.u., but it has already been shown that a; is 
of the order of 10 a.u. For the circular model at external 
contact and i= 90°, 


r,=sind./(1+R), 


yr, being in terms of unit radius of the orbit a. However, 
if the orbit is eccentric the relation for r, becomes 


sind. a(1—e?) 
VER 1—ecosy 


tg 


where a is the semi-major axis of the relative orbit and 
vy is the true anomaly at the time of eclipse. 

In the case of the astrometric orbit it can be shown 
that 7,~0.29a. If the mass ratio is 0.8, then 7,~1200 
Ro. In the case of the spectroscopic orbit we find 
r,~0.12a or rz~580 Ro. This latter radius corresponds 
to an effective temperature of approximately 1800° 
according to Stefan’s law when a bolometric magni- 
tude of —4 is used. The astrometric case under the 
same conditions requires an effective temperature of 
1200°. 

There is thus a contradiction between the astrometric 
and the spectroscopic information at two points. (1) 
the effective temperature is too low; (2) the absolute 


WOW (C18, IP Te 18) I 
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JD 2435600 
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Mid-eclipse 


3rd 4th 


Frc. 11. Color curve from Larsson-Leander’s observations. Times of contact are given below the curve. Slopes for the body eclipse are indicated. 
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magnitudes disagree. If the temperatures are to be 
brought into agreement without changing the parallax 
very drastically, then the radius of the M star must 
be decreased. In order to decrease the size of the M 
star and still satisfy the eclipsing data the mass ratio 
must increase. This is not at all impossible. For ex- 
ample, if the mass ratio were unity the M star would 
have a radius of 500 Ro, an effective temperature of 
about 2100°, and each star would have a mass of 
10 Mo. The mass ratio may well be more nearly 2 than 1. 

The question of the absolute magnitude might 
perhaps be explained by an argument due to M. 
Schwarzschild. The M component of VV Cephei is 
known to be a magnetic variable (Babcock 1958) with 
a field strength varying between approximately +600 
to —600 gauss. The ratio of the magnetic energy 
density E, to the thermal energy density EH; may be 
written as 


2 


En wA|H 
E;  SxpkT 


where p is the density, & the gas constant, 7 the 
temperature, » the molecular weight, A the unit 
atomic weight, and |H| the absolute value of the 
field strength in gauss. At the surface of the sun the 
magnetic field is so small (a few gauss perhaps) that 
the ratio cited plays no role. However, in a sunspot 
where the magnetic field is large (a few thousand gauss) 
the ratio is nearly unity and the magnetic energy 
density can no longer be ignored. For VV Cephei 
p~ 10-8, T~3000°, |H|~600, and p= or 1. If p=3, 
the value of the ratio is approximately 3. If w=1 and 
T is decreased slightly, then the ratio is slightly 
greater than 10. Thus the effects of the magnetic field 
cannot be ignored when discussing the spectrum or 
spectral classification of this star. The effect would 
be one of lifting the atmosphere and thereby sharpening 
the absorption lines without appreciably increasing the 
size of the star. 

In the previous analysis the spectroscopic value for 
periastron passage has been used. The latest data still 
prefer 1942 and though this differs considerably from 
the value preferred by the astrometric material the 
disagreement might be found in the first four years of 
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the astrometric series. During these years there is an 
abnormally large scatter in the values of & when 
compared with the later years. This situation will be 
improved as the orbit is covered for the second time. — 

The absolute magnitude of the M star and its size 
depend strongly on the parallax which depends upon 
the linear value of a given by the spectroscopy and the 
angular value of a given by the astrometry, which is 
in turn affected bythe light of the B star. Within 
limits these-values will not change except possibly 
toward a slightly larger value of the parallax. The 
real unknown in the solution is the mass ratio, which 
must always be assumed. If a reasonable estimate of 
the mass ratio can be made the dimensions become 
more rigidly defined. This estimate can be made using 
a beam interferometer. | 

On the assumption of a mass ratio of unity the 
maximum separation between the two stars will occur 
around 1972 and will be on the order of 0"08. Such a 
measurement would require a beam of the order of 
two meters and should be made around 4200 A, where 
the two sources are equally bright. A larger angular 
separation would indicate that the B star is the less 
massive of the two stars, where as a smaller or un- 
measurable separation would indicate that the B star 
is the more massive. 

The geometry of the next eclipse can be determined 
more readily if the intrinsic variations of the stars are 
better known. The M star certainly possesses two 
variations, one with a 13.7-year period and the other 
of slightly less than a year. The latter of these may be 
the resultant of two minor variations; this can be 
simply tested by seeing if the L wave remains linear 
with the wavelength or not. Little is known of the B 
star’s variations except that they appear to be quite 
erratic. If this is the case the accuracy of the solution 
will always be limited. Photoelectric observations now 
in progress in the ultraviolet light, where the light of 
the B star dominates, should show whether the vari- 
ations of the B star’s light is periodic or not. At the 
same time observations in the near infrared follow the 
M star’s variations essentially free of the B star’s light. 

Because the diameter of the M star is rather large, 
an analysis of the spectroscopy around periastron 
passage (1942 and 1962) will be of considerable interest. 
Even if the M star were as small as 300 Ro, its boundary 
would be rather close to the Lagrangian inner-contact 
surface at periastron passage and border upon an 
unstable condition. Wood (1957) has considered the 
question of instability for eclipsing binaries based 
upon the well-known cases of period changes and 
erratic light-curve behavior. In the case of VV Cephei 
it is doubtful whether a change in period can be 
detected in the next few centuries. Also, the light 
curve is so erratic that it is difficult to associate general 
instability with it. P. Merrill has stated that the shell 
about the B star is constantly collapsing upon the B 


| 
star. This material must come from the M star’s 
atmosphere, which indicates that an unstable condition 
must be present. The collapsing shell may be the cause 
of the light variations associated with the B star. 
Recent spectroscopic observations at Lowell Observa- 
tory confirm that this is not a “normal’’ shell star. 
The Paschen-y line seems to be present at 10940 A, but 
the He line at 10830 A is not. In P Cyg, for example, 
both lines are present, the He line being the stronger. 
The Paschen-y identification is very weak due to the 
fact that the M-star’s continuum is very strong there. 
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In the course of a more extensive investigation of late-type stars in the vicinity of the north galactic pole, 
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a number of relatively cool stars have been found. The list of M and carbon stars prepared for the present 

discussion contains the magnitude and spectral class for each star and is probably complete to a photographic f 
magnitude of 13.0 in an area of about 400 square degrees. A limited statistical study indicates that the ratio if 
of dwarfs to giant stars to this limiting magnitude in the area covered, is about one to three. Nineteen 

dwarf M stars were found. The space density of these dwarfs was found to be about 39 stars per 1000 cubic 

parsecs as against 36.3 per 1000 cubic parsecs for the known stars within five parsecs of the sun. 


INTRODUCTION 


CATALOGUE is given which contains the spectral 

class, luminosity class, and magnitude of 80 
late-type stars of photographic magnitude 13.0 and 
brighter in a region of about 400 square degrees 
surrounding the north pole of the galaxy. The region 
extends from 11530™ to 13" in right ascension and from 
+25° to +50° in declination, and is described by 
Slettebak and Stock (1959). 

The catalogue includes 75 M-type stars, five carbon 
stars and no §S stars. Since the entire region was 
surveyed uniformly, this list may be considered to be 
complete to an approximately uniform limiting magni- 
tude. An identification chart showing all of the stars in 
the catalogue is soon to be published in a more extensive 
survey of late-type stars in this region. For the time 
being the coordinates given in Table I should be 
sufficient for identification. 


DATA Richt 


The spectrographic plates used for this survey were 
taken with the 80 cm—120 cm Hamburg Schmidt tele- 
scope and the 4° objective prism, and have been 
described by Slettebak and Stock (1959). A portion of 
the region was also photographed with Ila-O and 
103a-O plates of exposure times 3 minutes and 24 
minutes, respectively, taken with the Case Schmidt 
telescope and the new ultraviolet transmitting objective 
prism, which gives a dispersion of 580 A/mm at Hy, 
a dispersion almost identical to that of the Hamburg 
equipment. 


THE CATALOGUE 


The stars which have been classified in this investi- 
gation are listed in Table I. The first three columns 
give a catalogue number, and the BD and HD numbers 
for the stars common to those catalogues. All of the 
stars in Table I which are also in the AGK2 catalogue 
have 1950 positions as given in that catalogue. The 
remaining BD stars have positions obtained by apply- 
ing precession to the 1855 positions. All remaining 
stars have coordinates estimated from their positions 


on the plates and are probably accurate to about three 
minutes of arc. } 

The photographic magnitudes listed in Table I, 
column 6, were taken from the AGK2 catalogue when 
available except for two stars (No. 65 and 70), which 
have a known blue magnitude in the UBV system as 
given by Stock and Wehlau (1956). All remaining 
magnitudes are indicated by a colon and were deter- 
mined by flyspanking charts of the Lick Aélas. Photo- 
electric sequences from the list of Stock and Wehlau 
were used wherever possible. Otherwise, sequences 
were taken from the AGK2 catalogue. 


. 


REMARKS 


Each plate was examined for stars of spectral clas: 
G5 or later, and all stars selected were classified twice 
on the Yerkes system. The two independent classifica: 
tions were averaged. The stars which were classified as 
M stars were reclassified by direct comparison of the 
spectral images with those of standard stars of the 
Yerkes system of classification on plates taken witk 
the Case Schmidt telescope. 

Two luminosity criteria were used to separate the M 
stars into giants and dwarfs. These are the intensity o: 
the \4227 line of neutral calcium and the intensity o! 
the continuum between this line and the G banc 
relative to the continuum on the short-wave side of the 
4227 line, criteria 1 and 2 of Nassau and Keenar 
(1946). The internal consistency between the twe 
independent classifications using these criteria wa: 
very high. 

Nine stars, all classified here as dwarfs, are found i 
the lists of M dwarfs of Vyssotsky et al. (1943, 1946 
1952, and 1956) and two of these are also in the Propei 
Motion Catalogue of Luyten (1955). In addition, eight 
known variables and two additional stars listed a: 
giants by Adams ef al. (1935) are classified as giant: 
or emission-line stars except for two of the variables 
which are carbon stars. The spectral classes given ir 
column 8 are those of the reference listed in the same 
column. 

Nineteen of the 75 M stars in Table I are classifiec 
as dwarf stars, and all of the other 56 are assumed t¢ 
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No. BD HD (1950) (1950) 
1 1153479 48°590/ 
2 B30) nea 
3 42°2233 101366 37.4 42 18 
4 101605 SSutn 88047 
5 451955 101585 39.0 44 28 
6 36 2216 102159 43.0 36 10 
7 33 2151 Asda 33) 13 
8 44 2126 44.6 43 45 
9 43 2149 44.6 42 39 

10 49 2086 45.6 48 48 
11 36 2219 48.4 3532 
12 Sie ou 41102 
13 Giles a il 
14 37 2228 SONO MNS 25 
15 29 2228 Beso 02 
16 37 2230 103500 Dot — SW 
17 45 1977 103796 SAT AS) 17 
18 55.8 38 06 
19 40 2491 59.7, 39°50 
20 30 2217 104710 1201.0 2957 
21 29 2261 07.9 28 36 
22 08.6 41 19 
23 09.4 4959 
24 09.6 40 14 
25 OD. AE ay 
26 OG SO Sy 
27 35 2325 OS 
28 10.9) 31 58 
29 27 2104 lil Dyp WB 

30 49 2126 12.9 4901 

31 44 2175 16-2) 4212 

32 29 2271 13530 29105 

33 42 2287 IB AT 

34 28 2097 = 106814 TATA S101 
35 1508) 146) 53 

36 29 2279 Alans r40 

37 42 2296 1ORAD ee 42) 25, 

38 28 2110 20.1 «27 53 

39 50 1915 20.4 5019 
40 20.7 44 35 
41 Milaey VES) 
42 41 2292 107905 21.4 4100 
43 46 1784 DD 45 52 
44 47 1958 22.7 4659 
45 BB BESO) 
46 DE 37738 
47 25.4 4509 
48 30 2279 26.1 29 52 
49 42 2310 DSL LOO 
50 42 2312 DT A216 
51 28.0 31 47 
52 44 2199 29.6 43 46 
53 46 1792 29.8 45 46 
54 32) 36829 
55 32.8" (34.35 
56 33.7 4200 
57 27 2152 Saran 721 
58 33 2256 34501 3253 
59 36.8 2454 
60 38 2355 30.2) 37 30 
61 ANal ero 22 
62 40.1 43 19 
63 42 2334 110687 DDS EGR 
64 461817 110914 42.8 45 43 
65 27 2167 = 110964 HE BEDE 
66 48 2055 111129 44.3 47 39 
67 111223 45.6 38 39 
68 48.9 40 38 
69 47 1998 112082 51.0 4656 
70 29 2336 = 112172 51.9 28 43 
71 47 2003. ~—-: 112264 SO An 28 
72 36 2321 53y8) 3537 
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Remarks 


RU UMa 
Adams gM3 
TV UMa 


Vyssotsky No. 135 M2 


Vyssotsky No. 136 MO 


Adams gM5 
Vyssotsky No. 636 MO 


Vyssotsky No. 640 MO 


Vyssotsky No, 644 M2=LFT 896 
Vyssotsky No. 649 MO=LFT 898 
Vyssotsky No. 134 MOp 


Vyssotsky No. 654 MO 


T CVn 


Y CVn 
Stock and Wehlau No. 22 M3.5 


U CVn 


Stock and Wehlau No. 157 M4 
TU CVn 
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a 6 
No. BD HD (1950) (1950) Mptg Sp Remarks 
73 12554™1 43°45’ 12758 R 
74 43°2283 54.5 43 19 125 M4 Ill 
75 36 2322 Bpati 35 30 11.9 MO V Vyssotsky No. 298 MO 
76 56.9 40 50 12.5 MO0:V 
77 45 2076 57.0 44 50 OA M3 Il 
78 38 2389 112869 57.0 38 05 10.8 R TT CVn 
79 44 2243 58.0 44 10 19s M5 Ill TW CVn 
80 58.9 45 55 LOS MS II 


be giants, since the few for which no luminosity class 
is given in Table I are all either known variable stars 
or have a spectral class of M7 or later, probably too 
late a class for a dwarf star. No supergiants were 
detected although it is difficult to separate such stars 
from giants at the spectral dispersion employed here. 

Four of the five carbon stars were found to be of 
spectral type R and one of type N. Three of the four 
R stars appear to be previously unknown. All three 
are very faint, hence they lie at a considerable distance 
from the galactic plane. If the value —0.4 given by 
Sanford (1944) is taken for the absolute visual magni- 
tude for these stars, and +3.0, a mean of nine carbon 
stars listed by Johnson (1955), for the B-V color, the 
distances from the galactic plane are about 1000 
parsecs for stars No. 54, 61, and 73. 


STATISTICAL DISCUSSION 


The statistical results of this section are provisional 
in nature as they are based on magnitudes of limited 
accuracy. However, care has been taken to oo for 
scale and internal systematic errors. 

Shown in Fig. 1 are cumulative Wolf PEN for 
the giants and for the dwarfs listed in Table I. The 
abscissa for Fig. 1 gives the photographic magnitude, 
and the ordinate the number of stars brighter than that 
magnitude. The solid and open dots represent the 
total numbers of giants and dwarfs, respectively, to 
the indicated limiting apparent photographic magni- 
tude, and have been plotted for each half-magnitude 
increment. The theoretical values of the cumulative 
Wolf diagram are shown by the curved lines in Fig. 1. 
These curves have been computed on the assumption 
of a uniform space density of both giants and dwarfs 
and no space absorption. The curves have been fitted 
arbitrarily to the observed points at magnitude 7.5 for 
the giants, and 11.0 for the dwarfs, since these are the 
brightest limiting half-magnitudes which include at 
least one star. 

Only in the case of the dwarfs does the curve continue 
to show a good fit up to a photographic magnitude of 
13.0, the limit for this survey. This close agreement 
suggests that the limiting magnitude of 13 for this 
investigation is approximately correct. The space 
density of the giants according to Fig. 1 appears to 


decrease sharply for magnitudes fainter than 9.0, 
indicating a drop in the number of giants with in- 
creasing distance from the galactic plane. Those gian 
stars fainter than magnitude 9.0 belong primarily to 
the galactic halo. Although no definite mean absolute 
magnitude can be assigned to these stars, a rough 
estimate of —0.4, the value given by Keenan and 
Morgan (1951) for the mean visual absolute magnitude, 
and a color index of +1.6 (Allen 1955) yield a distance 
of about 400 parsecs, a reasonable value for th 
beginning of a transition between the edge of the disk 
and the start of the galactic halo. There is no coke 
difference in the concentration toward the galactic. 
plane between the early and the late M-giant stars. 
Thus, although many of the late M-giants may be 
intrinsic variables with magnitudes of questionable 
value, their inclusion in this statistical discussion does 
not appreciably alter the results. 


e GIANTS 
© DWARFS 


FREQUENCY 


Fic. 1. Cumulative Wolf diagrams for giants and dwarfs. 
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An estimate of the space density of the 19 M-dwarf 
tars shows a value in general agreement with that of 
he solar neighborhood. According to Gliese (1957) 
here are 19 single stars or double systems within five 
yarsecs of the sun which have spectral classes of MO to 
V4. If a mean absolute visual magnitude of +9.8 and 
. mean color index of +1.4 are assumed for the 19 
lwarf stars listed in Table I [the values are those for 
Mi2V stars as given by Allen (1955) ] then the resulting 
nodulus for the limiting distance for dwarf stars in 
his investigation is 1.8, and the corresponding limiting 
listance is 23 parsecs. 

The volumes are 524 and 487 cubic parsecs, respec- 
ively, for the space within five parsecs of the sun and 
he volume covered in this survey. The density of M 
lwarfs based on these figures is 36.3 and 39.0 M-dwarfs 
xer 1000 cubic parsecs, respectively (for the two 
volumes). The agreement between these figures in- 
licates that the number of dwarfs found in this search 
s about as expected. 
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FURTHER REPORTS OF OBSERVATORIES 


Lockheed Solar Observatory, Lockheed Aircraft 
Corporation, Burbank, California. 


A study of the solar atmosphere by means of high- 
speed Hq cinematography was initiated by the Lock- 
heed Aircraft Corporation, California Division, in the 
fall of 1958. This report encompasses the 18-month 
period commencing with inception of routine obser- 
vations January 1, 1959, closing June 30, 1960. 

Personnel. The Solar Physics Project has been 
under the directorship of G. E. Moreton, previously 
Principal Investigator, Convair Radio Astronomy 
Project. H. E. Ramsey, formerly Sacramento Peak 
Observatory Chief Observer, joined the project Sep- 
tember 1959. In addition, the research group consists 
of J. W. Harvey, N. W. Christie, and G. F. Ander- 
son. G. A. Carroll is in charge of instrumentation 
design and construction. 

Some recent visitors have been Drs. R. Grant 
Athay, J. W. Warwick, Richard Hansen, and Nor- 
man Macdonald, of the High Altitude Observatory ; 
Dr. James W. Dungey, Atomic Weapons Research 
Establishment, England; Professor Robert Leighton, 
California Institute of Technology; Dr. Eugene N. 
Parker, University of Chicago; Dr. I. G. Poppoff, 
Stanford Research Institute; and Professor John R, 
Winckler, University of Minnesota. 

Facility. The observing site is located at an eleva- 
tion of 1500 feet, in the Hollywood Hills, Los An- 
geles. The well-known atmospheric inversion, re- 
sulting in local “smog,” has proved beneficial: in- 
hibited convective air flow apparently contributes to 
prolonged periods of good seeing. Long periods of 
cloud-free sky allow excellent observational con- 
tinuity: in spite of equipment-down time during the 
initial weeks, observations have been conducted 85% 
of the days in this reporting interval. The. facility 
consists of a heliograph incorporating a three-inch 
f/33 objective lens, secondary reimaging optics, Hallé 
0.5 A bandpass Hq birefringent filter, and Acme pin- 
registration camera. The system is photoelectrically 
guided, and mounted on the original Climax corona- 
graph loaned by HAO for this work, and modified to 
serve as a heliograph by Lockheed. The 18-mm solar 
image is photographed each 10 sec on 200-foot 
lengths of Eastman 4-E 35-mm film. A clock is 
imaged adjacent to the sun; photometric standards 
are impressed on each roll. 

Observing Program. Six frame per minute cine- 
matographic observations are conducted daily from 
about 1600 hr UT to sunset. The 35-mm negatives 
are scanned daily with an automatic film reader, and 
occurrences of flares and associated phenomena re- 
ported in accordance with IGY-IGC procedures. 
Daily flare reports are made by telegraph to the 


Boulder World Data Center ; other World Data Ce 
ters receive tabulated data monthly. Concomitant v 
ual observations are conducted daily in the aft 
noons and continuously on weekends ; occurrences 
major flares are reported by telephone to the | 
Belvoir World Warning Agency; such observatic 
have resulted in five world-wide geophysical aler 
During summer 1959°a continuous visual patrol ; 
sulted in successful rocket firings from Pt. Argue 
by the Naval Research Laboratory team measuri 
flare-enhanced x rays. 

In order to handle the flare data as routinely 
possible, extensive measurements for each flare < 
recorded on IBM punch cards and data compilati 
is accomplished by IBM tabulation. Monthly flz 
patrol results appear regularly in those publicatic 
concerned with solar data. 

Beginning January 1960 cinematography was < 
complished with the filter band pass displaced < 
proximately 0.5 A to the blue side of Ha. Futt 
observational research will be conducted with thi 
sequentially programed exposures each 10 sec at J 
and at wavelengths displaced to the blue and red. 

Since May 1960 the IGY-IGC activities ha 
received partial support from the National Scier 
Foundation. 

Research Summary. An immediate consequence 
the high time-resolution cinematography, good s« 
ing, clear sky, and modern data-processing tec 
niques has been the detection of a large numl 
of small-scale short-duration chromospheric even 
During the 2932 observing hours, a total of 40 
flares have been detected and recorded on pun 
cards ; 94% of these events are subflares. The month 
flare frequency averaged about 1.6 flares per hor 
Of the flares recorded by Doppler-shifted radiati 
in 1960, 50% have an associated surge or filame 
activation. 

Studies of the influence of flares on filaments 
Moreton have resulted in detection of seven eve 
where dark filaments were unambiguously activat 
by flares separated by distances up to 1Ro. The fi 
ments undergo “abrupt” activation; initial opac 
changes can be determined to +10 sec. The respon 
ble flares appear to have a distinct “explosive” phe 
suggestive of hydrodynamic effects independent 
other observed phenomena. If a disturbance leat 
the flare at the explosive phase and travels the sho 
est route to the filaments, then the inferred veloc 
of propagation for the disurbances is ~1000 km/s 
The velocity suggests close association with hi 
velocity disturbances assumed to cause type-II rac 
emission and associated with geomagnetic storms. 

Further investigation of this new class of evet 
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gan January 1960 by operating the heliograph at 
»—0.5 A, enhancing Doppler-shifted radiation from 
aments and moving features. This has resulted in 
nfirmation of inferred 1000 km/sec disturbances 
9m flares: three events have been photographed 
nere direct measurement of faint propagating dis- 
rbances confirm the velocities obtained from sub- 
quent filament activation. 
This report interval culminated with a singularly 
mificant observation: on June 25, 1960 an exten- 
ve disturbance was photographed propagating from 
flare at ~2500 km/sec. The disturbance was visi- 
e over 400 000 km from the flare. 
H. E. Ramsey has been analyzing film records, 
udying loop prominences seen in absorption on the 
sk, and the occurrences and properties of preflare 
ament activations. J. W. Harvey has developed a 
1otographic technique which may make possible 
sy and accurate measurement of high-velocity flare 
cpansions. N. W. Christie has investigated data- 
indling techniques and is collaborating with Harvey 
1 statistical analysis of subflares, surges, and radio 
vise. 
A discussion of flare-activated filaments has been 
ympleted by R. Grant Athay, of HAO, and More- 
nn, and has been submitted for publication. A study 
is been completed by Moreton in collaboration with 
ugene Greenstadt, of Space Technology Labora- 
ries, relating Satellite Pioneer V magnetometer 
ata to solar activity; the results are being readied 
yr publication. 

G. E. Moreton, in charge 


[ational Aeronautics and Space Administration, 
Washington 25, D. C. | 


OBSERVATIONAL PROGRAM 


Vanguard III and Pioneer V were launched within 
lis period. The Lyman-alpha detectors on Vanguard 
[I saturated in the Van Allen Belt but the micro- 
1eteorite experiment confirmed the results from Ex- 
lorer I within a factor of 2. From these observations 
1e dust density in space near the earth is estimated 
» be approximately 10°°° g/cm*. 

Pioneer V contained experiments for measuring 
1e interplanetary magnetic field and ion densities in 
iterplanetary space. An intense solar storm occurred 
n April 1, during the active life of this probe. It 
vas found that the magnetic field disturbance follow- 
ig a solar flare is not localized in the earth’s mag- 
etic field as had previously been assumed. Some 
f the interesting results are 

(1) An interplanetary magnetic field exists with 
steady magnitude of slightly more than one gamma. 
‘his field fluctuates in intensity up to ten gammas 
uring periods of solar flare activity. 

(2) The planar angle of the interplanetary mag- 
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netic field forms a large angle with the plane of the 
ecliptic, perhaps as much as 90 deg. 

(3) The exospheric ring current about 2500 miles 
in diameter and at a distance of about 40 000 miles 
from the earth causes a westward moving current 
estimated at five million amperes. 

(4) The geophysical magnetic field extends at 
times 65 000 miles from the earth. 

(5) The sudden decrease in galactic cosmic rays 
associated with large solar flares does not depend on 
the presence of the earth’s magnetic field. 

Goddard Space Flight Center launched three night- 
sky rockets in the spring of 1960 for stellar and 
nebular photometry in the ultraviolet. The results 
have not yet been reduced but the flights were suc- 
cessful and the program will be continued. Develop- 
ment of detectors for the ultraviolet has been under- 
taken and will be continued and expanded during the 
coming year. 

A program for the scientific study of the moon, 
planets, and interplanetary space by means of space 
probes was formulated and instrumentation develop- 
ment has been started for the early flights in this 
program. The first lunar flights in this program will 
carry a vidicon camera and gamma-ray spectrometer 
on the spacecraft and a single-axis seismometer and 
decelerometer on a capsule to be impacted on the 
moon, The first planetary probe will contain an ul- 
traviolet spectrograph and radiometer. All flights will 
carry an interplanetary package consisting of radia- 
tion experiments in various energy ranges and a 
magnetometer. 

Instrumentation is partially complete for a satellite 
to study the sun in the wavelength region below 
2000 A. Experiments for the Orbiting Solar Observ- 
atory Satellite will be furnished by Ames Research 
Center, Goddard Space Flight Center, University of 
California, University of Colorado, University of 
Minnesota, and the University of Rochester. 

The George C. Marshall Space Flight Center was 
preparing a satellite to observe 100 Mev gamma radi- 
ation from space. The Massachusetts Institute of 
Technology is responsible for the scientific instru- 
mentation and the Marshall Space Flight Center for 
the satellite structure, power supplies, telemetry sys- 
tem, and auxilliary instrumentation. 

The Marshall Space Flight Center is studying the 
motions of a satellite about its center of mass due to 
interactions with the space environment and_ is 
analyzing the observed motions of Explorers IV and 
V. They are also studying the applicability of image 
intensifier and image orthicon techniques. Negotia- 
tions are underway for a test using minor planets to 
evaluate these techniques. 

The Ames Research Center carried out an exten- 
sive study of a satellite to be used for astronomical 
observations with a telescope up to 36 inches in di- 
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ameter. Practical problems associated with control, 
guidance, and data transmission of such a satellite 
were studied, particularly under the direction of R. 
M. Crane. Removal of the primary component of the 
earth’s gravitational field may permit particularly 
delicate control of the instruments. Among disturb- 
ing factors remaining are 


(1) The gradient of the earth’s gravitational field. 
(2) Solar radiation pressure. 

(3) Interaction with the earth’s magnetic field. 
(4) Aberration of light. 


THEORETICAL AND LABORATORY STUDIES 


Jet Propulsion Laboratory. Y. Hiroshige and F. 
Yagi worked on the theory of the interaction between 
a dipole magnetic field and a conducting fluid issuing 
from the center of the dipole field. This has obvious 
application to the sun. H. Lass studied dynamical 
problems associated with a satellite of the moon. He 
concluded that little knowledge of the moon’s in- 
ternal structure would be gained from one. He also 
formulated for a steady-state universe a special rela- 
tivity cosmology which satisfies all observations. T. 
Kreiter studied lunar crater distribution extensively. 
R. Carr has been investigating the use of seismo- 
logical techniques for determining the structure of the 


moon, particularly viewed as a layered medium con-, 


sisting of concentric spherical shells. C. Barth and 
A. Hildebrandt measured interaction processes be- 
tween atoms, molecules, ions, and electrons that occur 
in the upper atmosphere of the planets. Electron para- 
magnetic resonance spectrometry and optical spec- 
troscopy were used. 

The Ames Research Center. The Ames Research 
Center has long specialized in studies of high-speed 
dynamic phenomena. Under the direction of H. Julian 
Allen, studies of the interaction of meteors with the 
earth’s atmosphere were made. Comparison..of the 
results with several well-recorded meteor falls indi- 
cated that the phenomenon of deceleration and abla- 
tion of material in iron meteorites can be accounted 
for. 

Goddard Space Flight Center. Under the direction 
of the Goddard Space Flight Center some minitrack 
stations were used to track the stronger radio stars 
(Cygnus and Cassiopeia) and the sun. The stations 
used are near the same longitude and 108-Mc radio 
bursts could thus be recorded simultaneously at sev- 
eral locations. Combining these observations nearly 
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eliminates atmospheric effects and positions can be 
obtained which are accurate to two minutes of are. 
Those solar bursts which can be observed with the 
minitrack systems are time-correlated with solar 
H-alpha flares. S. Huang studied the possibility of 
habitable planets in other solar systems and pre- 
dicted. regions in which a fourth body can be stable 
in the earth-moon or similar systems. He also started 
computations of stellar models by a perturbatio 
technique. B. Donn showed that ion clusters of the 
type hypothesized by V. I. Krassovsky cannot occur 
in interstellar space. In addition, he applied various 
theoretical and observational results from other fields 
to the origin and properties of interstellar grains and 
studied properties of an icy-conglomerate comet 
nucleus. 

In the Theoretical Division of the Goddard Spas 
Flight Center, J. O’Keefe, Ann Eckels Bailie, and 
K. Squires completed their analysis of the Vanguard 
orbit for the determination of the earth’s gravita- 
tional field. P. Musen developed the theory for luni- 
solar effects on an artificial satellite and more re- 
cently, the theory of radiation pressure effects from 
the sun. The effect on the Vanguard orbit was found 
to be measurable ; that on Echo is part of the present 
Echo observational program. As a result of B. Strom- 
gren’s visit, a_project of astrophysical calculations in 
the field of stellar interiors was initiated. This in- 
volved much of the junior staff. C. Hayashi and 
R. Cameron studied the evolution of giant stars with 
special reference to the double cluster h and x Persei. 

R. Jastrow and W. Cameron in conjunction with 
Professor Urey of the University of California in- 
itiated an investigation of the early stages of forma- 
tion of the planetary bodies. G. J. F. MacDonald 
continued his study of resonances in the earth result- 
ing from major earthquakes. G. J. F. MacDonald, 
R. Jastrow, and P. Musen have been redeveloping the 
lunar theory to improve the values of the moments of 
inertia of the moon. J. O’Keefe and P. Lowman ‘in- 
vestigated tektites. G. J. F. MacDonald is studying 
free oscillations in the earth’s atmosphere as a guide 
to large-scale effects produced by thermal sources. 
I. Harris and R. Jastrow in collaboration with J. 
Chamberlain of Yerkes Observatory have been cal- 
culating the thermal equilibrium of the upper atmos- 
phere and the effect of temporal variation in solar 
radiation. 

Nancy G. Roman, Chief, 
Astronomy and Astrophysics 
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ABSTRACT 


The Luminosity Functions of Galactic Star 
vlusters.* SIDNEY VAN DEN BercH AND Davip 
HER, David Dunlap Observatory—The luminosity 
unctions of 20 galactic clusters down to m,,~20 
lave been determined by means of star counts on 
70 plates obtained with the Palomar 48-inch 
schmidt. It is found that striking differences exist 
etween the main-sequence luminosity functions of 
ndividual clusters. Also it appears that the faint 
nds of the luminosity functions of galactic clusters 
iffer systematically from the van Rhijn-Luyten 
uminosity function for field stars in the vicinity of 
he sun in the sense that (with one exception) all 
he clusters which were investigated had luminosity 


functions which either decrease or remain constant 
below M,,=+5. The differences between individual 
clusters and the differences between the luminosity 
functions of clusters on the one hand and field stars 
on the other show that the luminosity function of 
star creation is not unique. This result is taken to 
indicate that the luminosity function with which stars 
are created probably depends on the physical condi- 
tions prevailing in the region of star formation. 

* This abstract of a paper presented at the 106th Meeting 
of the American Astronomical Society was omitted, through 
an editorial oversight, from the abstracts of the meeting 


published in the November 1960 No. 1284 issue of the 
Astronomical Journal. 


NOTICE 


Graduate Laboratory Development Program 


HE National Science Foundation announces 
that March 1, 1961 is the next closing date for 
eceipt of proposals in the Graduate Laboratory De- 
elopment Program. Proposals received after March 
will be reviewed following the next closing date, 
september 1, 1961. This program requires at least 
0% participation by the institution with funds 
lerived from non-Federal sources. ~ 
Purpose of the grants is to aid institutions of 
iigher education in modernizing, renovating, or ex- 


panding graduate-level basic research laboratories 
used by staff members and graduate students. Only 
those departments having an on-going graduate train- 
ing program leading to the doctoral degree in science 
at the time of proposal submission are eligible at 
present. 

Proposals, as well as requests for additional infor- 
mation, should be addressed to: Office of Institu- 
tional Programs, National Science Foundation, 


Washington 25, D. C. 


ERRATUM 


Astron. J. 65, 354. 


Left column, line —8: for 900 read 800. 

Left column, line —5: for IV read III-IV. 

Right column: delete last 11 lines; substitute “The 
maximum in the space velocity frequency distribu- 
tion for the thirty weak CN stars occurs at 
approximately 50 km/sec, with only 20% of the 
velocities exceeding 65 km/sec.” 
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CQ And. 65: 381—1960 

U Adil. 65: 473—1960 

VY Aql. 65: 381—1960 

FF Aql. 65: 473—1960 

HH, V345 Aql. 65: 381—1960 

AW, AY Aur. 65: 381—1960 

e Aur. 65: 497—1960(A) 

¢ Aur. 65: 497—1960(A) 

SU, SW, WY, YZ Cam. 65: 381—1960 

R CMa. 65: 326—1960 

DL CMa. 65: 381—1960 

Y CVn. 65: 503—1960(A) 

SU Cas. 65: 473—1960 

AO Cas. 65: 127—1960 

DI, IW, KO Cas. 65: 381—1960 

VV Cep. 65: 628—1960 

AL, CU Cep. 65: 381—1960 

TZ Cr A. 65: 23—1960 

X, SU Cyg. 65: 473—1960 

XY Cyg. 65: 381—1960 

AG, AT, BB, CZ, DI Cyg. 65: 381—1960 

DT Cyg. 65: 473—1960 

FG, GQ, V429, V686 Cyg. 65: 381—1960 

RW, RY, SS, SZ, WX, AG, AN, BB Del. 65: 381—1960 

SY Dra. 65: 381—1960 

SY, UV Eri. 65: 381—1960 

AS Eri. 65: 139—1960 

BK, DF Her. 65: 381—1960 

FQ, FT Hya. 65: 381—1960 

XY Leo. 65: 3741960 

B Lyr. 65: 84—1960 

BC, BD, IK Mon. 65: 381—1960 

U Oph. 65: 55—1960 

Y Oph. 65: 473—1960 

Al, V374, V389 Oph. 65: 381—1960 

V839 Oph. 65: 79—1960 

U Peg. 65: 88—1960 

TU, TV, AP, EG Peg. 65: 381—i960 

GG Per. 65: 381—1960 

S, Y Ser. 65: 473—1960 

AH Ser. 65: 381—1960 

GO1GXx Hk HP. Vs Exe TES IM AV NOLO: 
LW, V945 Ser. 65: 101—1960 

HV 9461, 9486; Maria Mitchell a,b,c,d, in Sgr. 65: 101— 
1960 

V, SV Sct. 65: 381—1960 

WW, BG Ser. 65: 381—1960 

AH Vir. 65: 358—1960 

T Vul. 65: 473—1960 

WX Vul. 65: 381—1960 

SN (1954) NGC 5668, SN (1956) NGC 3992, SN (1957) 
NGC 4374. 65: 54-1960(A) 

Small Mag. Cld. 
HV 837, 840, 843, 844, 847, 848, 850, 856, 1695, 1744, 
1768, 1785, 1787, 1793, 1809, 1825, 1827, 1836, 1850, 1855, 
1858, 1868, 1873, 1877, 1884, 1891, 1892, 1898, 1903, 1905, 
1923, 1925, 1929, 1933, 1934, 1945, 1950, 1954, 1966, 1967, 
1974, 1978, 1981, 1987, 1994, 2000, 2002, 2014, 2017, 2019, 
2021, 2022, 2027, 2035, 2046, 2049, 2060, 2063, 2064, 11182, 
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